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In the beginning of the 1970's, Wilson developed the concept of a fully non-perturbative renor- 
malization group transformation. Applied to the Kondo problem, this numerical renormalization 
group method (NRG) gave for the first time the full crossover from the high-temperature phase 
of a free spin to the low-temperature phase of a completely screened spin. The NRG has been 
later generalized to a variety of quantum impurity problems. The purpose of this review is to give 
a brief introduction to the NRG method including some guidelines of how to calculate physical 
quantities, and to survey the development of the NRG method and its various applications over 
the last 30 years. These applications include variants of the original Kondo problem such as the 
non-Fermi liquid behavior in the two-channel Kondo model, dissipative quantum systems such as 
the spin-boson model, and lattice systems in the framework of the dynamical mean field theory. 
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I. INTRODUCTION 

The last decades saw a steadily increasing interest in a 
wide range of physical systems involving quantum impu- 
rities. The expression 'quantum impurity system' is used 
in a very general sense here, namely a small system (the 
impurity) with only a few degrees of freedom coupled to 
a large system (the environment or bath) with very many 
degrees of freedom, and where both subsystems have to 
be treated quantum mechanically. 

The use of the terminology 'impurity' is due to histori- 
cal reasons. In the Kondo problem, the small system is a 
magnetic impurity, such as an iron ion, interacting with 
the c onduction elect rons of a nonmagnetic metal such as 
gold (|Hewsonl . [l993h . Other realizations are for example 
artificial impurities such as quantum dots hosting only 
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a small number of electrons. Here, the environment is 
formed by the electrons in the leads. The term 'quantum 
impurity systems' can also be used for what are tradi- 
tionally called 'dissipative systems'. As an example, let 
the impurity correspond to a spin degree of freedom and 
the environment be built up by a bosonic bath; this de- 
scribes the well-known spin-boson model experimentally 
relevant, for example, for dissi pative two-level syste ms 
like tunneling centers in glasses ( Leggett et all . [l987h . 

Any theoretical method for the investigation of quan- 
tum impurity systems has to face a number of serious 
obstacles. First of all, because the environment typically 
consists of a (quasi-)continuum of quantum-mechanical 
degrees of freedom, one has to consider a wide range of 
energies - from a high-energy cut-off (which can be of 
the order of several eV) down to arbitrarily small excita- 
tion energies. On the other hand, because the impurity 
degrees of freedom usually form an interacting quantum- 
mechanical system, their coupling to a continuum of ex- 
citations with arbitrarily small energies can result in in- 
frared divergencies in perturbational treatments. A well- 
known example for this difficulty is the Kondo problem: 
Its physics is governed by an energy scale, the Kondo 
temperature 7k, which depends non- analytically on the 
spin-exchange coupling J between the impurity and the 
condu ction band of the host, lnTk cx 1/J (see iHewsonl 
(1993) :br a detailed description of the limitations of the 
perturbational approach for the Kondo model and the 
single-impurity Anderson model) . One is thus faced with 
the task to perform non-perturbative calculations for an 
interacting quantum-mechanical many-body system with 
a continuum of excitations covering a broad spectrum of 
energies. 

A very efficient way to treat systems with such a broad 
and continuous spectrum of energies is the renormaliza- 
tion group approach. It allows, in general, to go in a cer- 
tain sequence of renormalization group steps from high 
energies, such as the bandwidth, to low energies, such as 
the Kondo temperature. General introductio ns of renor- 
malization group c oncepts have bee n given in Goldcnfeld 



mahzation group c oncepts have bee n given m Hjpldenteldl 
(l!992l): lMal ll976l):ISalmhofeij (I1999D (se e also the origin al 
papers: IWilson and KogutI (|l974 ) and lWilsonl (|l975bl )). 
Here we focus on a specific implementation of the renor- 
malization group i dea: Wilson's n umerical renormaliza- 
tion group method ( Wilsonl . ll975at) . referred to as 'NRG' 
in the remainder of the review. This approach is dif- 
ferent from most renormalization group methods as it is 
non-perturbative in all system parameters; however, the 
price one has to pay is that the renormalization group 
steps have to be performed numerically. 

The general strategy of the NRG is the following (more 
details are given in Sec. [TTJ) . As specific example, let 
us consider the Kondo model which describes a mag- 
netic impurity with spin S coupled to the electrons of 
a conduction band, assumed to be non-interacting, via 
an interaction of the form JS ■ s, with s the spin of the 
electrons at the impurity site. The NRG starts with 
a logarithmic discretization of the conduction band in 



intervals [A~^ n+1 ^uj c , A""o; c ] and [-A""o; C) -A-(" +1 )w c ] 
(n = 0, 1, 2, . . .). We shall call A > 1 the NRG discretiza- 
tion parameter. After a sequence of transformations, the 
discretized model is mapped onto a semi-infinite chain 
with the impurity spin representing the first site of the 
chain. The Kondo model in the semi-infinite chain form 
is diagonalized iteratively, starting from the impurity site 
and successively adding degrees of freedom to the chain. 
Due to the logarithmic discretization, the hopping pa- 
rameters between neighboring sites fall off exponentially, 
i.e. going along the chain corresponds to accessing de- 
creasing energy scales in the calculation. 

In this way, Wilson achieved a non-perturbative de- 
scription of the full crossover from a free impurity spin at 
high temperatu res to a screened spin at low temperatures 



high temperatures to a screened spin at low temperatures 
(IWilsonl . TT975ah . thus solvin g the so-called K ondo prob- 
lem as discussed in detail in IHewsonl ( 19931 ). After this 
first application more than 30 years ago, the NRG has 
been successfully generalized and applied to a much wider 
range of quantum impurity problems. The first exten- 
sion was th e investiga t ion o f the single-impurity Ander- 
son model ([Andersonl Il96ll ). which extends the Kondo 
mode l by including charge fl uctuat ions at the impurity 
site. iKrishna-murthv et all (|l980al ibl) described all the 
technical details, the analysis of fixed points, and the 
calculation of static quantities for this model. 

In the following, the development of the NRG con- 
centrated on the analysis of more complicated impurity 
models, involving either more environment or impurity 
degrees of freedom. For example, in the two-channel 
Kondo model the impurity spin couples to two conduc- 
tion bands. This model, which has a non-Fermi liquid 
fixed point with associated non-Fermi liqui d properties, 
has be en first investigated with the NRG bv lCragg et all 
(1980). The numerical calculations for such a two- 
channel model are, however, much more cumbersome 
because the Hilbert space grows by a factor 16 in each 
step of the iterative diagonalizati on, instead of the fac - 
tor 4 in the singl e -chan nel case. iPang and Col (Il99lh 
and lAffleck et all (Il992ft later analyzed the stability of 
the non-Fermi liquid fixed point with respect to various 
perturbations such as channel anisotropy. 

The two-impurity Kondo model as paradigm for 
the competition of local Kondo screening and non- 
local inB£netic__correlations_was studied w ith NRG 



Jones and Varmal (Il987l ) 
Sakai and Shimizul i|l992allbl ): 



by 



iJones et all 



.Sakai et all 

Silva et al. 1 1996?) . Here, the focus was on the question, 



(1988); 
19901 ): 



if the two regimes are connected by a quantum-phase 
transition or rather by a smooth crossover. Later on, 
such studies were ex t ended to the two-channel situation, 
too (jlngersent et all , Il992l ). 

Originally, the NRG was used to determine thermo- 
dynamic properties of quantum impurity systems. The 
calculation of dynamic quantities with the NRG started 
with the T = absorption and photoemission spectra 
of th e x-ray Hamiltonian ( Oliveira and Wilkinsl . Il98ll 
Il985l ) , followed by the T = single-particle spectral func- 



3 



tion for the or bitally non-degenerate and degenerate An- 
derson model (Brit o and Frot I Costi and Hewsonl . 

1990L Il992ri IFrota and Oliveiral Il986t ISakai et all 
1989T ). The resulting spectral functions are obtained on 
all energy scales, with a resolution proportional to the 
frequency as discussed in detail in Sec. lIII.Fll Calculation 
of finite-temperature spectral functions A(u, T) are more 
problematic since all excitations can, in principle, con- 
tribute. Nevertheless, the NRG has been shown to give 
accurate results for A(u>, T), t oo, which also allows to cal- 
culate trans port p r operties (ICosti and Hewsonl . Il992ri 
ICosti et al. I. Il994at Isirzuki et all Il996l) . A subsequent 
development is the introduction of the concept of the 
reduced density matrix, which allows to calculate dy- 
namic quanti ties in equilibriu m in the presence of ex- 
ternal fields (jHofstetterl . l2000f ) . The calculation of non- 
equilibrium transi ent dynamics requires a multiple-shell 
NRG procedure (ICostil . Il997al) and has been accom- 



plished with the aid of a complete basis se t and the re- 
duced density matrix ( Anders and Schillerl . [2005) . The 
first applications of this approach show very promis- 
i ng results, b oth for fermionic and bosonic s ystem s 
(|Anders et all 120061 : lAnders and Schillerl . I2005L l2006h . 



Another recent generalization of the NRG approach is to 
quantum impurities couple d to a bosonic bath (bosonic 
NR G, see Bulla et al\ (l2005[) ; for ea rly related approaches 
see lEvangelou and Hewsonl ( 19821 )). The bosonic NRG 
has already been successfully applied to the sub-Ohmic 
spin-boson model which shows a line of quantum crit- 
i cal points separa ting localized and delocalized phases 
dBulla et aLl . l2003h . 

Additional motivation to further improve the NRG 
method came from the develop ment of the dynam- 
ical mean- field theory (DMFT) (jGeorges et al. ,199G; 
Met zner and Vollhardtl I 19891 ) in which a lattice model 
of correlated electrons, such as the Hubbard model, is 
mapped onto a single-impurity Anderson model with 
the impurity coupled to a bath whose structure has to 
be determined self-consistently. This requires the NRG 
to handle impurity models with an arbitrary density of 
states of the conduction electr ons and to calcula te di- 
rectly the impurity self-energy ( Bulla et~al\ . [1998) . The 
first applications of the NRG within the DMFT frame- 
work concentrated on the Mott transition in the Hub- 
bard model a nd accurate results could be obtained for 
both T = ([Bulla! . Il999[ ISakai and Kuramotol Il994h 
and finite temperatures ([Bulla et all l200lh . Within 
DMFT, the N RG has been applied to the periodic An- 
derson mod el (fPruschke et all |2000|) . the Kondo lat- 
tice model dCosti and Maniml. |2002|), m ulti-band Hub- 
bard models (|Pruschke and Bullal |2005[ ). the ferromag- 
netic Kondo lattice mode l with interactions in the band 
( Liebsch and Costi |2006[ ). and to lattice models with a 
couplin g to local phonon modes such as the Holstein 
model dMever et all 2002 ) and the Hubbard- Holstein 



model (|Koller et all [2004b) 



leads can give rise to Kondo-like features in the trans- 
port characteristics, has also led to a resurgence of in- 
terest in quantum impurity systems, both experimen- 
tally and theoretically. An important feature of quan- 
tum dot systems is their flexibility and a number of 
different set-ups have been realized so far, and inves- 
tigated theoretically by various methods including the 
NRG. Applications of t he NRG in t h is field include the 
standard Kondo effect dBorda et all. 120051 : ICostil l200ll 
iGerland et aZl.l2000Hlzumida et al.l.ll998l). coupled quan- 
tum d o ts ( Borda et all 20031 ; Cornaglia and Grempel , 



2005bl: IGalpin et all 12006 ri IHofstetter and Schoellei 



r 

120021 ; IHofstetter and Zarand. 12004), q uantum dots in 
a superconductor ( Choi et all 12004a), and quantum 
dots coupled to ferro magnetic leads ( Choi et all , l2004ri 
iMartinek et all 120031) . 

From this brief overview one can see that the range 
of applicability of the NRG has widened considerably 
since Wilson's original paper, covering physical phenom- 
ena such as the Mott transition, quantum dot physics, lo- 
cal criticality, dissipative quantum systems, etc. Further 
applications are still lying ahead and various optimiza- 
tions of the technique itself are still being developed - we 
shall come back to this point in the summary section. 

This paper is the first review of the NRG approach 
(since Wilson's original paper on the Kondo problem) 
which attempts to cover both the technical details and 
all the various applications. In this way, the reader 
should get an overview over the field, learn about the 
current status of the individual applications, and (hope- 
fully) come up with ideas for further calculations. This 
review can only be a start for a deeper understanding of 
the NRG. The following shorter reviews on sele cted top- 
ics are also helpful: section 4 in iHewso n (1993) contains 
a pedagogical int roduction to the NRG as applied to the 
Kondo problem, iGonzalez-Buxton and Ingersentl (Il998l) 
discus s the soft-gap Anderson and Kondo models. ICostil 
(1999) gives a general overview of the key concepts, in- 
cluding the applic ation to the anisotropic Kondo model, 
Bulla et al\ (|2005l ) present a detailed introduction to the 
bosonic NRG, and, finally, the two papers on the first 
calculations for the si n gle-impu rity Anderson model by 
iKrishna-murthv et all (1980a. b) are still valuable read- 



The observation that the coupling between electronic 
degrees of freedom in quantum dots and the surrounding 



ing for an overview of the method and the details of the 
analysis of fixed points. 

The review is organized as follows: in Sec. [IT] we start 
with an introduction to the basic concepts of the NRG 
approach. The single-impurity Anderson model serves 
as an example here, but the strategy applies to quan- 
tum impurity systems quite generally. At the end of this 
section, we discuss the flow of the many-particle eigenen- 
ergies and the appearance of fixed points in this flow. 
This analysis already gives important insights into the 
physics of a given model, but the calculation of physical 
quantities needs some extra care, as described in Sec. lIIII 
This section is divided into Sec. IIII.AI (thermodynamic 
and static quantities, such as entropy, specific heat and 
susceptibilities) and Sec. 1111.51 (dynamic quantities, both 
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in equilibrium and non-equilibrium). 

The following two sections deal with the various ap- 
plications of the NRG and we distinguish here between 
quantum impurity systems (Sec. IIV[) and lattice mod- 
els within DMFT (Sec. [VJ). Section ITVl covers most of 
the work using the NRG which has been published so 
far. We shall present results for systems which show 
the standard Kondo effect (Sec. HV.Ai this also includes 
most of the NRG-results on quantum dots), the two- 
channel Kondo problem (Sec. IIV.B[) . models displaying 
impurity quantum phase transitions (Sec. IIV.C[) . quan- 
tum impurity systems with orbital degrees of freedom 
(Sec. lIV.Dp . and, finally, impurities coupled to bosonic 
degrees of freedom (Sec. IIV.E[) . The section on lattice 
models within DMFT (Sec.|Vj covers calculations for the 
Hubbard model (Sec. IV. A|) . the periodic Anderson and 
Kondo lattice models fSec. lV.Bj) . and lattice models with 
coupling to phonons (Sec. IV. Cj) . 

In the summary we shall discuss open problems as well 
as possible directions for future developments of the NRG 
approach. 

Let us finish the introduction with a few remarks on 
the selection of the material presented and the references: 
Due to the flexibility of the NRG, the review covers 
a broad range of physical phenomena, in particular in 
Sees. QV] and |V] We shall, however, only give very brief 
introductions to these phenomena and refer the reader 
to the references given in the individual subsections, in 
particular reviews or seminal books. Furthermore, due to 
lack of space, we shall mostly not review the results from 
other theoretical approaches which have been applied to 
quantum impurity systems, such as Bethe ansatz, quan- 
tum Monte Carlo, resolvent perturbation theory, local- 
moment approach, etc, unless these appear crucial for an 
understanding of relevant NRG results. Comparisons be- 
tween the results from NRG and these approaches are, 
in most cases, included in the papers cited here (see also 
iHewsonl fl993> ). This means that we shall focus, almost 
completely, on references which use the NRG. 



II. INTRODUCTION TO THE NUMERICAL 
RENORMALIZATION GROUP METHOD 

The NRG method can be applied to systems of the 
following form: a quantum mechanical impurity with 
a small number of degrees of freedom (so that it can 
be diagonalized exactly) coupling to a bath of fermions 
or bosons, usually with continuous excitation spectrum. 
There is no restriction on the structure of the impurity 
part of the Hamiltonian; it might contain, for example, 
a Coulomb repulsion of arbitrarily large strength. The 
bath, however, is required to consist of non-interacting 
fermions or bosons, otherwise the various mappings de- 
scribed below cannot be performed. 

Whenever we discuss, in this review, models of a differ- 
ent kind, such as the Hubbard model, these models will 
be mapped onto impurity models of the above type. For 



the Hubbard model and other lattice models of correlated 
electrons this is achieved via the dynamical mean-field 
theory, see Sec. LY1 

Before we start with the technical details of the NRG 
approach, let us give a brief overview of the general strat- 
egy. For basically all NRG applications, one proceeds as 
follows: 

a) Division of the energy support of the bath spectral 
function into a set of logarithmic intervals. 

b) Reduction of the continuous spectrum to a discrete 
set of states (logarithmic discretization). 

c) Mapping of the discretized model onto a semi- 
infinite chain. 

d) Rerative diagonalization of this chain. 

e) Further analysis of the many-particle energies, ma- 
trix elements, etc., calculated during the iterative 
diagonalization. This yields information on fixed 
points, static and dynamic properties of the quan- 
tum impurity model. 

Parts a),b) and c) of this strategy are sketched in Fig. [TJ 
where we consider a constant bath spectral function 
within the interval [—1,1]. The NRG discretization pa- 
rameter A defines a set of discretization points, ±A _n , 
71 = 0,1, 2,..., and a corresponding set of intervals. The 
continuous spectrum in each of these intervals (Fig. [T^) 
is approximated by a single state (Fig. [TJj). The result- 
ing discretized model is mapped onto a semi-infinite chain 
with the impurity (filled circle) corresponding to the first 
site of this chain. Due to the logarithmic discretization, 
the hopping matrix elements decrease exponentially with 
distance from the impurity, t n oc A~ n / 2 . 

While the various steps leading to the semi-infinite 
chain are fairly straightforward from a mathematical 
point of view, the philosophy behind this strategy is prob- 
ably not so obvious. 

Quite generally, a numerical diagonalization of Hamil- 
tonian matrices allows to take into account the various 
impurity-related terms in the Hamiltonian, such as a lo- 
cal Coulomb repulsion, non-perturbatively. Apparently, 
the actual implementation of such a numerical diagonal- 
ization scheme requires some sort of discretization of the 
original model, which has a continuum of bath states. 
There are, however, many ways to discretize such a sys- 
tem, so let us try to explain why the logarithmic dis- 
cretization is the most suitable one here. As it turns out, 
quantum impurity models are very often characterized 
by energy scales orders of magnitudes smaller than the 
bare energy scales of the model Hamiltonian. If the ratio 
of these energy scales is, for example, of the order of 10 5 , 
a linear discretization would require energy intervals of 
size at most 10 -6 to properly resolve the lowest scale in 
the system. Since for a finite system the splitting of ener- 
gies is roughly inversely proportional to the system size, 
one would need of the order of 10 6 sites, which renders 
an exact diagonalization impossible. 
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a) 



A(w) 



1-+ 




A(w) 




•— o— o— o— o- 

FIG. 1 Initial steps of the NRG illustrated for the single- 
impurity Anderson model in which an impurity (filled circle) 
couples to a continuous conduction band via the hybridization 
function A(w); a) a logarithmic set of intervals is introduced 
through the NRG discretization parameter A; b) the contin- 
uous spectrum within each of these intervals is approximated 
by a single state; c) the resulting discretized model is mapped 
onto a semi-infinite chain where the impurity couples to the 
first conduction electron site via the hybridization V; the pa- 
rameters of the tight-binding model (see Eq. (J26J ) are e„ and 
t n . 



Apparently, the logarithmic discretization reduces this 
problem in that the low-energy resolution now depends 
exponentially on the number of sites in the discretized 
model. Of course, the accuracy of such an approach has 
to be checked by suitable extrapolations of the discretiza- 
tion parameters, in particular a A — > 1 extrapolation, 
which recovers the original continuum model. Very often 
it turns out that already for A of the order of 2 the re- 
sults are accurate to within a few percent and a A — * 1 
extrapolation indeed reproduces exact results, if these are 
available. 

However, this argument in favor of the logarithmic dis- 
cretization does neither explain the need for a mapping 
to a chain Hamiltonian as in Fig. [Ik, nor how the problem 
of an exponentially growing Hilbert space with increas- 
ing chain length is resolved. As far as the first point is 
concerned, an iterative diagonalization for the discretized 



model with a star geometry as in Fi g. Hb has b e en im - 
plemented for the spin-boson model ([Bulla et ali . l2005h . 
For reasons which are not yet completely clear, such a 
'star-NRG' is only partly successful. Let us just men- 
tion here that for a fermionic model such as the single- 
impurity Anderson model, the iterative diagonalization 
of the model in the semi-infinite chain form is much more 
convenient since one site of the chain can be added in 
each step without violating particle -hole symmetry (f or 
a detailed discussion of this point see. Bu lla et ali ( 2005f )). 

The quantum impurity model in the semi-infinite chain 
form is solved by iterative diagonalization, which means 
that in each step of the iterative scheme one site of the 
chain is added to the system and the Hamiltonian matri- 
ces of the enlarged cluster are diagonalized numerically. 
As already pointed out, without taking further steps to 
reduce the size of the Hilbert space this procedure would 
have to end for chain sizes of « 10. Here the renor- 
malization group concept enters the procedure through 
the dependence of the hopping matrix elements on the 
chain length, t n oc A - "/ 2 . Adding one site to the chain 
corresponds to decreasing the relevant energy scale by 
a factor ^fk. Furthermore, because the coupling t n to 
the newly added site falls off exponentially, only states 
of the shorter chain within a comparatively small energy 
window will actually contribute to the states of the chain 
with the additional site. This observation allows to in- 
troduce a very simple truncation scheme: after each step 
only the lowest lying N s many-particle states are retained 
and used to build up the Hamiltonian matrices of the next 
iteration step, thus keeping the size of the Hilbert space 
fixed as one walks along the chain. 

All these technical steps will be discussed in detail in 
the following. Let us briefly remark on the general set- 
up of this section. We would of course like to keep this 
section as general as possible, because it should serve 
as an introduction to the NRG technique, whose appli- 
cation to a variety of problems is then the subject of 
the remainder of this review. This quest for general- 
ity is, however, contrasted by the large variety of pos- 
sible impurity-bath interactions. Instead of presenting 
explicit formulae for all possible quantum impurity mod- 
els, we restrict ourselves to the single-impurity Ander- 
son model as a specific - and most important - example 
here. The origi nal introduction s to the technique for the 
Kondo model (IWilsonl. Il975al ) and the single- i mpurity 
Anderson model ( Krishna- murthv et ali , 1980a. b) were 
restricted to a constant bath density of states (or bet- 
ter, a constant hybridization function A(w) as defined 
below). Here, we consider a general frequency dependent 
A(a>) from the outset. This generalization is essential 
for various applications of the NRG (the soft-gap mod- 
els, see Sec. HV.C.21 and lattice models within DMFT, 
see Sec. [V} where the physics is largely determined by 
the frequency dependence of A (a;). If the hybridization 
function is non-zero for positive frequencies only, the ma- 
nipulations of the bat h degrees of freedo m equally hold 
for a bosonic bath, see Bulla et al. I (|2005l) . 
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In this section we cover the steps a), b), c) and d) 
of the list given above. Concerning the analysis of the 
data [step e) in the list] , we discuss the flow of the many- 
particle spectra and all related issues here. The calcula- 
tion of static and dynamic quantities will be described in 

secnrn 



A. Structure of the Hamiltonian 

The Hamiltonian of a general quantum impurity model 
consists of three parts, the impurity i?i mp , the bath 
-Hbath, and the impurity-bath interaction H i mp _bath : 



H — -ffimp + #bath + -ffirnp — bath 



(1) 



For the single- impurity Anderson model (SIAM) 
(Anderson), Il96lh with the Hamiltonian H = -ffsiAM, 
these three terms are given by: 

(J 

■Hbath = y ^^ £ kc\ (T Ck ( T , 
ka 

^imp-bath = Vfc (^fa c kcr + c \ a f° 



(2) 



In this Hamiltonian, the fermionic operators cj}' corre- 
spond to band states with spin a and energy Ek , and the 
fa^ to the impurity states with energy e{. The Coulomb 
interaction between two electrons occupying the impu- 
rity site is parametrized by U and both subsystems are 
coupled via a /c-dependent hybridization 

The influence of the bath on the impurity is completely 
determined by the so-called hybridization function A(u>): 



AH 



Efe) 



(3) 



Thus, if we are only interested in the impurity contri- 
butions to the physics of the SIAM, we can rewrite the 
Hamiltonian in a variety of ways, provided the manipula- 
tions involved do not change the form of A(u>). Without 
loss of generality, we assume that the support of A(w) 
completely lies within the interval [—D,D], with D > 
chosen suitably. Henceforth, we will use D = 1 as energy 
unit. 

One such possible reformulation is given by the follow- 
ing Hamiltonian: 



H = H imp + J2 [ deg(e) 

r, •> — 1 



W (4) 



Here we introduced a one-dimensional energy representa- 
tion for the conduction band with band cut-offs at ener- 
gies ±1, a dispersion g(s) and a hybridization h(e). The 



band operators fulfill the standard fermionic commuta- 
tion rela tions: [q^, a e > a >\ = 8 (e — e')S aa i . It has been 
shown in lBulla et all ( 1997bl ) that the two functions g(e) 
and h(e) are related to the hybridization function A(w) 
via 



AM 



de(w) 2 
dui 



(5) 



where s(lu) is the inverse function to g(e), g[e(u>)] = lu. 
For a given A(w) there are obviously many possibilities 
to divide the w-dependence between e(u) and h(e(u>)). 
This feature will turn out to be useful later. 

For a constant A(w) = Ao within the interval [—1, 1], 
Eq. ([5]) can be fulfilled by choosing s(u>) — u> (this corre- 
sponds to g(e) = e) and h 2 (e) = A /ir. 



B. Logarithmic discretization 

The Hamiltonian in the integral representation Eq. ^ 
is a convenient starting point for the logarithmic dis- 
cretization of the conduction band. As shown in Fig. [lji, 
the parameter A > I defines a set of intervals with dis- 
cretization points 



x n = ±A~ n , n = 0,1,2,. 
The width of the intervals is given by 
d n = A-™(1 - A^ 1 ) . 



(6) 



(7) 



Within each interval we now introduce a complete set of 
orthonormal functions 



^_ e ±iu n pe for Xn+i <±£< Xri 



^ np ( ^ 1 (; outside this interval . 



(8) 



The index p takes all integer values between — oo and 
+oo, and the fundamental frequencies for each interval 
are given by uj n — 2n/d n . The next step is to expand the 
conduction electron operators a ea in this basis, i.e. 



(9) 



which corresponds to a Fourier expansion in each of the 
intervals. The inverse transformation reads 



(10) 



The operators (Xnpa and bnpa defined in this way fulfill 
the usual fermionic commutation relations. The Hamil- 
tonian Eq. ((J) is now expressed in terms of these discrete 
operators. 
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In particular, the transformed hybridization term (first 
part only) is 



de h{e)fla e 



ft E 

up 
+ b 



^npa 



+,n 



d£h(e)ip+ _(e) 



npcr 



deh(e)ip (e) 



,(H) 



where we have defined 

c+,n 



de 



d£ 



de 



— X„ + l 



de . (12) 



For a constant h(e) = h, the integrals in Eq. (fTTj) filter 
out the p = component only 



±,n 



(13) 



In other words, the impurity couples only to the p = 
components of the conduction band states. It will be- 
come clear soon, that this point was essential in Wilson's 
original line of arguments, so we would like to main- 
tain this feature (h(e) being constant in each interval 
of the logarithmic discretization) also for a general, non- 
constant A(w). Note that this restriction for the func- 
tion /i(e) does not lead to additional approximations for 
a non-constant A(w) as one can shift all the remaining 
e-dependence to the di spersion g(e), see Eq. dp . 

As discussed in IChen and Javaprakashl (jl995al ) 
in the context of th e soft -gap model (see also 
IChen and Javaprakashl ( 1995bl )). one can even set 
h(e) = h f or all e. Here we follow the proposal of 
iBulla et all (Il997bl) . that is we introduce a step function 
for h{e) 



He) 



L n+1 



< ±e < x n 



(14) 



with given by the average of the hybridization func- 
tion A(u>) within the respective intervals, 



= — 

n d 



±,n 



de - A(e) . 



(15) 



This leads to the following form of the hybridization term 
f 1 1 

/ de h(e)fla ea = -=fl V [7+ a n0rT + j~b n0rT ] , 

J-l ^ n 

(16) 

with 7 ± 2 = / ± '"deA(e). 

Next, we turn to the conduction electron term, which 
transforms into 

J ^ de g(e)a\ cr a e(X = ^ (C4 P <r a «P<x + C & n P <A; 



npa a rip'a 



n,p^p' 



a n(P'P') b lp<r b np>a) 



The first term on the right hand side of Eq. (fl~7f is diag- 
onal in the inde x p. The discrete se t of energies can 
be expressed as ( Bulla et~al\ , Il997bl ) 



f'" deA(e)e 
J ± '"deA(e) 



-A- n (l + A- x ) 



(18) 



where we added the result for a constant A(e) in brack- 
ets. The coupling of the conduction band states with 
different p,p' (the second term) recovers the continuum 
(no approximation has been made so far, Eq. I|17p is still 
exact). For the case of a linear dispersion, g(e) — e, the 
prefactors a~(p,p') are the same for both sides of the 
discretization and take the following form 



<xt(p>P') 



I A 



A" 



2ni 



P -P 



■ exp 



2iri (p' — p) 



1- A- 



(19) 



The actual discretization of the Hamiltonian is now 
achieved by dropping the terms with p ^ in the expres- 
sion for the conduction band Eq. (fl"7|) . This is, of course, 
an approximation, the quality of which is not clear from 
the outset. To motivate this step we can argue that (i) 
the p 7^ states couple only indirectly to the impurity 
(via their coupling to the p = states in Eq. (JTTJ)) and 
(ii) the coupling between the p = and p ^ states has a 
prefactor (1 — A -1 ) which vanishes in the limit A — ► 1. In 
this sense one can view the couplings to the states with 
p 7^ as small parameters and consider the restriction to 
p = as zeroth order step in a perturb ation expansio n 
with respect to the coefficients a^(p,p r ) i Wilsonl . 1975a). 
As it turns out, the accuracy of the results obtained from 
the p = states only is surprisingly good even for values 
of A as large as A = 2, so that in all NRG calculations 
the p states have never been considered so far. 

Finally, after dropping the p =/= terms and relabel- 
ing the operators a n o a = aw, etc., we arrive at the dis- 
cretized Hamiltonian as depicted by Fig. [T]b 

H = H imp + [C a LoW + €n b tur b no] 
1 la 

+ 7y^ E ft E ( 7 « Qna + 7 « bna ) 



—y 



(20) 



Before we continue with the mapping of the Hamiltonian 
Eq. (|20[) onto a semi-infinite chain, Fig. [It, let us make a 
few remarks on alternative discretizations of the contin- 
uous bath spectral function. 

The above procedure obviously applies for general 
asymmetric A(w), also for different upper and lower cut- 
offs, D u and D\. A spe cial case is D\ = , which occurs 
for a bosonic bath, see Bulla et al. (2005); here the log- 
arithmic discretization is performed for positive frequen- 
cies only, and the operators bnl in Eq. (|2"D|) are no longer 



(17) present. 
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In Sec. Mil we shall discuss that the discreteness of 
the model Eq. (|2"Uj) can be (in some cases) problematic 
for the calculation of physical quantities. As it is not 
possible in the actual calculations to recover the contin- 
uum by taking the limit A — -> 1 (or by including the 
p =/= terms), it has been sugg ested to average over vari- 
ous discretizations for fixed A dFrota and Oliveiral [19861: 
lOliveira and Oliveiral 1 1994 lYoshida et all Il99dh . The 
discretization points are then modified as 



1 : n = 

A -(n+Z) . n>1 



(21) 



where Z covers the interval [0,1). This 'Z- trick' is, in- 
deed, successful as it removes certain artificial oscillations 
(see Sec. IIII.A.4[> . but it should be stressed here that the 
continuum limit introduced by integrating over Z is not 
the same as the true continuum limit A — * 1. 

Another shortcoming of the discretized model is that 
the hybridization function A(w) is systematically under- 
estimated. It is therefore convenient to multiply A (a;) 
with the correction factor 



A\ = - In A 



A + l 
l A^T 



(22) 



which accelerates the convergence to the continuum limit. 
For a recent derivation of th is correction factor see 
ICampo. Jr. and Oliveiral (|2005l ) , where it was also shown 
that by a suitable modification of the discretization pro- 
cedure, the factor can be taken into account from the 
outset. 



C. Mapping on a semi-infinite chain 

According to Fig. [TJd and c, the next step is to trans- 
form the discretized Hamiltonian Eq. (|2"D|) into a semi- 
infinite chain form with the first site of the chain (filled 
circle in Fig. [IJ;) representing the impurity degrees of 
freedom. In the chain Hamiltonian, the impurity directly 
couples only to one conduction electron degree of freedom 
with operators CqJ , the form of which can be directly read 
off from the second and third line in Eq. (j20|) . With the 
definition 



1 



COc 



v $o 



+ 7, 



n brio] 



(23) 



in which the normalization constant is given by 



(similarly for the hermitian conjugate term). Note that 
for a fc-independent hybridization, Vk = V in Eq. (J2J), the 
coupling in Eq. (j2"5| reduces to v^o/tt = V. 

The operators CqJ represent the first site of the con- 
duction electron part of the semi-infinite chain. These 
operators are of course not orthogonal to the operators 
flri2, bna- Constructing a new set of mutually orthog- 
onal operators Cni from CqJ and ana-, bnl by a stan- 
dard Gram-Schmidt procedure leads to the desired chain 
Hamiltonian, which takes the form 



H = H\ 



a 

oo 



an— 



(26) 



with the operators Cnd corresponding to the nth site of 
the conduction electron part of the chain. The parame- 
ters of the chain are the on-site energies s n and the hop- 
ping matrix elements t n . The operators cltr in Eq. (|26p 
and the operators {a^tr, &$} in Eq. ® are related via 
an orthogonal transformation 



m— Tci— 

oo 

Cna" — ^ [^nra^ma ~t~ ^nm^mfj] • (^^) 
tci— 

From the definition of coo- in Eq. (|2"3")) we can read off the 
coefficients uo m and UQ m 



7ro 7m IOQ\ 

~HF= , "Om — —j= ■ (28J 

For the remaining coefficients u nm , v nm , as well as for 
the parameters e n , t n , one can derive recursion relations 
following the sc heme described in detail in, for example, 
Appendix A of iBulla et al\ ( 2005f) . The starting point 
here is the equivalence of the free conduction electron 
parts 



[£n a na a n(T + £, n Hitybna] — 

oo 



(771 — 



The recursion relations are initialized by the equations 



1 f 1 

£0 = 7-/ deA(e)e , 
50 J-l 



£0 = J2 f + (7,7) 2 ) = I* deA(e) , (24) ° 1 



the hybridization term can be written as 



( 7n a na- 



- 7., b 



<n u tuj , 



flcoa 



(25) 



Ulr, 
Vlr, 



T(t.m — e o) U Om 7 

to 



to 



(£m ~ £o)vOr, 



(30) 
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For n > 1, the recursion relations read 



r = 



(£+u 2 +£~v 2 ) 

\Srn nrn 1 Sm ^nm ) ' 

\(£ + \ 2 v 2 + f£"W I-/ 2 -f 2 

LVSm/ ^nm ~ \Sm/ u nmJ °n— 1 °n ' 



"-n+l,m 



r [(^-£n)t 



Cn— l^n— l,r, 



, [(^m ^ n') / V nrn tfi— 1 ^n— 1 , r 
Cri. 



(31) 



Note that for a particle-hole symmetric hybridization 
function, A(u>) = A(—u>), the on-site energies e n are zero 
for all n. 

For a general hybridization function, the recursion re- 
lations have to be solved numerically. Although these 
relations are fairly easy to implement, it turns out that 
the iterative solution breaks down typically after about 
20-30 steps. The source of this instability is the wide 
range of values for the parameters entering the recursion 
relations (for instance for the discretized energies 
In most cases this problem can be overcome by using ar- 
bitrary precision routines for the numerical calculations. 
Furthermore, it is helpful to enforce the normalization of 
the vectors u nm and v nm after each step. 

Analytical solutions for the recursion relations have so 
far been given only for few special cases. Wilson derived 
a formula for the t n for a constant density of states of the 
conduc tion elec t rons in the Kondo version of the impurity 
model (|Wilsonl . Il975al) : this corresponds to a constant 
hybridization function A(u) in the interval [—1, 1]. Here 
we have e„ = for all n and the expression for the t n 
reads 



tn 



1 + A- 



1-A 



-71-11 



2 VI - A- 2n - 1 Vl 3 A = 



(32) 



(Similar ex pressions h a ve been given for the soft-gap 
model, see iBulla et Ml (ll997blU In the limit of large 
n this reduces to 



tn 



^(l + A- 1 ) A-™/ 2 . 



(33) 



The fact that the t n fall off exponentially with the dis- 
tance from the impurity is essential for the following dis- 
cussion, so let us briefly explain where this n-dependence 
comes from. Consider the discretized model Eq. ([20]) with 
a finite number 1 + Mj 2 (M even) of conduction electron 
states for both positive and negative energies (the sum 
over n then goes from to M/2). This corresponds to 
2 + M degrees of freedom which result in 2 + M sites 
of the conduction electron part of the chain after the 
mapping to the chain Hamiltonian. The lowest energies 
in the discretized model Eq. (|2U)l are the energies £ M / 2 
which, for a constant hybridization function, are given by 



± l A -M/2 ( - 1 + A -l^ gee Eq (pj) Thig energy 



M/2 — ""2" 

shows up in the chain Hamiltonian as the last hopping 
matrix element t^j, so we have tyi ~ Cm/2 equivalent to 
Eq. 



Equation (|26[) is a specific one-dimensional representa- 
tion of the single-impurity Anderson model Eq. ^ with 
the special feature that the hopping matrix elements t n 
fall off exponentially As mentioned above, this represen- 
tation is not exact since in the course of its derivation, 
the p terms have been dropped. We should stress 
here that the dimensionality of the chain Hamiltonian is 
not related to that of the original model which describes, 
for example, an impurity in a three-dimensional host (ap- 
parently, this only holds for a non-interacting conduction 
band). Nevertheless, the conduction electron sites of the 
chain do have a physical meaning in the original model 
as they can be viewed as a sequence of shells centered 
around the impurity. The first site of the conduction elec- 
tron chain corresponds to the shell with t he maximum of 
its wav e functio n closest to the impurity (jHewsonl . Il993t 
IWilsonl . Il975al) : this shell is coupled to a shell further 
away from the impurity and so on. 



D. Iterative diagonalization 

The transformations described so far are necessary to 
map the problem onto a form (the semi-infinite chain, 
Eq. (|26p) for which an iterative renormalization group 
(RG) procedure can be defined. This is the point at 
which, finally, the RG character of the approach enters. 

The chain Hamiltonian Eq. (j2l)|) can be viewed as a 
series of Hamiltonians Hm (N = 0,1,2,...) which ap- 
proaches H in the limit N — > oo. 



H= lim A-^-^/'H 



N- 



N 



(34) 



with 



N 



N-X 



(771 — 



an—O 



(35) 



The factor A^ -1 )/ 2 in Eq. (I35[) (and, consequently, the 
factor A - ^ -1 )/ 2 in Eq. (f3"4")) ) has been chosen to can- 
cel the iV-dependence of ijv-i, the hopping matrix ele- 
ment between the last two sites of Hpj. Such a scaling 
is useful for the discussion of fixed points, as described 
below. For a diff erent n-dependenc e of t n , as for the 
spin-boson model (jBulla et all 12005). the scaling factor 
has to be changed accordingly. (The n-dependence of e n 
is, in most cases, irrelevant for the overall scaling of the 
many-particle spectra.) 

Two successive Hamiltonians are related by 

H N+ i = VAH N + A N/2 ^ ew+icjv+icrCw+icr 

c jV<x C JV+l<T + c ^N+la C N<j] i (36) 
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and the starting point of the sequence of Hamiltonians is 
given by 



H = A" 1 / 2 



Himp + £o c o<t C 0<t 



(37) 



This Hamiltonian corresponds to a two-site cluster 
formed by the impurity and the first conduction electron 
site. Note that in the special case of the single-impurity 
Anderson model, one can also choose H-i = K^ 1 H lrav 
as the starting point (with a proper renaming of param- 
eters and operators) since the hybridization term has the 
same structure as the hopping term between the conduc- 
tion electron sites. 

The recursion relation Eq. (|36|) can now be understood 
in terms of a renormalization group transformation R: 



H 



N+l 



R(H 



N) 



(38) 



In a standard RG transformation, the Hamiltonians are 
specified by a set of parameters K and the mapping R 
transforms the Hamiltonian H{K) into another Hamilto- 
nian of the same form, H(K'), with a new set of parame- 
ters K' . Such a representation does not exist, in general, 
for the Hn which are obtained in the course of the iter- 
ative diagonalization to be described below. Instead, we 
characterize Hn, and thereby also the RG flow, directly 
by the many-particle energies En{t) 



H n \t)n = E N {r)\r) N 



1,. 



(39) 



with the eigenstates |r)jy and N s the dimension of Hn- 
This is particularly useful in the crossover regime be- 
tween different fixed points, where a description in terms 
of an effective Hamiltonian with certain renormalized pa- 
rameters is not possible. Only in the vicinity of the fixed 
points (except for certain quantum critical points) one 
can go back to an effective Hamiltonian description, as 
described below. 

Our primary aim now is to set up an iterative scheme 
for the diagonalization of Hn, in order to discuss the 
flow of the many-particle energies En{t). Let us assume 
that, for a given N, the Hamiltonian Hn has already 
been diagonalized, as in Eq. (f3"9"| . We now construct a 
basis for Hn+i, as sketched in Fig.O 



\r;s) N +i = |r)jv ® \s(N + 1)) . 



(40) 



The states |r; s)n+i are product states consisting of the 
eigenbasis of Hn and a suitable basis \s(N + 1)) for the 
added site (the new degree of freedom). From the basis 
Eq. (|40[) we construct the Hamiltonian matrix for Hn+i- 



H N+ i(rs,r's') = N+i(r; s\H N+1 \r'; s') N+1 



(41) 



For the calculation of these matrix elements it is useful 
to decompose Hn+i into three parts 



H 



N+l 



AHjy + Xn,N+1 + Y/v+1 



(42) 



|r,s) 




® |S(N+1)> 



H - : O— O- 



-o-o 



FIG. 2 In each step of the iterative diagonalization scheme 
one site of the chain (with operators c N \ j and on-site energy 
EN+i) is added to the Hamiltonian Hn- A basis \r; s)n+i for 
the resulting Hamiltonian, Hn+i, is formed by the eigenstates 
of Hn, |r)jv, and a basis of the added site, \s(N + 1)). 



MO 



b) 



a E„(r) 



d) 



after truncation 



FIG. 3 (a): Many-particle spectrum J5jv(r) of the Hamilto- 
nian Hn with the ground-state energy set to zero, (b): The 
relation between successive Hamiltonians, Eq. (|36p . includes 
a scaling factor y/X. (c) Many-particle spectrum Sjv+i(r) 
of Hn+i, calculated by diagonalizing the Hamiltonian ma- 
trix Eq. (|41|l . (d) The same spectrum after truncation where 
only the N s lowest-lying states are retained; the ground-state 
energy has again been set to zero. 



(see, for example, Eq. (|3"fj)0 where the operator Yn+i 
only contains the degrees of freedom of the added site, 
while Xjv,iv+i mixes these with the ones contained in 
Hn- Apparently, the structure of the operators X and Y, 
as well as the equations for the calculation of their matrix 
elements, depend on the model under consideration. 

The following steps are illustrated in Fig. [3J In Fig. [3^ 
we show the many-particle spectrum of Hn, that is the 
sequence of many-particle energies En(t). Note that, 
for convenience, the ground-state energy has been set to 
zero. Figure [SJd shows the overall scaling of the energies 
by the factor vA, see the first term in Eq. (|36|) . 

Diagonalization of the matrix Eq. (|4"Tj) gives the new 
eigenenergies En+i(w) and eigenstates \w)n+i which are 
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related to the basis \r; s) jv+i via the unitary matrix U : 



\w)n+i =^2u(w,rs)\r;s) 



N+l 



(43) 



The set of eigenenergies Ejy + i(w) of H^+i is displayed 
in Fig. [3fc (the label w can now be replaced by r). Ap- 
parently, the number of states increases by adding the 
new degree of freedom (when no symmetries are taken 
into account, the factor is just the dimension of the basis 
|s(iV + l))). The ground-state energy is negative, but will 
be set to zero in the following step. 

The increasing number of states is, of course, a prob- 
lem for the numerical diagonalization; the dimension of 
-ff/v+i grows exponentially with N, even when we con- 
sider symmetries of the model so that the full matrix 
takes a block-diagonal form with smaller submatrices. 
This problem can be solved by a very simple truncation 
scheme: after diagonalization of the various submatrices 
of Hn+i one only keeps the N s eigenstates with the low- 
est many-particle energies. In this way, the dimension 
of the Hilbert space is fixed to N s and the computation 
time increases linearly with the length of the chain. Suit- 
able values for the parameter N s depend on the model; 
for the single-impurity Anderson model, N s of the order 
of a few hundred is sufficient to get converged results for 
the many-particle spectra, but the accurate calculation 
of static and dynamic quantities usually requires larger 
values of N s . The truncation of the high energy states is 
illustrated in Fig. [3U. 

Such an ad-hoc truncation scheme needs further ex- 
planations. First of all, there is no guarantee that this 
scheme will work in practical applications and its quality 
should be checked for each individual application. Impor- 
tant here is the observation that the neglect of the high- 
energy states does not spoil the low-energy spectrum in 
subsequent iterations - this can be easily seen numeri- 
cally by varying N s . The influence of the high-energy on 
the low-energy states is small since the addition of a new 
site to the chain can be viewed as a perturbation of rel- 
ative strength A -1 / 2 < 1. This perturbation is small for 
large values of A but for A — > 1 it is obvious that one has 
to keep more and more states to get reliable results. This 
also means that the accuracy of the NRG results is get- 
ting worse when N s is kept fixed and A is reduced (vice 
versa, it is sometimes possible to improve the accuracy 
by increasing A for fixed N s ). 

From this discussion we see that the success of the 
truncation scheme is intimately connected to the special 
structure of the chain Hamiltonian (that is t n oc A~™/ 2 ) 
which in turn is due to the logarithmic discretization of 
the original model. A direct transformation of the single- 
impurity Anderson model to a one -dimensional chain re- 
sults in t n — > const (|Hewsonl . Il993f ). and the above trun- 
cation scheme fails. A similar observation is made when 
such a truncation is applied to the one-dimensional Hub- 
bard model, see the brief discussion in Sec. W\ 

Let us now be a bit more specific on how to construct 
the basis \t;s)n+i- For this we have to decide, first of 



all, which of the symmetries of the Hamiltonian should 
be used in th e iterative diago nali zation. In the original 
calculations o f Wilson (|l975al ) and lKrishna-murthy et all 
(1980a, b) the following quantum numbers were used: to- 
tal charge Q (particle number with respect to half- filling) , 
total spin S and z-component of the total spin S z . It has 
certainly been essential in the 1970's to reduce the size 
of the matrices and hence the computation time as much 
as possible by invoking as many symmetries as possi- 
ble. This is no longer necessary to such an extent on 
the modern computer systems, i.e. one can, at least for 
single-band models, drop the total spin S and classify the 
subspaccs with the quantum numbers (Q, S z ) only. This 
simplifies the program considerably as one no longer has 
to worry about reduced matrix elements and the corre- 
sponding Clebsch-Gordan coe fficients, see, for example 
iKrishna-murthv et all ( 1980al) . As we use this represen- 
tation in Sec. IIII.A| let us here explicitly state the form 
of \r; s) N+1 : 



\Qi S z , f; 1) jv+i 
\Qi S z , t; 2) N+1 

\Q: S z , T ', 3) jv+l 

\Q,Sz,r;4) N+1 



\Q + l,S z 

c jv+it 



'A' 



J 

J o 

-JV+lT^jV+U 



\Q,S Z 2 ' r ) 

\Q,S z + \,r) 

1 u \Q-i,S z 



N 



N 



(44) 



' N 



Note that the quantum numbers (Q,S Z ) on both sides 
of these equations refer to different systems, on the left- 
hand side they are for the system including the added 
site, and on the right-hand side without the added site. 
We do not go into the details of how to set up the Hamil- 
tonian matrices Eq. (|4"Tj) , as this procedure is described 
in grea t detail in Appendix B in IKrishna-murthv et all 
(|l980af )). 

For fermionic baths, the discretization parameter A 
and the number of states -/V s kept in each iteration are 
the only parameters which govern the quality of the re- 
sults of the NRG procedure. As discussed in more detail 
in Sec. IIV.E1 for the case of a bosonic bath the infinite 
dimensional basis \s(N + 1)) for the added bosonic site 
requires an additional parameter At,, which determines 
the dimension of \s(N + 1)). 



E. Renormalization group flow 

The result of the iterative diagonalization scheme are 
the many-particle energies E^{r) with r — 1, . . . , N s (ap- 
parently, the number of states is less than N s for the very 
first steps before the truncation sets in). The index N 
goes from to a maximum number of iterations, A max , 
which usually has to be chosen such that the system has 
approached its low-temperature fixed point. 

As illustrated in Fig. [31 the set of many-particle ener- 
gies cover roughly the same energy range independent of 
N, due to the scaling factor A^' 1 ^ 2 in Eq. ([35]). The 
energy of the first excited state of Hm is of the order 
of A^-V/Hff-i, a constant according to Eq. ([55]). The 
energy of the highest excited state kept after truncation 
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FIG. 4 Flow of the lowest-lying many-particle energies of the 
single-impurity Anderson model for parameters e/ = —0.5 ■ 
10" 3 , U = 10~ 3 , V = 0.004, and A = 2.5. The states are 
labeled by the quantum numbers total charge Q and total 
spin 5*. See the text for a discussion of the fixed points visible 
in this plot. 



depends on N s - for typical parameters this energy is 
approximately 5-10 times larger as the lowest energy. 

Multiplied with the scaling factor A - ^ -1 )/ 2 , see 
Eq. J331), the energies En(t) are an approximation to 
the many-particle spectrum of the chain Hamiltonian 
Eq. (|26[) within an energy window decreasing exponen- 
tially with increasing N. Note, that the energies for 
higher lying excitations obtained for early iterations are 
not altered in later iteration steps due to the truncation 
procedure. Nevertheless one can view the resulting set 
of many-particle energies and states from all NRG iter- 
ations N as approximation to the spectrum of the full 
Hamiltonian and use them to calculate physical proper- 
ties in the whole energy range, see Sec. IIII1 

Here we want to focus directly on the many-particle 
energies En(t) and show how one can extract informa- 
tion about the physics of a given model by analyzing their 
flow, that is the dependence of Ejsr(r) on N. 

As a typical example for such an analysis, we show 
in Fig. @] the flow of the many-particle energies for the 
symmetric single-impurity Anderson model, with param- 
eters e f = -0.5 • 10- 3 ,~ U = 10~ 3 , V = 0.004, and 
A = 2.5 (the same p arame ters as used in Fig. 5 in 
iKrishna-murthv et all ( 1980al ); note that we show here 
a slightly different selection of the lowest-lying states). 
The energies are plotted for odd N only, that is an odd 
total number of sites (which is N + 2). This is neces- 
sary, because the many-particle spectra show the usual 
even-odd oscillations of a fermionic finite-size system (the 
patterns for even TV look different but contain, of course, 
the same physics). The data po ints are connected by 
lines to visualize the flow. As in IKrishna-murthv et all 
(1980a), the many-particle energies are labeled by total 
charge Q and total spin S. 

What is the information one can extract from such a 
flow diagram? First of all we note the appearance of 



three different fixed points of the RG transformation for 
early iteration numbers N < 10, for intermediate values 
of N and for N > 60 (strictly speaking, because we look 
at N odd only, these are fixed points of R 2 , not of R). 
The physics of these fixed points cannot be extracted by 
just looking at the pattern of the many-particle energies. 
This needs some further analysis, in particular the di- 
rect diagonalization of fixed point Hamiltonians (which 
usually have a simple structure) and the comparison of 
their spectrum with the numerical data. An excellent 
account of this procedure for the symmetric and asym- 
m etric single-impurity Anderson model has been given 
bv lKrishna-murthv et all (|l980al lbl). and there is no need 
to repeat this discussion here. The analysis shows that 
for N w 3 — 9, the system is very close to the free-orbital 
fixed point, with the fixed point Hamiltonian given by 
Eq. (IB) for U = and V = 0. This fixed point is unsta- 
ble and for N ss 11 — 17, we observe a rapid crossover to 
the local- moment fixed point. This fixed point is char- 
acterized by a free spin decoupled from the conduction 
band. The local-moment fixed point is unstable as well 
and after a characteristic crossover (see the discussion 
below) the system approaches the stable strong-coupling 
fixed point of a screened spin. Note that the terminology 
'strong-coupling' has been introduced originally because 
the fixed point Hamiltonian can be obtained from the 
limit V — > oo, so 'coupling' here refers to the hybridiza- 
tion, not the Coulomb parameter U. 

The NRG does not only allow to match the structure 
of the numerically calculated fixed points with those of 
certain fixed point Hamiltonians. One can in addition 
identify the deviations from the fixed points (and thereby 
part of the crossover) with appropriate perturbations of 
th e fixed point Hamilton i ans. A gain, we refer the reader 
to lKrishna-murthv et all ( 1980al rj) for a detailed descrip- 
tion of this analysis. The first step is to identify the 
leading perturbations around the fixed points. The lead- 
ing operators can be determined by expressing them in 
terms of the operators which diagonalize the fixed point 
Hamiltonian; this tells us directly how these operators 
transform under the RG mapping R 2 . One then proceeds 
with the usual classification into relevant, marginal, and 
irrelevant perturbations. The final results of this analysis 
perfectly agree with the flow diagram of Fig. 2) There is a 
relevant perturbation which drives the system away from 
the free-orbital fixed point, but for the local-moment 
fixed point there is only a marginally relevant perturba- 
tion, therefore the system only moves very slowly away 
from this fixed point. Note that this marginal perturba- 
tion - which is the exchange interaction between the local 
moment and the spin of the first conduction electron site 
- gives rise to the logarithms observed in various physi- 
cal quantities. Finally, there are only irrelevant operators 
which govern the flow to the strong-coupling fixed point. 
These are responsible f or the Fermi-li quid properties at 
very low temperatures ( Hewsonl . Il993|) . 

Having identified the leading operators for each fixed 
point, it is possible to calculate physical properties 
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close to the fixed points perturbatively. We do not 
want to go into the calculational details here, see 



Crishna-murthv et al\ (Il980at) and a lso Se c. 4 in 
1993[) . Recently. iHewson et all |2004) and 



Hewson 



Hewson 



(2005) developed an alternative approach based on the 
renormalized perturbation theory. This approach is 
much easier to implement, has been used to describe the 
physics close to the strong-coupling fixed point and is, 
in principle, app licable also on all en ergy scales and for 
non-equilibrium (jHewson et all , [20051 ) . 

Flow diagrams as in Fig. 3] also give information about 
the relevant energy scales for the crossover between the 
fixed points. For example, an estimate of the Kondo tem- 
perature Ik (the temperature scale which characterizes 
the flow to the strong-coupling fixed point) is given by 
Tk ~ uj c Ar N f 2 , with N « 55 for the parameters in Fig. [4] 

The discussion of flow-diagrams as in Fig. 5] concludes 
our introduction to the basics of the NRG approach. An 
important part is still missing, of course, that is the cal- 
culation of physical quantities from the flow of the many- 
particle energies (and from certain additional matrix el- 
ements). This is the topic of the following section. 

In Sec. IIVI we will come back to the discussion of flow 
diagrams and the structure of fixed points when study- 
ing various other quantum impurity systems, in particu- 
lar the two-channel Kondo model which displays a non- 
Fermi liquid fixed point, see Sec. IIV.B1 and the soft-gap 
Anderson model which has a quantum critical point sep- 
arating the strong-coupling and local-moment phase, see 
Sec. IIV.C.2I 



III. CALCULATION OF PHYSICAL PROPERTIES 

In the previous section [TT] we discussed the information 
that can be gained from the the low-lying energy levels 
during the RG flow. Apparently, a lot can already be 
learned on this level about the physical properties of the 
system. However, an obvious aim of any method is also 
to calculate thermodynamic quantities like specific heat, 
susceptibilities or even dynamical properties. 

Let us start by reminding the reader that the co- 
efficients t n appearing in the transformed Hamiltonian 
Eq. (|26[) decay like A - ™/ 2 for large n. This aspect can be 
used in the following way (iKrishna-murthy et aLl . ll~980at 
lOliveira and Ohveiral I1994T Wilson . 1975al) : Diagonaliz- 
ing the Hamiltonian (|3"5|) for a given chain length N yields 



a set of eigenvalues rj\ 



(JV) 



±A Z . Obviously, eigenval- 



ues A ( N 1 )/ 2 ^ JV - ) ^> fc B T, to will not contribute signif- 
icantly to the calculation of physical properties anyway. 
On the other hand, for those I where K~^ N ~ l ^^m N ' <C 

k^T, lj one can safely approximate 77. ~ 0, which 
means that for the calculation of impurity properties 
these contributions will drop out. With (3 = (/cbT) -1 , 



^A 



-(N-l)/2 _ 



=:[3 



(45) 



Hn instead of the full Hamiltonian ([261) . 

Provided we can keep enough states in the truncation 
scheme introduced in Sec.[TT]to ensure convergence of the 
partition function on the scale k^T, it is thus permissible 
to use the truncated Hamiltonian on the level N obtained 
from the iterative diagonalization to calculate physical 
properties for the impurity on the temperature or energy 
scale A-^^^/p. 



A. Thermodynamic and static properties 

1. Entropy, specific heat and susceptibility 

The simplest physical quantities related to the impu- 
rity degrees of freedom are the impurity contribution to 
the entropy, Simp, specific heat, Ci mp , and magnetic sus- 
ceptibility, Ximp- 

The entropy and specific heat are the first derivative 
of the free energy F = -feTlnZ and internal energy 
U = (H) with respect to temperature, i.e. 

Q __dF 

,b - w 



and 



C = 



dU 
df 



From a numerical point of view, performing differentia- 
tions is something to avoid if possible. For the numerical 
implementation of the NRG another complication arises. 
To avoid an exponential increase of energies, it is neces- 
sary to subtract the ground state energy at each NRG- 
level N, i.e. one would have to keep track of these sub- 
tractions. Apparently, a much more convenient approach 
is to evaluate the derivative analytically, yielding 

S/k B = (3(H) + In Z , 

for the entropy and 

C/k B = 1 [(H 2 ) - (Hf\ , 

for the specific heat. 

The prescription to calculate the impurity contribu- 
tion to the magnetic susceptibility requires some more 
thought. The standard definition for the magnetic sus- 
ceptibility is (we set g/iB = 1) 



X (T)= I {S z [r]S z )dT - (3{S z y 




(46) 



with 



and (3 chosen properly, it will thus be sufficient to use 



(S z [r]S z ) = -Tr [e-f >a e rH S,e- TH S,] . 

However, the evaluation of the latter expectation value 
is equivalent to the calculation of a dynamical correla- 
tion function. This is in general a much more complex 
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task and will be discussed in detail in the next section. 
Here, we employ a different approach, which in turn is 
also closer related to the experimental definition of this 
quantity. 

In general, experiments address the susceptibility of 
the whole system. Since the total spin commutes with 
the Hamiltonian, the expression (I46p simplifies to 



obtained from the bare conduction Hamiltonian 



Xtot(T) = (3 [<$ 



2 ) 

tot, zl 



(S t , 



in this case. From this, one subtracts the susceptibility 
of a reference syste m, i.e. without impurity, leading to 
Wilson's definition ( W ilsonl . Il975at) of the impurity con- 
tribution to the susceptibility 



Ximp(T) = Xtot(T) - xl°ol{T) 



(47) 



Since 5tot,z is a quantum number used to classify the 
states in the calculation, the expectation values in |47|) 
can be evaluated straightforwardly. 

Similarly, the impurity contributions to the entropy 
and specific heat can be calculated as 



and 



where s[„l (T) and (T) are again entropy and specific 
heat of a suitable reference system. 

Let us discuss the details of the actual calculation for 
the entropy as specific example. Following the introduc- 
tory remarks, we can - for a given temperature k B T - 
restrict the Hilbert space to the NRG iteration L fulfill- 
ing (|45|). If we denote the corresponding Hamiltonian by 
H^ N \ we can introduce the quantity 



Sim V {T) — S to t{T) - S t °l(T) , 



C imp (T) = C tot (T)-C^(T) 
,(o), 



(48) 



(49) 



SW/fe :=/3(ff (Ar) ) (i) +hii? W , 



(50) 



where, using the notation of Sec. |TT] (see, for example, 
Eq. p])), 



EE 



,-/3E L (Q,S z ,r). 



Q,S Z 



N(Q,S z ,r\ . . . \Q,S z ,r) N 



and 



Z {N) ._ ^ e -/3£«(QA,r) 

Q,S 2 r 



(51) 
(52) 



The impurity contribution to the entropy for a tempera- 
ture /cbT/v := A - ^ -1 )/ 2 / '(3 can then be obtained as 



SW{Jj»)/*BWS W /*B-"cb 

Here we introduced the "free entropy" 



(53) 



(54) 



H. 



(N) 



an—Q 



(55) 



Similarly, for k B T N = A"^- 1 )/ 2 /^ the specific heat 
and magnetic susceptibility are obtained as 



C imp (T N )/k B ^C^-C^ , 



and 



Ximp(Tjv) 



(N) 
Xtot 



(JV) 
Xcb 



(56) 



(57) 



Since the Hamiltonian (|55p is a non-interacting sys- 
tem, these quantities o^P etc. can be expressed via the 



eigenenergies r\i a of (|55[) in standard fashion. 

For T — > the behavior of S\ mp (T) and Xi mp (T) given 
by Eqs. (|56"|) and (|57p can be obtained analytically from 
the fixed point sp ectra. We refer t he reader interested in 
this derivation to IWilsonl (|!975ah and concentrate here 
on the actual numerical calculations. 

Another aspect is that the fixed points and the flow 
to them are different for N even and odd. This in turn 
means, that one in principle has to calculate thermody- 
namic properties either for N even or odd only and thus 
loose half of the temperature values. One can, however, 
use all information by properly averaging odd and even 
steps: 

• For a given N, calculate the quantities O^ -1 ), 
OWandO^ 1 ). 

• Approximate O(Tjv) as 



0(T N ) 



1 r 

2 



(N) +0 (N-1) + 

Tn+i — Tjf-i 



(TV — Tjv-i) 



The first term in the square bracket is the observ- 
able calulated at step N. The second and third 
term are a linear interpolation of the values at N— 1 
and N + 1 to iteration N. 



Conti 



with N+ 1. 



As a positive side effect, this averaging also improves the 
accuracy of the thermodynamic quantities calculated. 

At this point some remarks about potential numeri- 
cal problems one can encounter should be made. The 
arguments given in the introduction to this section rely 
on the assumption, that one can keep states with high 
enough energy to ensure (i) the accuracy of the states 
at medium and low energies and (ii) the convergence of 
the partition function and expectation values. Depend- 
ing on the actual quantity to be calculated, the latter 
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point can in principle lead to problems. As an exam- 
ple, consider (H) and (H 2 ). While for a given energy 
cut-off E cut the contribution (3E CVLt e~ l3Ecut to (H) can 
already be small enough to use the sum up to E cut as 
approximation to (H), this must not be necessarily true 
for (H 2 ). Thus, the resulting values for the specific heat, 
(3 [(H 2 ) — (-H") 2 ], can be rather poor although entropy 
and susceptibility come out much more accurate. 



2. Other local properties 

While entropy, specific heat and impurity suscepti- 
bility can be obtained directly from the spectra of the 
Hamiltonian, other local quantities require the calcula- 
tion of the corresponding local matrix elements. As an 
example, we want to discuss here the local occupancy 
n <y = {ft fa) and double occupancy D = (f\f^f\f{) for 
the single- impurity Anderson model Eq. I[2]). Both quan- 
tities are of interest in actual applications. Expectation 
values of other local operators can be calculated in a sim- 
ilar manner. 

As before, on a given temperature scale fceT/v = 
A - ^ -1 " 2 //3, we approximate the expectation values by 



n a {T L 



(N) EE 



zw 



,S Z r 



N 



-0E N (Q,S g ,r) 



(Q,S z ,r\ttf*\Q,S z ,r) 



(58) 



for the occupancy and a corresponding expression for the 
double occupancy. The matrix elements 

n„{Q,S z ,r,r';N) := N {Q,S z ,r\flf a \Q,S z ,r') N , (59) 

at a given step N can be calculated from those of 
the previous step N — 1 with the help of the basis 
transformation (|4"5|) for the step N. The same scheme 
works for the matrix elements of the double occupancy 
D(Q, SztWjW 1 ; — 1) and the matrix elements of general 
local operators - like /J needed in the calculation of the 
single-particle Green function (see Sec. IIII.Bj) . 

All we are left to specify are the initial values for 
n a {Q, S z ,w, w'] —1) and D(Q, S z ,w, w'; —1) on the level 
of the impurity. For the Anderson model Eq. @ they 
are explicitly given as 



n a (0, 0,0,0 


-1) 


= o , 


n CT (l, a, 0, 


-1) 


= 1 , 


n CT (2, 0,0,0 


-1) 


= 2 , 


£1(0,0,0,0 


-1) 


= o , 


L>(l,cr,0,0 


-1) 


= o , 


£(2,0,0,0 


-1) 


= 1 . 



(60) 



With these prerequisites we are now in the position to do 
actual calculations for the thermodynamic properties of 
quantum impurity models using NRG. 



3. Example: The Kondo model 

As example for the method let us present results for 
the Kondo model Eq. (|128[) . Depending on the number of 
bands coupling to the local spin, one observes a conven- 
tional Kondo effect with the formation of a local Fermi 
liquid or a non-Fermi liquid fixed point with anomalous 
temperature dependencies of specific heat and suscep- 
tibility as well as a residual entropy S(0) = ^ l n2 at 
T = (|Cragg et all Il980t iNozieres and Blandirj, Il980(l 
(this will be discussed in more detail in Sees. IIV.AI and 
HV.Bp . In Fig. [5]we show the entropy 5i mp (T), suscepti- 




T/T„ 



FIG. 5 Entropy Sim P (T), susceptibility Ximp(T), Sommerfeld 
coefficient 7i mp = Ci mp (T)/T, and Wilson ratio Rw for the 
single-channel Kondo model. The Kondo temperature is de- 
fined by the Wilson relation Ximp (0) = 4j- • 

bility Xi m p(T), Sommerfeld coefficient 7i mp = Ci mp (T)/T 
and Wilson ratio Rw := 47r 2 Xi mp (T)/(37i mp (£)) as func- 
tion of T/Tk for the single-channel Kondo model. As 
Kondo coupling we choose J = 0.05D, where D is the 
half-bandwidth of the conduction band, for which we as- 
sume a density of states p c b( e ) = NpQ(D — |e |). The 
value o f Xk is obtained from Wilson's definition ( Wilsonl . 
Il975al ) 4T K Xim P (0) = 0.413. The calculations are per- 
formed with a discretization parameter A = 4, keeping 
400 states at each NRG step. Although this value of A 
seems to be fairly large, experience tells that for static 
properties such large values of A are still permissible, con- 
siderably reducing the number of states one has to keep 
in the truncation procedure. 

One nicely sees in Fig. [5] the quenching of the local mo- 
ment by the Kondo effect for temperatures of the order 
of Tk- Also the high-temperature values for the entropy 
ffjmpfX — > oo ) = hi 2 and the Wilson ratio Rw = 2 
(|Wilsod . fl975al) are obtained with high precision. 

If one adds a second screening channel to the Kondo 
model, one arrives at the so-called two-channel Kondo 
model. The rather interesting physics of this model will 
be discussed in detail in Sec. IIV.BI Here we merely want 
to demonstrate that NRG calculations for this model are 
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FIG. 6 Entropy S l i mp (T), susceptibility Xim P {T), Sommerfeld 
coefficient 7imp = d mp (T)/T, and Wilson ratio 7?w for the 
two-channel Kondo model. The value for the Kondo temper- 
ature is the same as in Fig. [5] 



possible, too; however, the additional bath degrees of 
freedom, which lead to Hilbert spaces larger by a fac- 
tor of four, make calculations more cumbersome and for 
some quantities also less accurate. In Fig. [5] we show 
as before the impurity contributions to the entropy, sus- 
ceptibility and Sommerfeld coefficient as well as the Wil- 
son ratio as function of T/Xk- The impurity parameters 
are the same as in Fig. [5] for the NRG we again choose 
A = 4 but keep 8100 states per iteration. The value 
of 2~k is that of the corresponding single-channel model. 
As emphasized before, entropy and susceptibility come 
out quite accurately, in particular the residual entropy 
5(0) = i In 2 is obtained as well as the l ogarithmic in 



crease of Y inip (T) e x In (T/T K ) for T < Ty (lAffleck et all . 
Il992t ICragg et all \l98& iPang and CoIll99lh . The~ 



spe- 
cific heat, however, is less accurate, but also shows the 
logarithmic increase as expected, although with strong 
oscillations superimposed. Fitting both quantities with 
a logarithmic form, one can recover the correct Wilson 
ration Rw = 8/3 for T — > 0. Note, however, that the 
latter value is approached only logarithmically. 



4. Improving the accuracy: The Z-averaging 

For more complex quantum impurity models, like the 
two-channel Kondo model (discussed briefly in the pre- 
vious section and in more detail in Sec. HV.Bp . or multi- 
orbital models, the Hilbert space per NRG step increases 
more strongly than for single-channel models. Conse- 
quently, the fraction of states kept in the truncation pro- 
cedure has to be reduced. As has been pointed out by 
lOliveira and Oliveiral (1994), this leads to an exponential 
decrease of accuracy, which can however be compensated 
by an increase of the discretization parameter A. How- 
ever, the use of a large A, (i), takes one further away 



from the continuum limit A — > 1 of interest, and, (ii), 
introduces oscillations into the thermodynamic expecta- 
tion values. A way out of this dilemma, proposed by 
lOliveira and Oliveiral ( 19941 ) . is as follows: 



• Instead of the discretization Eq. ^ choose 

x n = \- n+z , n > 1 , Z G [0, 1) . (61) 

The mapping to a semi-infinite chain is done as 
before (see Sec. III. CI) with the replacement hr n — > 

A -n+Z for n > 1 

• For fixed Z S [0, 1) perform a NRG calculation for 
a fixed set of temperatures Tl = A~( i_1 )/ 2 //3 as 
before. 

• Average over several calculations for different Z. 
This averaging is meant to r eintroduce the contin- 
uum limit to some extent ( Oliveira and Oliveiral 
Il994h and also can be shown to remove oscillations 
introduced by the use of a large A> 1. 

Already for two different values of Z this procedure re- 
moves spurious oscillations in thermodynamic quanti- 
ties and reproduces the exact result with good accu- 
racy for A as large as A = 10. This technique can be 
incorporated str aightforwardly into the NRG code (fo r 
applications , see Campo. Jr. and Oliveiral (120031. 2004): 



ICosta et al\ dl997ft7 iPaula et all (Il999fk iRamos et al\ 
H (|1996I 



(l2003r ): ISilva et all fjl99&\ ) 



B. Dynamic properties 

1. Equilibrium dynamics and transport 

We consider now the application of the NRG to the cal- 
culation of dyna mic and transport propertie s of quantum 
impur i ty models dCosti and Hewson L 1992ri|Costi et all , 
Il994at iFrota and Oliveiral , Il986t ISakai et ad . Il989h . For 
definiteness we shall consider the Anderson impurity 
model and illustrate the procedure for the impurity spec- 
tral density A a (oj,T) = -±lmG a (w,T), with 



dit-t'y^-^G^t-t') ,(62) 



G a (t - O = -iQ{t - t')([Ut)Ji{t')] + ) e , (63) 

with g the density matrix of the system. Suppose, for the 
moment, that we have all the many-body eigenstates, |r), 
and eigenvalues, E r , of the Anderson impurity Hamilto- 
nian, H, exactly. Then the density matrix, g{T), and 
partition function, Z(T), of the full system at tempera- 
ture kftT —1/(3 can be written 



Q{T) 
Z(T) 



-l3E r 



r)(r\ 



(64) 
(65) 
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and the impurity spectral density, A a , can be written in 
the Lchmann representation as 

M»>T) = ^yE|M„,| 2 (e-^ /fcBT + e-^/ fc ^) 



xS(uj - (E r > - E r )). 



(66) 



with M r y — (r\f a \r') the relevant many-body matrix 
elements. 

Consider first the T ~ case (T > is described 
below), then 

A a (u,T = 0) = E \ M r,o\ 2 H" + (E r - E )) 



(67) 



with Eq — the ground state energy. In order to eval- 
uate this from the information which we actually obtain 
from an iterative diagonalization of H, we consider the 
impurity spectral densities corresponding to the sequence 
of Hamiltonians Hn, N — 0, 1, . . ., whose characteristic 
scale is con = |(1 + A^A-^" 1 )/ 2 , 



A-a (w, T = 0) = £ \M? >Q \ 2 5{u + E, 



(68) 



Here, £f and |r) a? are the eigenvalues and eigenstates 
of H n ? i-e. 



iV, 



and, 



N 



N, 



(69) 
(70) 



are the relevant many-body matrix elements, whose cal- 
culation will be outlined below. Since the spectrum of 
Hn is truncated, the range of excitations it describes 
is limited to < u> < K(A)ujn, where K(A) depends 
on both A and the actual number of states retained at 
each iteration and is typically 5 — 10 for A = 1.5 — 2.0 
for N s — 500 — 1000 retained states. Moreover, excita- 
tions and eigenstates below the characteristic scale u>n of 
Hn will only be approximations to the excitations and 
eigenstates of the infinite system described by H. These 
excitations and eigenstates are refined in subsequent it- 
erations. Hence, for each N = 1, 2, . . . we can evaluate 
the spectral density from at a frequency uj chosen to 
lie in the window lon < u> < K(A)lon, 



A a (u>,T = 0)ttA»(u>,T = 0) 



(71) 



A typical choice, for A = 1.5 — 2.0, is uj = 2ujn- 

The above procedure still only yields discrete spectra. 
For comparison with experiment, smooth spectra are re- 
quired, so we replace the delta functions 5{uj ± E^) ap- 
pearing in (|6"8)) by smooth distributions P(ui ± E^). A 
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FIG. 7 The spin-resolved Kondo resonance at high magnetic 
fields calculated with NRG and perturbative RG. The large 
Gaussian broadening used at ui = <?/ib£> = 100Tk reduces the 
height of the sharp peak and overestimates its width, as is 
evident on applying the NRG b roadening procedu re to the 
analytic perturbative RG result dRosch et aU 12003). 



natural choice for the width rj^ of P is ujn, the character- 
istic scale for the energy level structure of H^- Two com- 
monly used choices for P are the Gaussian, Pg, and the 
Logarithmic Gaussian, Pr,a, distributions ( Bulla et all 
120011: ICosti ' et al. |, Il994at ISakai et all Il989h : 



Pg(u±R 



p lg {w±e 1 ;) = 



i 



- b 2/4 ln(^/E~) 



bE?^ 



(72) 



, (73) 



For the G aussian, a width tin = 0.3w7v — 0.8wjv is typi- 
cally used ( Costi et~aZl . ll994ah . whereas, for the logarith- 
mic Ga ussian, a typical width parameter b = 0.3 — 0.7 
is used (|Bulla et aZ.I . l200lfc ISakai et a/.l . ll989t) . Note that 



the logarithmic Gaussian gives little weight to excitations 
below lun and more weight to the higher energy excita- 
tions. Due to the logarithmic discretization, this might 
appear to be the better choice. In practice, the difference 
to using a Gaussian is small. 

In general, spectra for even and odd N differ by a 
few % at most as a result of finite-size effects (see also 
the discussion of the reduced density matrix approach 
in Sec. IIII.B.2[) . so generally either even N or odd N 
spectra are calculated (as for thermodynamics). It is 
also possible to combine informatio n from shell 7V an d 
N + 2 by an appropriate weig hting ({Bulla et all l200lh . 
We note also, that since the broadening is proportional to 
energy, a peak of intrinsic width Y at frequency Qq will 
be well resolved by the above procedure provided that 
Hq <C r, which is the case for the Kondo resonance and 
other low energy resonances. In the opposite case, the 
low (logarithmic) resolution at higher frequencies may 
be insufficient to resolve the intrinsic widths and heights 
of such peaks, although their weights are correctly cap- 
tured. In cases where the width of such high energy 
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peaks is due to single-particle effects, e.g. the resonant 
level in the empty orbital regime of the Anderson model, 
one can use the representation of the spectral density 
in terms of the correlation self-energy, as described in 
the next section, with the single-particle broadening be- 
ing put in explicitly so that essentially the correct peak 
widths and heights is obtained. In other cases, when 
the width of such peaks is due to correlations, one in- 
evitably has some over-broadening. An extreme example 
is the spin-resolved Kondo resonance at high magnetic 
fields, B > Tk, which is sharply peaked at u> — B and 
is highly asymmetrical, as shown in Fig. [7J The extent 
of the problem is quantified here by comparison with an- 
alytic perturbative results with and without the NRG 
broadening procedure. 

A procedure for obtaining smooth spectra, which re- 
solves finite frequency peaks without broadening the dis- 
crete spectra, involves an averaging over many different 
discretizations of the band (the Z-averaging discussed in 
the previ ous section on therm odynamics). We refer the 
reader to lYoshida et al\ (|l990h for details. 

In calculating the impurity spectral density, one re- 
quires also the matrix elements M^ r , at each iteration. 
These are obtained recursively by using the unitary trans- 
formation Eq. ((33]) yielding 



M r,r> = E E U N (r,ps N )U N (r ,p's' N ) 



P,SN p',sL 



N-l 



(74) 



Hence, the matrix elements M^ r , can be evaluated re- 
cursively from a knowledge of the eigenstates of all fi- 
nite size Hamiltonians up to Hn starting from the initial 
matrix elements -\{r\f a 

Himp = £f X/er fa fa 



_i of the isolated impurity 
C//|/|/|/|. Similar considerations 
apply to other local dynamical quantities such as dy- 
namical spin and charge susceptibilities. Figure [5] shows 
T = spectral densities for single-particle, magnetic and 
charge excitations calculated using the above procedure. 
These NRG calculations have been shown to satisfy exact 
Fermi liquid relations, such as the Friedel sum rule for the 
single-particle spectral density and the Shiba relation for 
the magnetic excitation spectrum, to within a few per- 
cent irrespective of the interaction strength U/irA in the 
Anderson model o r the value of the excha nge J in the 
Kondo model, see ICosti et all (|l994al ) and ICostil (|l998l ) 
for a discussion. 

The case of finite temperature dynamics is more com- 
plicated. Contributions to the spectral density at fre- 
quency uj ~ lun now arise from excitations between arbi- 



Eq = E r 



En 



trary excited states, i.e. uj = E r 
E r n - E p " = . . . with E = < Epf < E v „ < . . . 
see Fig. |H1 Consequently, the finite-T spectral density 
at uj ~ lun will have contributions from all energy shells 
n = 1, . . . , N. These need to be summed up, as in the cal- 
culation of transient quantities described in Sec. IIII.B.3I 
It is clear, however, that in the case of equilibrium spec- 
tral densities, the contributions from shells n < N will 
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FIG. 8 T = spectral densities for single-particle (solid line), 
magnetic (dashed line) and charge (dot-dashed line) excita- 
tions of the spin degenerate symmetric Anderson model versus 
energy to/D for U = 0.6D, D = 1.0 and A/vr = M = 0.03D 
(|Sakai et al.l,ll989h . 



be suppressed by Boltzmann factors. This motivates the 
following approximation: at u) = 2uj n > k^T one can 
calculate A a (oj, T) as in the T = case 

4(«n,r)«^K,T) 



Z n (T) ^ 

xS(oj n -(E?,-E?)). 



(75) 



In the other limit, uj — 2uj n < T, there is no com- 
pletely satisfactory procedure. One approach assumes 
that the main contribution to the spectral density for 
uj = 2uj n < T comes from the energy window containing 
thermal excitation s O(ksT) (|Costi and Hewsonl Il992bl: 
ICosti et aLl . ll994ah . In this case, the relevant shell, M, 
is determined by temperature via ujm ~ Pk^T, as in the 
evaluation of thermodynamic properties in the previous 
section, so that, for uj = 2uj n < /cbT, we use 

A a {uj n ,T)^ Af K,T) 

= Y^Tn E \M r A U 2 (e- E " /kBT + e-^l^ T ) 

&M{-L ) 



x6(uj n -(E™-EM)) 



(76) 



In practice, this procedure gives a smooth crossover as 
uj is lowered below k^T for temperatures comparable to 
the Kondo scale and higher, but becomes less reliable at 
uj < k B T < T K . 

Once the finite-T spectral density is known, one can 
also calculate transport properties, since the transport 
time, Tt r , for electrons scattering from a small concentra- 
tion, rii , of magnetic impurities is given in terms of the 
spectral density by 



— = -±AA a (uj,T), 
T tI (uj,T) N F 



(77) 
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FIG. 9 Excitations of Hn contributing to the spectral func- 
tion at frequency uj for, (a), T = 0, and, (b), T > 0. 



where Np is the conduction electron density of states 
and A is the hybridization strength. For example, the 
resistivity R(T) due to Kondo impurities in a clean metal 
is given by 



R(T) 



1 



(78) 



and the conductance through a quantum dot (or the re- 
sistivity of Kondo impurities in a dirty metal) modeled 
by an Anderson impurity model is given by 



(79) 



Figure ITD1 compares the scaled resistivity R(T)/R(0) for 
the Kondo and Anderson models with the scaled conduc- 
tance G(T)/G(0) for the Kondo model. The conductance 
and resistivity are seen to be almost identical universal 
functions of T/Tk. At finite magnetic field, the two quan- 
tities deviate from each other in the region T w B (|Costil . 
2000). The NRG results can be compared to analytic re- 
sults at low and high temperature. The resistivity of 
the Anderson model in the low- temperature Fe rmi liquid 
regime is given by the result of iNozierei ( 1974) . 



/?(7> /?(.M = (-V(7) •<-;(!))= I -< 



,T«T K , 



(80) 

where c = 7r 4 /16 = 6.088 and Tk is the low-temperature 
Kondo scale defined from the static spin susceptibility 
via 



X (T = 0) = ( gf , B ) 2 /4k B T K . 



(81) 



At high temperatures, T > Tk, H a mann used the 
Nagaoka-Suhl approximation ( Hewsonl . [T993) to obtain 
for the resistivity of the Kondo model 



R(T)/R(0) - 1(1 - {HT/TK ^ T ll2sl s + w /2 



), 
(82) 
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FIG. 10 Scaled resistivity and condu c tance of the S — 1/2 
Kondo model. Adapted from ICostil |2000). For compar- 
ison the resistivity of the symmetric Anderson m odel for 
U/nA — 4 is also shown (|Costi and Hewsonl . Il993f) and is 
seen to be identical to that for the Kondo model, up to 
non-universal corrections arising from charge fluctuations at 
higher temperatures (for U/nA = 4 these corrections occur 
for T > IOTk). 



where S is the impurity spin and Tkh is a Kondo scale 
defined by 



R(T = T] 



KH 



i?(0)/2. 



(83) 



iMicklitz et all (|2006a^ found numerically that Tkh ~ 
0.91T K . We see from Fig. [TO] that the NRG result for the 
resistivity of the Kondo model agrees with the Hamann 
result for T > Tk. The T 2 Fermi liquid behavior at 
low-temperature T <C Tk is also recovered. In contrast, 
the Hamann result violates Fermi liquid behavior and 
cannot be trusted for T < Tk- A numerical determi- 
nation of the coefficient c in Eq. (|80p requires obtaining 
Tt r (uj,T) accurately up to second order in both u> and T 
( Costi et al. L ll994ah . Typical errors for c can be as large 
as 10-30% so there is room for further improvement of 
the finite T dynamics in the Fermi liquid regime T <C Tk- 
For a discussion of other transport properties of Kondo 
systems, such as thermopow er and thermal conductiv - 
ity, the reader is r e ferred to ICosti and Hewsonl (|l993T ); 
ICosti et all (|l994al) ; Iziatic et all (|l993l ). 



2. Self-energy and reduced density matrix approach 

We now describe two improvements to the calculation 
of dynamical quantities. The first of these, a direct cal- 
culation of the correlation part of the s elf-en ergy of the 
Anderson impurity model (jBulla etUll , Il998h , is partic- 
ularly important for applications to DMFT, where the 
impurity self-energy plays a central role. The second, 
the introduction of the reduced density matrix into the 
calculation of dynamics, is important, for example, in 
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correcting large finite-size errors in spin-resolved spectra 
of the Anderson and Kondo models when a m agnetic field 
perturbs the ground state ( Hofstetterl . l2000h . 

The correlation part of the self-energy for the Anderson 
impurity model, S^, is defined via 



uj — Ef + iA — T) ' 



(84) 



and can be express ed, via the equation of motion for G a 
( Bulla et all Il998f ), as the ratio of a two-particle and a 
one-particle Green function 



Mu,T) = U((Uf_J-y,ft))/((U;fi)) 



(85) 



Evaluating the spectral densities of the two Green func- 
tions in (|85p as in the previous section, and calculating 
from these, via a Kramers-Kronig transformation, the 
corresponding real parts of the Green functions one ob- 
tains the self-energy. Using this in l]84p one is able to 
obtain the impurity spectral density with improved reso- 
lution of high-energy peaks, since in this procedure, the 
single-particle broadening A is included exactly. In par- 
ticular, this scheme recovers the limit U — > exactly. 
It is also found to improve the spectral sum rule with 
typical errors as low as 0.1% or less. 

The evaluation of spectral densities described in the 
previous section is subject to systematic errors due to 
neglect of high-energy states in constructing Hn ■ These 
are the same as for thermodynamic properties, and they 
can be controlled by increasing A and the number of re- 
tained states N s . Another source of error, specific to 
the method used to calculate dynamics, is that while we 
chose the frequency u in evaluating spectra from Hn 
carefully so that to > lun, nevertheless the eigenstates in 
the range < < ojn, which for small N are only 
crude approximations to the eigenstates of H, are also 
used in the evaluation. They enter the calculation di- 
rectly, as can be seen from Fig. [HI and, also, via the 
density matrices (e.g. at T — via qn = |0)jviv(0|) 
which are used to arrive at (|68|) and (|75ti76[) . As a result, 
the spectral density is subject to errors for small N, i.e. 
for high energies, due to the use of l ow- lying states which 
are not converged ()Hofstetteil . [2000h . With increasing A, 
i.e. lower ene rgies, this error w ill decrease. An improve- 
ment, due to iHofstetterl (|2000h . is to use in place of qn 
the reduced density matrix, g N ' Icd , of Hjy, obtained from 
the density matrix of the largest finite size Hamiltonian 
diagonalized, ffjv mal! i-e., 



JV.rod 



[Q 



(86) 



where s 



N+l- ' 



are the extra degrees of freedom 



contained in ffjv majt but absent in H^- As g"> red is not 
diagonal in the eigenbasis of iZjVj the resulting spectral 
function takes on a more complicated form 

A»(u,0) = *£c r ,, r M» r ,5(u;-(E»-E?)) (87) 

r' ,r 

C r ',r = ^ Py^ CAM r',y + Pr'fp M p,r- ( 88 ) 




FIG. 11 Comparison of the spectral function for the Ander- 
son impurity model calculated in a magnetic field with and 
without the reduced d ensity matr i x: A = 0.01, U = 0.1, e/ = 
-0.05, B = 0.001 from lHofstetteil ^OOOD . 



The reduced density matrices g n - red are calculated iter- 
atively backwards starting from the density matrix of 
_Hjv max . One situation where the reduced density matrix 
is important is in obtaining correctly the spin-resolved 
spect ral density of the Anderson model in a magnetic 
field (IHofstetterl . |2000). A magnetic field comparable 
to 7k changes the magnetization and therefore the oc- 
cupation of up/down states by 0(1) so large shifts in 
spectral weight occur at high energies in the impurity 
spectral density which are only captured correctly by us- 
ing the reduced density matrix, see Fig. 1111 Results for 
dynamical susceptibilities in a magnetic field u sing the 
reduce d density matrix have been discussed bv iHewsonl 
(|2006t ). The reduced density matrix also eliminates to 
a large extent the difference between the spectra calcu- 
lated for even/odd N and allows a correct description of 
the as ymptotics of the Kondo resonance in high magnetic 
fields (|Rosch et aLl . l2003l ). 



3. The x-ray problem and transient dynamics 

We consider here the calculation of the response of a 
system to a sudden local perturbation, such as the ex- 
citation of an electron from a core-level to the conduc- 
tion band of a metal in the x-ray problem, or, the time- 
dependent response of a spin in the spin-boson model fol- 
lowing an initial-state preparation. The NRG app roach 
to the x-ray problem (Ol iveira and Wilkinsl . Il98lh was 
the first application of the method to dynamical quanti- 
ties. In common with pre vious treatmen ts of the x-ray 
problem, for a review see (lMaha4ll975l) . the approach 
developed bv lOliveira and W ilkins ( 1981 ) calculates the 
absorption spectrum within linear response theory, and, 
therefore, belongs logically to Sec. IIII.B.ll We include 
it here for two reasons, (i), because the response of the 
electrons to the appearance of the core-hole potential is a 
problem of transient dynamics, and, (ii), because it uses 
the idea of formulating the calculation of the absorption 
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spec trum in terms of initial and fina l state Hamiltoni- 
ans ( Nozieres and De Dominicisl . Il969l ) , which is also in- 
herent to the recent NRG approach to transient dynam- 
ics beyond linear r esponse theory ( Anders and Schillerl . 
120051: ICostiLll997ah . 

A simple model for describing the x-ray absorption 
spectra in metals, is given by the following spinless 
Hamiltonian 



H 



k 



E d dU ■ 



k,k' 



U dc c\c k ,dd\ (89) 



where d) creates a core-electron with energy E d and the 
attractive screening interaction, U dc , acts only when the 
core-level is empty {dd^ = 1). The core-level lifetime is 
assumed infinite, and the interaction with the x-ray field 
is taken to be of the form 



wic^de 



+ H.c. 



(90) 



where cq = J2 k Cfc - The x-ray absorption spectrum, h{uj), 
is obtained using linear response theory from the imagi- 
nary part of the optical conductivity, \cd — ((cod; d^co)). 
At zero temperature, one finds for the absorption spec- 
trum a power law singularity of the form 



n(w) ~ (to — Et) 



(91) 



where, Et is the absorption threshold, and, a, is an expo- 
nent which depends on the strength of the core-hole po- 
tential. The exponent a has t wo contr i bution s a — a' — e, 
an excitonic part, a', due to iMahanl (|l96?t ) and an or- 
thogonality part, e, which follow s from Anderson 's or- 
thogonality catastrophe theorem ( Anders on. 1967]). An 
ex act solution of the 'x-ray problem' has been obtained 
by INozieres and De Dominicii ( 19691 ) by using the de- 
composition of (|89|) into single-particle initial-state, if/, 
and final-state, Hp, Hamiltonians, corresponding to the 
situations before (defl = 0) and after (dd^ = 1) a core- 
electron is excited to the conduction band: 



Hi = ^ £ k c\c k + E d, 

k 

Hp = e k c\c k + ^ Ud c c\c k , 



(92) 
(93) 



k, k ' 



For th e spinless model, (|89")1 , INozieres and De Dominicisl 
( 19691 ) found the exponents 



a' = 25/ir, e = (S/n) 2 , 



(94) 



where the phase-shift S = arctan(— irNpUdc) is that for 
conduction electrons scattering from the additional po- 
tential created by the core-hole, and, Np is the conduc- 
tion electron density of states at the Fermi level. In ad- 
dition to the absorption spectrum, the core-level pho- 
toemission spectrum, Ad(uj) = —Im[((d;d^))]/TT, is also 
of interest. In the core- level photoemission spectrum, 
only the orthogonality effect is operative and the core- 
electron spectral function, which has the quasiparticle 
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FIG. 12 Absorption spectrum normalized by the Nozieres- 
De Domenicis result (Eq. (|9ip ) for sever al screening potentials 
of the form U dc (k,k') = G +Gi(k+k') l|Qliveira and Wilkinsl . 
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form Ad{oj) = 5(lu — Ed) in the absence of screening, is 
replaced by an incoherent spectrum of the form 



A d (u) ~ 9(E d - co)(E d - lu)- {1 -^ 



(95) 



in the presence of scree n ing dDoniach and Suniid Il970l ) . 

lOliveira and Wilkinsl (|l98ll ) applied the NRG to the 
initial-state and final-state Hamiltonians (|92j|93p and cal- 
culated the zero temperature linear response absorption 
spectrum, 



= 2nw 2 ^ \{mp\c\d\m^Gs)Y 

1 11 p 

xS(uj — (Em F — E mi GS )), 



(96) 



with Hj t p\mi_p) = E mi F \mi^p) and \mj^Gs) the ground 
state of Hj. In evaluating truncated Hamiltoni- 

ans, Hf F , were used and the spectrum was evaluated 
at an appropriate frequency ui = lun as in Sec. MLB. II 
but with a box broadening function on a logarithmic 
scale. They w ere able to recover the exact th resh- 
old exponent of Nozieres and De Dominicisl ( 19691 ) (see 
iLibero and Oliveiral ( 1990a ) for similar calculations of 
the photoemission spectra). They also extended the cal- 
culation of absorption spectra to core-hole potentials of 
finite range Udc — > U dc (k,k'), finding (see Fig. [T2]) that 
the threshold exponent remains universal, i.e. it depends 
only on the phase shift 8 at the Fermi level, but that 
the asymptotic scale for the onset of power-law behavior 



depends on both the range (|Qliveira and Wilk ins. 1981) 
and the strength of the core-hole potential (|Cox et 



119851) . It also de pends on any ener gy dependence in the 
density of states (jChen et~a/lll995l ). These results reflect 
the fact that the crossover scale to the low energy fixed 
point, determining universal properties, depends on de- 
tails of the density of states and the core-hole potential. 

X-ray singularities also play an important role in 
the dynamics of auxiliary particles in the slave-boson 
approach to the infinite U Anderson impurity model 
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(|Colemaij . 1 1984 iMuller-Hartmannl . 1 1984) . A NRG cal- 
culation of the T = photoemission spectra for slave- 
bosons, A;, (a;), and pseudo-fermions, A/o^w), showed 
that these diverge with the exponents given above for 
the photoemission and abs orption spec t ra re s pective ly, 
generalized to include spin, (jCosti et aLl . ll99fiLll994bh : 



with 



-atf 

5 



A frT (to) ~ (u-Er) 
A b {uj) ~ (uj-Et) 

7T z — ' 7T 

« 6 = 1-E(-) 2 ' 



(97) 
(98) 

(99) 
(100) 



with the phase shift 5 a = ixnf a and nf a the occupancy 
per spin of the local level. Here, also, the scale for the 
onset of power-law behavior is determined by the relevant 
low energy crossover scale, e.g., the Kondo scale in the 
Kondo regime. 

We turn now to a problem which is formally similar to 
the x-ray problem, namely the dynamics of a spin subject 
to an initial state preparation, as in the case, for exam- 
ple, of the dynamics of the spin in the spin-boson model 
Eq. p39[) . Further discussion of the effects of screening 
on the spectra of impurity models is given in Sec. lIV.A.ll 
In the spin-boson model one is interested in the dynam- 
ics of a spin, <r z , described by P(t) = (cr z (t)) ei , following 
an initial-state preparation o f the system describ ed by 
an initial density matrix gi (jLeggett et all Il987t ). For 
example, the spin a z in (|139[) could be prepared in state 
X) at t < by an infinite bias e = oo with the environ- 
ment fully relaxed about this state, and the bias could 
subsequently be switched off at t — allowing the spin to 
evolve. The time evolution of the spin for t > is then 
described by 



i 



Trg_ 



■Tr[ gi e- lHFt/h a z e 



iH F t/hl 



(101) 



where gj = Tr[e P Hl ] and the initial and final state 
Hamiltonians are given by 



Hi = H SB (e = oo) 
H F = H SB (e = 0) 



(102) 
(103) 



where Hqb is the spin-boson Hamiltonian (|139|) . 

This appro ach has been investigated within NRG by 
ICostil ((l997a) using, for the Ohmic case, in place of -ffsB, 
the equivalent anisotropic Kondo model. Despite the for- 
mal similarity with the x-ray problem, the exact formula- 
tion indicates the difficulties that have to be overcome in 
calculating transient dynamical quantities beyond linear 
response. Consider the spectral function 



P(u) 



1 



E 



-/3E„ 



' (mi\mp) (m' F \mi) 



x (m F \a z \m' F )5{uj - (E mF - E m <)), (104) 



P(t) = / 

Jo 



P{lo) cos(ujt)duj. 



(105) 



We see that even at T = 0, no ground state enters the 
delta functions in (|104[) , in contrast to the linear-response 
expression (|96|1 for h(uj) in the x-ray problem. This im- 
plies that in evaluating P{u>) at a frequency u> ~ u>n, con- 
tributions will arise from all energy shells n — 1,2, ... ,N 
as discussed previously for finite-T dynamics. In the 
present situation, however, the contributions from higher 
energy shells (i.e. n < N) are not suppressed by Boltz- 
mann factors, so it is not clear a priori that using a single 
shell approximation will give meaningful results. Such 
an approximation shows that the short-time dynamics of 
the spin-boson model can be recovered and that in or- 
der to obtain the long-time d ynamics one h as to sum up 
contributions from all shells ( Costil . [l997al ). Adding up 
such contributions using the retained states of succes- 
sive Hamiltonians Hfp, Hfp 1 , ... is problematical due 
to the overlap of the spectra at low energies. An elegant 
solution of this problem, allowing multiple-shell NRG cal- 
culations to be carr i ed ou t, has recently been found by 
lAnders and Schillerl (|2005l ). Their idea, was to recognize 
that the set of states discarded, |r) m)( jisc) at each NRG 
iteration m, supplemented with the degrees of freedom 



»m+i7 



I s N„ 



for m = 1, 



with N max the largest Hamiltonian diagonalized, forms a 
complete basis set 



Nmax 

E 

m=l rS{disc} 



E \r,e;m)(r,e;m\ 



= 1 



(106) 



Using this identity, their result for P{t) in the basis of 
final states is 



N trun 

p{t) - E E e l(K 

in — 1 mp ,m' F 
/ i i / \ m.red 

x(m F \a z \m F )g m • , , 



-E m , )t/h 



(107) 



where Q™/ 1 ^ are the matrix elements of the reduced den- 
sity matrix qi for HJ 1 introduced in the previous sec- 
tion and the sum Y^mp m' F implies that at least one 
of the states mp,m! F is in the discarded sector for it- 
eration to. Rotating g^f^ F to the initial state ba- 
sis gives overlap matrix elements (mi\mp), (mi\m! F ) 
as in (1104[) above. Within this approach the time- 
dependent transient dynamics of a number of models 
has been investig ated, including the Anderso n and reso- 
nant lev el models (I Anders and Schillerl . l2005f ). the Kondo 
model ( Anders and Schillerl |2006 | ) an d the sub-ohmic 
spin-boson model ([Anders et aUfeOOd ). 

The use of a complete basis set (|106l) has poten- 
tial to improve also the finite-T calculation of spec- 
tral densi ties, particularly in the pro b lemat ical range 
lu < T (jWeichselbaum and von Delftl . l2006h . A fur- 
ther improvement in using a complete basis set is that 
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the sum rule for spectral densities J ™ dcuA a (uj, T ) = 1 
is, by construction, fulfilled exactly (jPeters et al .1. 1200fil : 
IWeichselbaum and von Delfd . l2006h . 



IV. APPLICATION TO IMPURITY MODELS 

In this section we review the application of the NRG 
to a range of quantum impurity models. Section IIV.AI 
reviews work on models with conduction electron screen- 
ing (Sec. lIV.A.ip . underscreened and fully screened 
Kondo models (Sec. IIV.A.2j) . and models exhibiting the 
Kondo effect in nanostructured devices (Sec. IIV.A.3[) . 
Section IIV.BI deals with the prototypical overscreened 
Kondo model, the two-channel Kondo model, which 
is often encountered as an effective model describing 
the quantum criticial point of more complex quan- 
tum imp urity models, e .g. ce rtain double quantum dot 
systems ( Zarand et all, 120061) or two- impurity systems 
([Jones and Varmal . ll987r ). Impurity quantum phase tran- 
sitions are reviewed in Sec. IIV. Cl in the contexts of multi- 
impurity systems (|IV.C.1[) , soft-gap systems pV.C.2|) and 
in the context of magnetic impurities in superconductors 
PV.C.3[) . Sec. IIV.DI reviews work on multi-orbital sys- 
tems, including the effects of crystal-field splittings in 
Anderson impurity models. Finally, models with bosonic 
degrees of freedom are reviewed in Sec. IIV. El Note 
also, that a number of models of nanostructured devices, 
for example, the single-electron box, quantum dots with 
phonons, or multi-orbital quantum dots, are to be found 
in Sec. EHS1 



A. Kondo effect and related phenomena 

1. Screening and photoemission 

Screening effects are important whenever an electron 
is excited from a localized core- or valence-state into the 
conduction band or removed completely, leaving behind a 
hole which attracts the conduction electrons. Such effects 
can have a drastic influence on the photoemission and 
absorption spectra of impurity systems. In this section 
we consider a number of extensions to the basic screening 
model ([M]) introduced in Sec. IIII.B.3I 



H = £kC l Ck + EddU + H U dc c\c k> dS. (108) 



k.k' 



We consider first the generalization of (|108|) to a sim- 
plified model of an atom ad sorbed on a metallic surface 
(|01iveira and Wilkinsl . [l985[) . ^n addition to a core-level, 
as in (|108p . the atom has a resonant level (created by tf) 
whose position, E b j or E b F, depends on the occupancy 
of the core-level according to the term 



and the level hybridizes with the conduction band via a 
term 



H iaix = Vy(clb + H.c). 



k 
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Note that Hdb just represents a screening of the core- 
hole by electrons in the resonant level. Excitation of 
an electron from the core-level by an x-ray can pro- 
ceed either directly or via the resonant level, leading 
to a Fano anti-resonance in the x-ray absorption at fi- 
nite energy (in addition to the usual edge singularity at 
u> = Et)- This Fano anti-resonance, present also with- 
out the core-hole potential, is found to be significantly 
narrowed and shifted in the p r esence of the core-hole 
potential I Oliveira and Wilkinsl . Il985t) . It would be of 
interest to investigate also the core-level photoemission 
spectrum of this model using the NRG, as both this, and 
the absorption spectru m, are accessible in experiments. 
iBrito and Frotel {l990) have carried out such a calcula- 
tion for an appropriately generalized spinfull version of 
the above model, i.e., the Anderson model ([2]) in the 
presence of both an interaction, 



Hdc = Udcdd) ^2 C ka C k'<J, 



(111) 



kk' 



between the core-hole and the conduction electrons, and, 
an interaction, 



(112) 



H db = E bI tfbdU + E bF tfb(l - dU), 



(109) 



between the core-hole and the valence level. Signatures 
of the valence states could be identified in the XPS spec- 
tra and their dependence on Udc was investigated in the 
mixed valence and empty orbital regimes. A reduction 
of the hybridization between the valence level and the 
conduction band, arising from orthogonality effects, was 
found with increasing Udc- The corresponding calcula- 
tion in the Kondo regime is still lacking. The x-ray ab- 
sorption spectrum of the sam e model has been investi- 
gated in (Hol mes et al. I. I2005D in the context of excitons 
in Kondo correlated quantum dots and the expected ab- 
sorption exponent from the Nozieres - De Domenicis the- 
ory was recovered. 

In the models discussed so far, the core-level was as- 
sumed to have infinite lifetime. Consequently, the screen- 
ing interaction gave rise to true singularities in the core- 
level absorption and photoemission spectra. These sin- 
gularities are cut off as soon as the core-level lifetime 
is finite. Another situation where the singularities due 
to screening are cut off, but where screening effects may 
nevertheless be important, is in the valence band photoe- 
mission spectra of heavy fermions within a local impurity 
approach, which we now address. 

It is often assumed that the effects of conduction elec- 
tron screening on the / electron photoemission spectra of 
heavy fermions can be taken into account by renormaliz- 
ing the bare parameters of an effective Anderson impurity 
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model. However, this is not a priori clear as the screening 
interaction in these systems can be an appreciable frac- 
tion of the local Coulomb repulsion. One of the merits of 
the NRG, which allows such questions to be investigated, 
is that it can deal with all local Coulomb interactions on 
an equal footing and in a non-perturbative manner, and 
some examples of this have already been given above. 
For the partic le-hole symmetric Anderson mod el, it 
was shown in ICosti and Hewsonl ( 199l[ Il992al) that the 
effect of a screening term of the form 



H f c = u fc flu 



ka 



Ck>. 



(113) 
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on the valence band photoemission spectrum could be 
well accounted for by a renormalization of the bare pa- 
rameters of the Anderson model, both the local level po- 
sition, e = —17/2, and the hybridization. An excitonic 
like enhancement of the hybridization was found with in- 
creasing Uf c . Similar effects are reflected in the STM con- 
ductance of a ma gnetic adatom modeled by the screened 
Anderson model ( Cornaglia and Balseirol 20031 ). Calcu- 
lations for the inf inite U Ander s on impurity model, for 
therm odynamics ( Alascio et al. [l986t IZhuravlev et all 



120051 ). and, dynamics ( Takavama and Sakail . Il993l ). are 



consistent with the above findings. 

The above model for screening in heavy fermions as- 
sumes that the largest contribution to screening of / holes 
arises from conduction electrons in states that hybridize 
with the / states. These hybridizing states are usually 
the p levels from neighboring ligand ions, so the screen- 
ing from these (denoted Uf c above) should be expected 
to be smaller than the onsite screening, Udf, from the d 
electrons of the rare-earth ion. By symmetry, the latter 
do not hybridize with the / states. Neglecting Uf c and 
representing the d electrons by a spinless s wave band, 
we may represent the screening of / holes by d electrons 
by adding to the Anderson model ^ the term 



Ha, 



J2e k d{d k 
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where nf = ^2 a ftfu- The resulting model is a two- 
channel Anderson model, in which one channel screens, 
but, does not hybridize, and, the other channel hy- 
bridizes, but, does not screen. Assu ming localized d 
electro ns gives the model studied by iBrito and Frotal 
(1990) and discussed above. For the full model, 
iTakavama and Sakail ( 19931 . fl997h calculated the valence 
band photoemission spectrum and, surprisingly, found 
that the effect of Udc in the Kondo regime could be ab- 
sorbed into a renormalization of the Anderson model pa- 
rameters. This result was for infinite U, but it should re- 
main valid in the Kondo regime for any finite U provided 
Udc remains smaller or comparable to U. In contrast to 
the model described above, where the screening occurs 
in the hybridizing channel, the effect of the screening in- 
teraction in the present model is to reduce the effective 
hybridization of the valence electrons to the conduction 



band, which can be understood as an orthogonality ef- 
fect. Nevertheless, we conclude from these NRG calcula- 
tions that in the Kondo regime, and with realistic values 
of Udf (Ufc), the valence band photoemission spectra of 
both the above screening models can be well accounted 
for by an Anderson model with renormalized parameters. 
The two-channel screenin g model above, ha s also 



been studied for f inite U ( Perakis and Varmal Il994t 
iPerakis et aLl .[l993). At particle- hole symmetry, increas- 
ing Udc reduces both U and the effective hybridization, 
resulting, for sufficiently large Udc, in an effective attrac- 
tive local Coulomb interaction and a charge Kondo ef- 
fect. For still larger Udc, a Kosterlit z-Thouless transition 
to a non-Fermi liquid sta te occurs ( Perakis and Varmal . 
11994 IPerakis et all , 1 19931) with a collapse of the Kondo 
resonance in the valence band photoemission spectrum 
(ICostil . ll997bh . 



2. Kondo effect in the bulk and underscreened models 

Real magnetic impurities in metals have both or- 
bital and spin degrees of freedom and the resulting low- 
energy ef f ective impurity models can be very complicated 
( Hewsonl . [l993). The NRG has been applied so far to 
models with at most three orbitals, see Sec. IIV.DI In 
cases where the ground sta te is an orbital singlet, e .g. for 
dilute Mn ions in metals, iNozieres and Blandi 3 (119801) 
have given a useful classification of the resulting effec- 
tive single-impurity Kondo models in terms of the size of 
the impurity spin 5 and the number of conduction chan- 
nels, n, which couple to the spin via the Kondo exchange. 
These multi-channel Kondo models are described by 



H 



^ £kC kaa 



C-kaoi 
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where a — 1, ...,n is the channel index, and the ex- 
change constant, J, is antiferromagnetic. For n = 25, 
complete screening of the impurity spin takes place lead- 
ing to a local Fermi liquid at low temperatures. The 
overscreened case, n > 25, exhibits non-Fermi liquid be- 
havior and is reviewed in Sec. IIV.BI In this section, we 
deal with some recent developments on the fully screened 
5=1/2 Kondo model, relevant to bulk Kondo impuri- 
ties, and also describe work on the the single-channel 
underscreened case, n = 1 < 25. 

One of the signatures of the Kondo effect is the appear- 
ance in the impurity spectral density of the Kondo reso- 
nance at the Fermi level. Point contact spectroscopy on 
Cu wires containing magnetic impurities, using the me- 
chanically controllable brake junction technique, show a 
zero bia s anomaly, wh i ch is a ttributed to the Kondo res- 
onance (jYanson et all Il995j ). In addition, these exper- 
iments show that the Kondo resonance splits in a mag- 
netic field. NRG calculations for the 5 = 1/2 Kondo 
model in a magnetic field do indeed show that the Kondo 
resonance splits in a magnetic field, B, provided t he Zee- 
man splitting gi/isB exceeds the Kondo scale Tk ( Costil . 
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,2000). Here gi,^B are the impurity g- factor and Bohr 
magneton and Tj< is the Kondo scale defined from the 
half-width at half-maximum of the T = Kondo reso- 
nance. The latter is obtained from the imaginary part of 
the many-body T-matrix, T kk / a , for spin a, defined by 



G kk > a (u>) - 6kk<G W(T + G kka T kk , a G k , k , c 



(116) 



where Gkk'o{w) = ((cfc CT ; %'<?)) is the full conduction elec- 
tron Green function and G kk , is the corresponding unper- 
turbed Green function. From the equations of motion for 
Gkk' one finds for the orbitally isotropic Kondo model 

W") = ^{S z ) + (^) 2 ({0„;Ol)), (117) 

(118) 



' to- C k<r , 



with t the Pauli matrices. From the T-matrix one can 
also extract the transport time and thereby the magne- 
toresistivity. The latter is found to agree well w ith ex- 
perim ental data on diluted Ce impurities in LaAl 2 ( Costil . 
l2000h . 

A recent develop ment has been the realization by 
IZarand et ali ([2004) that one can use the NRG to extract 
from the many-body T-matrix both elastic and inelastic 
scattering rates and cross sections. The total scattering 
cross section, <7 t ot(w), is related to the imaginary part of 
the T-matrix by the optical theorem 

2 

0tot(w = £ k) = Im[T kkrT (u!)], (119) 

Vk 

with Vk the velocity of electrons with wavevector k. Con- 
sequently, by using the expression for the elastic scatter- 
ing cross section 



<T e i(cJ = £k) = — ^ S(Sk> 



e k )\T kk >a(uj)f 



(120) 



IZarand et all ( 2004 ) were able to calculate the inelastic 
scattering cross section (7i ne i = Utot — fci an d the inelastic 
scattering time, T; nc i ~ Oj^ B i- I n order to shed some light 
on the expression for <7i ne i, consider the Anderson model 
for a flat band with density of states N F = J2k' H £ k'~ £ k) 
and resonant level width A = TrNpV 2 . We have, T kk > — 
V 2 G d , with G d = (u - e d + iA - S(w)) -1 and E the cor- 
relation part of the self-energy. Th e inelastic scatterin 
cross section for u = e k reduces to (jZarand et a*]. 1200 



Cinel = 



V 2 T,"{lo) 



v k (^-£ d -E'( w )) 2 + (A-E"( w )) 2 ' 



which shows that the inelastic scattering rate vanishes 
for electrons at the Fermi lev el due to the Fermi liquid 
properties of the self-energy. IZarand ~et Ml (|2004 eval- 
uated tJinoi for the S — 1/2 Kondo model via the NRG 
using the T-matrix in (| 1 1 Z|) at T = and for both zero 
and finite magnetic fields. The maximum in the inelastic 
scattering rate occurs close to the u> « Tk (see Fig. IT51) . 
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FIG. 13 Elastic, inelastic and total scattering rates for the 
S = 1/2 fu l ly scr eened Kondo model at T = (from 
IZarand et all (|2004h V The Kondo scale T K is that from the 
half-width at half-maximum of the Kondo resonance and is 
approximately twice that from the T — susceptibility de- 
fined in |[8T]). 
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FIG. 14 Universal dephasing rate for the S = 1/2 fully 
screened Kondo model calculated via NRG fo r Kondo impu- 
rities in d — 1, 2 dimensional conductors ffrom.Mickl itz et al\ 
(2006a)). The solid line in the inset is the analytic T 2 result 
from Fermi liquid theory valid for T <JC Tk, where Tk is the 
scale defined in (181 II . 



A quantity closely related to the inelastic scattering 
time, Tind ~ crjj, is the dephasing time, r^, for elec- 
trons scattering from magnetic impurities and measured 
i n weak-localization experiments on diffusive conductors 



( 121 ) (iBauerle et a/J. 120051 iMohantv et aR Il997t IPierre et ali 
et al\ ,\ 



2003: ISchopfer et all |2003) . The two quantities are, how- 
ever, not identical. An exact expression for the dephasing 
rate of electrons scattering from a dilute concentration 
of Kondo impurities in a weakly disordered metal has 
recently been derived (jMicklitz et all l2006al ). Fig. [Ml 
shows the dephasing rate as a universal function of T/Tk 
for the S = 1/2 Kondo model, obtained by using the 
NRG for finite temperature dynamics. The maximum 
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dephasing rate occurs at T w Tk and decreases at first 
linearly with temperature below Tk and eventually as 
T 2 in the Fermi liquid region T <C 7k- The magnetic 
field dependence of the dephasing time, t^B, T), has also 
been calculated and the expression for the dephasing rate 
has been generalized to arbitrary dynamical scatterers 
( Micklitz et aLl . l2006bt ). Recent experiments on Fe impu- 
rities in Ag wires show better than expected agreement 
with the theoretical predictio ns for the dephasing rate 
of the S = 1 /2 Ko ndo model (jAlzoubi and BirgeL Bool 



iMallet et all [2006) . Fe impurities in Ag will have both 
an orbital moment and a spin S = 2 in the absence of 
crystal field and spin-orbit interactions. Inclusion of the 
latter, may, however, result in an effective 5=1/2 single- 
channel Kondo model at the low temperatures T s» IK 
of the experiments, thereby helping to explain the good 
agreement with the S — 1/2 theory. At the very lowest 
temperatures, T < O.lTk, a slower decay of the dephas- 
ing rate has been reported in these experiments, as com- 
pared to that expected from a fully screened model. One 
possible explanation for this is that a small fraction of 
Fe impurities is only partially screened. Underscreened 
Kondo models, to which we now turn, are known to give a 
much slower decay o f the dephasing rate below the Kondo 
sc ale, see below and IVavilov and Glazmanl ( 20031 ). 

ICragg and Llovdl (11979T ) investigated the single- 
channel S = 1 underscreened Kondo model and showed 
that its low-energy fixed point corresponds to the spec- 
trum of the ferromagnetic S' — 1/2 Kondo model. The 
deviations from the fixed point at iteration N are of the 
form J(N)S' ■ s with J(N) being 0(-A/(N + C(J))) 
with A being a constant and C(J) depending on J, i.e. 
the deviations are margi nally irreleva n t. Th ese calcu- 
lations were extended by iKoller et al. I (j2005bl ) to mod- 
els with S — 1, . . . , 5/2. They also determined C(J) 
explicitly for the different cases. Using the relation of 
N to energy u> ~ Ar N / 2 , the effective coupling can be 
written as J( u>) ~ l/lnfy/Th) w ith To an appropriate 
Kondo scale (|Koller et all l2005bT ). Consequently, there 
are logarithmic corrections to thermodynamic quantities 
at low temperature, instead of the power-law corrections 
characteristic of fully screened Kondo models. Non- 
analyt ic corrections are al s o found in dynamic al quan- 
tities (|Koller et all l2005ri iMehta et all 120051) . so un- 
derscreened Ko ndo models have b een termed "singular 
Fermi-liquids" (|Mehta et all |2005() . For, example, the 
spectral density, pt{oS), obtained from the T-matrix (|117|) 
takes a finite value at the Fe rmi level, bu t the ap proach 
to this value is non-analytic (jKoller et all [2005b) : 



p t (u) = p t (0)-b/\n(uj/T o Y 



(122) 



Similarl y, the T = inelast ic cross section, also calcu- 
lated bv lKoller et all (|2005bl ). decays as l/ln(w/T ) 2 , at 
low energies, and consequently, the dephasing rate de- 
cays as t7 (T) ~ 1/ hi(T/Tb) 2 . As mentioned above, 
a small fraction of underscreened Fe impurities may ex- 
plain the excess dephasing obs erved at the l owest tem- 
peratures in the experiments of [Mallet et all ( 20061 ) and 



lAlzoubi and Birgel ( 20061 ). Calculations for the tem- 
perature dependence of the resistivity and dephasing 
rates of the spin S > 1/2 underscreened Kondo mod- 
els and t heir relevance to Fe impurities in Ag can be 
found in ([Mallet et all 120061 ). It is also interesting to 
note, that ca l culatio ns for ferromagnetic Kondo models 
(jKoller et all l2005hh show that all cross sections vanish 
at the Fermi level with the inelastic part contributing 
nearly all the scattering in this limit and with the elastic 
part being negligibly small. 

Finally we mention recent work on calculating spa- 
tial correlations such as spin-density correlations, C(x) = 
(S-s x ), around Kondo impurities, where s* x is the electron 
spin density at distance x from the impurity (for earlier 
wor k involving p erturbative aspects combined with NRG 
see (jChen et all Il987l 1 19921) ). iBordal (l2006f ) works with 
Wannier states centered at both the impurity and at x, 
thereby reducing the problem to a two- impurity type cal- 
culation (Sec. HV.C.ip . At T — and in one dimension, 
the decay of C{x) is found to change from 1/x to 1/x 2 
around x = £k = ^vp/T^, where the coherence length, 
£k> describes the size of the Kondo screening cloud. At fi- 
nite temperature, the expected exponential decay of C(x) 
for x > = hvp/ksT is recovered. 



3. Kondo effect in nanostructures 

Recent experimental work has demonstrated the 

importance of the Kondo effect in determining the 

low-tempcrature transport pro perties of nanosca l e size 



devices such as quantum d ots (ICronenwett et all 1998; 



Goldhaber-Gordon et all Il99a : van der Wiel et 
2000). An example of such a device, a quantum dot, is 
shown in Fig. [TS] More complicated devices, such as 
capacitively coupled double-dots or dots contained in one 
or two arms of an Aharanov-Bohm interferometer can be 
built up from this basic unit. A quantum dot consists of 
a confined region of electrons coupled to leads via tunnel 
barriers. It may be viewed as an artificial multi-electron 
atom, in which the different levels (filled, partially filled 
or empty) couple to electron reservoirs via one or more 
channels. A quantum dot can be described, in general, 
by the following multi-level Anderson impurity model 



H = -Hdot + H c + H t 



(123) 



H dot = J2 £ ^ d l^ + E c(N-My - J H S 2 



kacr 

-Htun = t a i\d icr C}-ac! + Cfcao-^tT 



where e; CT , i — 1, 2, . . ., are the dot level energies for spin a 
electrons, Af — (N) — Y^iaWedur) , is t ne dot occupancy, 
E c is the charging energy, S = \ J2i^u A^^ d ^ is tlie 
total spin of the dot and Jh > is the Hund's exchange 
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FIG. 15 Schematic top view of a lateral quantum dot, consist- 
ing of a confined region of typical size 10-100nm defined in the 
two-dimensional electron gas (2DEG) of a GaAs/AlGaAs het- 
erostructure. The dot is connected to left and right electron 
reservoirs. Gate voltages Vl,r control the tunnel barriers into 
and out of the dot, while Vg controls the dot level positions. 



coupling. In the above, a — L,R labels left/right lead 
degrees of freedom and k labels the wavevector of a single 
transverse channel propagating through the constriction 
between the 2DEG and the quantum dot. Electrons tun- 
nel into and out of the dot with amplitudes t a i and give 
rise to a single-particle broadening I\ of the levels. 

The model above is essentially the same model as the 
multi-orbital model of Sec. IIV.D1 used to describe bulk 
Kondo systems. The novel situation in quantum dots is 
that parameters such as the tunnel couplings and level 
positions can be controlled by gate voltages. This allows 
such models to be experimentally investigated in all phys- 
ically interesting regimes, such as spin and charge fluctu- 
ation regimes, and in principle also to be tuned through 
quantum phase transitions. In addition, different real- 
izations of quantum dots (nanotubes, vertical dots) may 
have level degeneracies or near level degeneracies, allow- 
ing the effects of Hund's exchange to be investigated. 
Finally, the devices described by (|123p can be driven 
out of equilibrium by a finite transport voltage, allowing 
the study of non-equilibrium effects in relatively "sim- 
ple" quantum many-body systems. This would be one 
motivation to further develop the NRG to steady-state 
non-equilibrium situations. 



a. Single-level quantum dots In the low-temperature 
limit, only one or two partially filled levels close to the 
Fermi level of the leads will be important for transport. 
The remaining levels will be either filled or empty, and, at 
the low temperatures of interest for quantum transport, 
they may be neglected. The simplest model, therefore, 
to describe low-temperature transport through a quan- 
tum dot, is the single- level Anderson impurity model ([2]) 
with level position E\ = £/ = —eVc controlled by gate 
voltage and Coulomb repulsion U given by the charging 
energy Ec — U/2. Only one conduction channel, the 



even combination of left and right electron states, aka, 
below, couples to the local level, as can be seen by using 
the canonical transformation 



aha = (tLCkLa +tRC kRa )/^Jt 2 L +t 2 R , (124) 

bka = {-tRCkLo + t L CkRa)/<Jt 2 L + t 2 R , (125) 



with t a — t a iSi\. We note that treating the 
Coulomb interaction classically implies that, for an in- 
teger number of electrons on the dot, transport is 
blocked for large U, since transferring electrons through 
the dot requires overcom i ng th e lar ge Coulomb repul - 
sion. iGlazman and Raikhl (| 19881) and lNg and Lee] (|l988l ) 
pointed out, however, that in the situation where the to- 
tal spin on the dot is finite, as happens for an odd number 
of electrons (i.e. for TV = 1 in the effective single-level 
model), one should expect, on the basis of @, an en- 
hancement of the conductance to its maximum possible 
value of G = 2e 2 /h via the Kondo effect in the limit of 
zero temperature. A device, represen ting a tunable An- 
derson impurity model, was realized ( Cronenwett et all . 
119981 : iGoldhaber-Gordon et all . 1 19981) and the predicted 
enhancement of the low-temperature conductance for 
dots with an odd number of electron s was measured and 
compared ( Goldhaber-Gordon et all. Il998h to quantita- 
tive NRG calculations ( Costi et aZj . ll994bj ) such as those 
shown in Fig.fTOlfor the conduct ance in the Kondo regime 
(see also llzumida et all ( 2001al )). Tuning the quantum 
dot to the mixed valence and empty orbital regimes, has 
enabled also comparisons with theory i n those regimes 
(ICostil . l2003 ISchoeller and Konid . 120001) . 

The frequency dependence of the linear conduc- 
tance, G'(ui), of a single-level quantum dot de- 
scribed by_J5^_Jias_been_considered by several au- 
thors (ICampo. Jr and" Oliveiral [2 003; llzumida et all. 



119971 : ISindel et all 120051 ). ISindel et all (|2005l ) calculated 
G'(u)), in the Kondo regime at T = 0, and extracted also 
the current noise 



C{u>) 



+ OC 



dte luJt [(I{0)I(t)} - (I)' 



(126) 



by making use of the fluctuation-dissipation theorem 



„hw/k B T _ ^ 



(127) 



The conductance and spin-resolved conductances of 
single-level quantum dots in a magnetic field have also 
been calculated and a strong spin-filterin g effect has been 
observed in the mixed valence regime (|Costil . l200ll) . For 
spin-filteri ng effects in quantum dots with fer r omag netic 
leads see lMartinek et all (|2003l ): ISimon et all (|2006l ). 

One of the hallmarks of the S = 1/2 single-channel 
Kondo effect is the flow of t he exchange coupling to 
strong coupling ( Wilsonl , ll975at) . This can be interpreted 
as resulting in a phase shift of the condu ction electrons a t 
the Fermi level, at T = 0, of S a = tt/2 (jNozieresl fl97i) . 
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A direct measurement of this phase shift is possible if one 
embeds a quantum dot in one arm of an Aharanov-Bohm 
interferometer. Assuming a single- level Anderson model 
for the quantum dot a nd a multi-terminal open geometry, 
iGerland et al. (2000) carried out NRG calculations for 
the interference term, Gab, whose measurement can be 
us ed to extract 5 a . A simil ar set-up has been investigated 
by iHofstetter et all ( 200ll ) for the flu x dependenc e G(4>) 
of the conductance at T = and by Kang et al\ (12005T) 
for the complex transmission. Ilzumida et all ( 1997t ) cal- 
culated G(</>) for two single-level quantum dots embedded 
in the arms of an Aharanov-Bohm interferometer. This 
model, reduces, in general, to a two-channel two-orbital 
Anderson model, which we discuss next. 



b. Two-level quantum dots A quantum dot with two ac- 
tive levels for transport introduces some new physics due 
to the competition between the level spacing 5 = £2 — £\, 
the charging energy Ec and the Hund's exchange Jh- In 
particular, a Kondo effect with an even number of elec- 
trons on the dot can be realized. This can occur when the 
dot is occupied with two electrons and S < 2Jh so that 
the ground state of the dot has S = 1. Such a two- level 
dot will, in general, couple to two channels so a S = 1 
Kondo effect will result, leading to a singlet ground state 
and an enhanced conductance G(T) at low temperatures. 
In the opposite case 5 > 2 Jh the dot will have S = 0, the 
Kondo effect is absent and the conductance will be low. 
This behavior is believed to have been measured in the 
experiments of Sa saki et all ( 2000f) on vertical quantum 
dots, where a magnetic field was used to decrease the en- 
ergy splitting, Ats = 8 — 2Jh, between the triplet and 
singlet states, thereby leading to the above mentioned 
crossover behavior in t he conductance at Af = 2. The- 
oretical calculations by Ilzumida et all (2001b), shown in 
Fig. [TH1 are consistent with the experimental results. 

The singlet-triplet crossover behavior in a two-level 
quantum dot can become a quantum phase transition for 
the special case where only one conduction channel cou- 
ples to the leads, e.g. when all lead couplings are equal 
(IHofstetter and Schoellerl . 120021 iPustilnik and Glazmanl . 
|200lD. In this case, for large Hund's exchange, an ef- 
fective single-channel 5 = 1 underscreened Kondo model 
results which has a doubly degenerate many-body ground 
state. For small Hund's exchange, a model with S = re- 
sults having a non-degenerate many-body ground state. 
A sharp transition separates these two different ground 
states. As discussed above, however, two channels will, 
in general, couple to the dot and this will result in perfect 
screening of the S — 1 so that the ground state is always 
a singlet. Nevertheless, proximity to the singlet-triplet 
transition can still be seen as signatures in various quan- 
tities, such as a non-monotonic dependence of the con- 
ductance as a functi on of magnetic field on the triplet 
side of the crossover ( Hofs tetter and Zarandll200l . Ex- 
periments on la teral quantum dots at th e singlet-triplet 
crossover point ( van der Wiel et ~al\ , l2002h show behavior 
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uiF.j/i: 

FIG. 16 The singlet-triplet crossover in the linear conduc- 
tance of the two-channel two-orbital Anderson model includ- 
ing a Hund's exchange Jh, intra- and inter-orbital Coulomb 
energy, U, levels Ed ± S/2 and temperature T. Top/bottom 
panels shows G(T) on the triplet /singlet side of the crossover 
and (fx - s d )/U = 0.5, 1.5, 2.5 correspond tojV= 1.0,2.0,3.0 
electro ns on the dot. Adapted from ISakai and Izumidal 
(|2003h . 



in the differential co nductance similar to pred ictions for 
the spectral density ( Hofstetter and Schoellerl . [2002h . 

The above is only a brief account of the simplest nanos- 
tructured devices studied using the NRG. Further appli- 
cations include numerous studies of double -dot systems , 
inclu ding realizing an SU(4) Kondo state ( Borda et all , 
|2003[ ) and qua ntum critic a l points of two-impurity 
Kondo models (iGarst et all , 12004 IZarand et all , [20061 
IZhu and Varmal 120061) (see Sec. HV.C.ll). st atic and 
dynamics of double-dots (jGalpin et all , 2006al lrJ). dou- 
ble dots with only one dot coupled to the leads 
(jCornaglia and Grempell . l2005rj ). applic ations to quan- 
tum t unneling in molecular magnets ( Romeike et all , 
l2006allbh . a novel Kondo effect ha a v = 1 integer 



quantum Hall system ( Choi et all . l2003bl ) and the con- 
ductance of ultrasmall tu nnel junctions ( Frotal . l2004t 
iFrota and Flensberd . fl992[) . 



B. Two-channel Kondo physics 



iNozieres and Blandinl (Il980l) have proposed a variation 
of the Kondo model in which the localized spin couples 
to two conduction bands. The Hamiltonian of this two- 
channel Kondo model is given by 



H 



^ £kC kaa 



(128) 



with a = 1,2 the channel index and S (s Q ) the spin 
operators of the impurity (the conduction band electrons 
at the impurity site with channel index a). 
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An important feature of this model is the overscreen- 
ing of the impurity spin: in the strong-coupling limit, the 
spins of both conduction bands try to screen the impu- 
rity spin, so that again a net spin 1/2 object is formed. 
In other words, the strong-coupling fixed point at J = oo 
(which gives rise to the Fermi-liquid fixed point in the 
single-channel case) is unstable and an intermediate- 
coupling fixed point is realized. This new fixed point 
shows a variety of non-Fermi liquid properties such as 

• a divergence of the specific heat ratio C/T = 7 cx 
In T and of the spin susceptibility x ^ l n T for T — > 
0; 

• an anomalous Wilson ratio R = \/l = 8/3, in 
contrast to the result for the standard Kondo model 
i? = 8/4 = 2; 

• a zero-point entropy of \ ln 2, indicating that 'half '- 
fermionic excitations (Majorana fermions) play a 
crucial role for the structure of the fixed point. 

We have discussed these features already in the section on 
the calculation of thermodynamic and static quantities, 
see Fig.[SJ An extensive review of the two-channel Kondo 
model, its physical properties and its relevance for non- 
Fe rmi liquid behavior in real materials has been given 
in ICox and Zawadowskil ( 1998| ). This paper also reviews 
the earlier NRG-calculations for this model. 

Historically, the two-channel Kondo model has been 
the first application of the NRG to a quantum impurity 
model in which the physics is not governed by the Fermi 
liquid fixed point of the standard Kondo m odel. In this 
sense, the early work of ICragg et al\ (|l98fj ) on the two- 
channel Kondo model opened the way for a variety of in- 
vestigations of more complex impurity models, displaying 
both Fermi liquid and non-Fermi liquid fixed points. Due 
to the importance of this and following work, we want to 
focus this section purely on two-channel Kondo physics 
and shift the discussion of other multi-band models to 
the section on orbital effects (Sec. IIV.DI) . 

As discussed already in Sec. HU the truncation of states 
within the iterative diagonalization scheme severely lim- 
its the applicability of the NRG to mu lti-band models. In 
the calculations of lCragg et all (1980), the iterations were 
observed to break down only after a few (approximately 
twelve) steps. The source of this problem is mainly the 
small number of states (N s w 400) used in this work, 
which would correspond to keeping N s sw V400 = 20 
states in a one-channel calculation. Specific symmetries 
of the two-channel Kondo model, such as the total ax- 
ial charge, have bee n used to reduce the matrix sizes in 
the diagonalization ( Pang and Coxl . IT99ll ). but later cal- 
culations showed that by simply increasing the number 
of states, the iterations can be stabilized sufficiently. In- 
dependent of the value of N s , it is important to avoid any 
symmetry breaking due to the truncation of states. 

In order to approach the no n-Fermi liquid fixed point 
within only a few i terations, ICragg et al\ ( 1980f) and 
Pang and Cox (|l99lh used large values of the exchange 



coupling J or large values of the discretization parame- 
ter (up to A = 9.0). Nevertheless, these calculations give 
the correct fixed point spectrum of the (isotropic) two- 
channel Kondo model with the characteristic structure of 
excitations at energies 1/8, 1/2, 5/8, 1, etc., at least for 
the lowest-lying excited states. Figure [T71 shows a typical 
flow diagram for parameters J — —0.051?, where 2D is 
the bandwidth of the featureless conduction band density 
of states, A = 4, and N s = 4900, for both even (dashed 
curves) and odd N RG iterations ( full cu rves) (for similar 
plot s, see Fig. 1 inlCragg et al\ (|l980l ) and Figs. 1 and 
2 in iPang and Coxl ( 199lf 0. After some initial even-odd 
oscillations, the flow reaches the non-Fermi liquid fixed 
point which does not show any even-odd effect. Note that 
this feature is by no means related to the non-Fermi liq- 
uid properties of the model; it just reflects the fact that 
in each iteration, two sites are added to the chain so that 
(for particle-hole symmetry) the number of electrons in 
the ground state is always even. 

Comparison with conformal field theory calculations 
(lAffleck et all Il992h gave an excellent agreement with 
the NRG for both the excitation spectrum and the cor- 
responding degeneracies. Such a comparison, however, 
requires the extrapolation of t he NRG fixed point spec- 
tra for A — > 1 (see Fig. 9 in lAffleck et al\ (|1992h ; the 
analysis is not quite satisfactory for A < 2 and it would 
be interesting to repeat these calculations using larger 
values of N s ). This work, and the previous paper by 
IPang and Coxl (|l99lh . also focussed on the stability of 
the non-Fermi liquid fixed point against various pertur- 
bations. As it turns out, the non-Fermi liquid fixed point 
is stable against anisotropy in the exchange interaction 
(J z 7^ JjJ but unstable against both the presence of a 
magnetic field and the lifting of the exchange symmetry 
between the two channels (J a ^ J&). In the latter case, 
a temperature scale T* cx (J a — J&) 2 for the crossover 
between the non-Fermi liquid fixed point at intermediate 
temperatures and the stable Fermi liquid fixed point at 
T — > has been found. These instabilities have been later 




FIG. 17 Flow diagram of the lowest lying many-particle levels 
for the isotropic two-channel Kondo model (even iterations 
dashed curves, odd iterations full curves). 
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invest igated in more detail in lYotsuhashi and MaebashU 
( 2002f ). via the calculation of the impurity entropy and 
the crossover temperature. 

Further investigations concerning the stability of 
the non-Fermi liquid fixed point have been per- 
formed in Pand (1 1994ft (flavor e xchange coupling) and 
iKusunose and Kuramotol ( 19991 ); IKusunose et all ( 19961 ) 
(effect of repulsion among conductions electrons and po- 
tential scattering). 

A 'pedestrian' approach for the understanding of 

the two-chann el Kondo model was introduced by 

IColeman et al. I (119951) . The authors of this work argued 
that the two conduction bands in the two-channel Kondo 
model can be replaced by a single conduction band, with 
a coupling between impurity spin S to both spin a and 
isospin r of the conduction band. The isospin r takes 
into account the charge degrees of freedom of the con- 
duction band and the compactificd 'cr-r' model takes the 
form 



H = J2 ekclrCk* + J (S • a + S • r) 



(129) 



ho 



This model can be related to an Anderson-type model 
(the 'Q(3)-symmetric' Anderson model) via a Schrie ffer- 
Wolff transformation ( Coleman and Schofieldl . IT995I ) . 

It has been later veri fied with the NRG approach 
( Bulla and Hewsonl . Il997l ) that these compactified mod- 
els indeed show many of the anomalous non-Fermi liq- 
uid properties of the two-channel Kondo model, although 
these models do not allow for an overscreening of the 
impurity spin. Furthermore, the structure of the non- 
Fermi liquid fixed point has been studied in detail. It 
turns out that the many-particle spectrum of this fixed 
point is composed o f single Majorana fermion excitations 
(iBulla et all . Il997a[) . This information can then be ex- 



tended to the fixed point structure of the two-channel 
Kondo model which can be described by two towers of ex- 
citations which are both composed o f Majorana fermions, 
see Sec. VIII in lBulla et all il997al) . 

Naturally, we expect that the non-Fermi liquid prop- 
erties of the two-channel Kondo model are also visible in 
its dynamic properties, but, unfortunately, detailed and 
comprehensive NRG calculations for the dynamics have 
not been performed so far. Apart from a brief sketch of 
the resu lts for the T-matrix and the magnetic suscepti- 
bility in lSakai et all ( 1993al) . the published data are only 
for models equivalent to the two-channel Kondo model 
in certain limits. 

It has been argued in iBradlev et all |l999) that the 
dynamical spin-susceptibility x( tJ ) °f t ne compactified 
models introduced above is exactly equivalent to that 
of the two-channel Kondo model, and that this equiva- 
lence holds for the full frequency range. The NRG-results 
show, for example, a m(u;)-diver gence of x'(u > ) for u> — > 0, 
in agreement with the results of lSakai et ofl (Il993al) . On 
the other hand, there is no co unterpart of t h e sing le- 
particle dynamics calculated bv lBradlev et all ( 1999f ) in 
the two-channel Kondo model. 



The two- ch annel Anderson model investigated in 
lAnders! (12005D is connected to the two-channel Kondo 
model via a Schrieffer- Wolff transformation (note that 
this only holds when the impurity degrees of freedom 
in the Anderson model are written in terms of Hub- 
bard operators which include the channel index). This 
connection is clearly visible in thermodynamic proper- 
ties, such as the zero-point entropy of ^ln2. Again, 
the single-particle dynamics (spectral function and self- 
energy) do not have a counterpart in the two-channel 
Kondo model. Concerning the results f o r the dynamic 
susceptibility x( w ) presented in lAndersl ( 20051 ). a com- 
parison to the corresponding results of the two-channel 
Kondo model has not yet been done. 

There is an ongoing discussion about the observability 
of two-channel Kondo physics in experiments for a vari- 
ety of systems. Let us stress here that the instability of 
the non- Fermi liquid fixed point itself does not exclude its 
observation. As for any system with a quantum critical 
point, the corresponding anomalous properties dominate 
a significant fraction of the finite-temperature phase di- 
agram (determined by the critical exponent) so that a 
precise tuning of the Hamiltonian parameters is not re- 
quired. Nevertheless, two-channel Kondo physics is now 
mainly discussed within systems in which alternative de- 
grees of freedom (such as orbital quantum numbers) take 
the role of spin or channel in the Hamiltonian Eq. (| 1281) ; 
one example is the quadrupolar Kondo m o del w hich is 
discussed in detail in I Cox and Zawadowskil ( 1998|) . 

Here we want to briefly discuss NRG calculations 
for two-chann el Kond o phys i cs in quantum d o t syste ms 
(|Anders et all , |2004 120051 : iLebanon et dl , l2003al lbl). 
Within a model of a quantum box coupled to the leads via 
a single -mode point contact (see Fig. 1 in lLebanon et all 
(|2003bl) ). the physics at the de generacy points of the 
Coulomb blockade staircase can be directly connected 
to that of the two-channel Kondo model. Here the two 
charge configurations in the box play the role of the im- 
purity spin and the physical spin of the conduction elec- 
trons corresponds to the channel index. For such a sys- 
tem, the NRG allows the non-perturbative calculation 
of the charge inside the box and the capacitance in the 
whole parameter regime. The results show, for example, 
that the shape of the charge steps is governed by the 
non-Fermi liquid fixed point of the two-channel Kondo 
model. 

To conclude this section let us mention that there are 
models involving a more complicated orbital structure 
of the impurity - including, for example, excited crys- 
talline electric field levels - which reduce to the two- 
channel Kondo model in certain limits or which dis- 
play non-Fermi liquid fixed points of the two-channel 
K ondo typ e. NRG-studies of such models can be found 



i nlHattoril(l2005l);lKoga and Coxl (ll999l);lKoga and Shibal 
(I1995L 119961 ); ISakai et all (|l997l) . Overscreening can also 



be realized in single-channel models when the conduction 
electron spin exeeds the i mpurity sp i n, for a discussion of 
this issue see for example iKim et oLl ()l997D . We note also 
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a recent study ( Kolf and Krohal . I2006T) showing an expo- 
nential dependence of the Kondo scale on —1/JNp and 
—JNp, for small and large coupling cases respectively, 
which may explain the absence of a broad distribution 
of Kondo scales in nanoconstrictions with two-channel 
Kondo impurities. 



C. Impurity quantum phase transitions 

In this section, we focus on models which, as a func- 
tion of one or more couplings in the model, give rise 
to a phase transition in the ground state. Typically, 
this is due to a competition between the Kondo effect 
on the one hand, which tends to favor a strong-coupling 
ground state with a screened or partially screened local 
moment, and some competing mechanism, which leads 
to a ground state with a free or almost free local mo- 
ment. In general, such phase transitions are termed 'im- 
purity quantum phase t r ansitions' (fo r recent reviews, see 
IBulla and Voital (|2003l) ; IVoital (|2006l) ). as they are only 
observable in the impurity contribution to physical prop- 
erties and not connected to possible phase transitions in 
the bulk to which the impurity couples. 

As impurity quantum phase transitions are usually 
associated with a vanishing low-energy scale, the NRG 
method is ideally suited to their investigation, allowing 
their detection and characterization with very high ac- 
curacy. This is most evident for continuous transitions 
where the critical exponents connected to the quantum 
critical point can only be calculated when a large range of 
energy or temperature scales is accessible. In this section 
we give an overview on NRG-results for multi-impurity 
models (Sec. HV.C.lj) . models with locally critical behav- 
ior (Sec. HV.C.2j) . and models with magnetic impurities 
in superconductors (Sec. HV.C.3|) . Note that impurity 
quantum phase transitions are also observed in models 
which arc discussed in other sections of this review: the 
non-Fermi liquid fixed point of the two-channel Kondo 
model Sec. IIV.BI can be viewed as a quantum critical 
point when the control parameter 'channel anisotropy' 
is tuned through zero; locally critical behavior is also 
connected to models with coupling to a bosonic bath as 
discussed in Sec. IIV.E1 



1. Multi-impurity physics 

An early extension of the NRG to more complex 
systems was the study of the two-impurity S = 
1/2 K ondo model ( Jones and Varmal . Il987l : IJones et~al\ . 
Il988l ). whose Hamiltonian is given by 



(130) 



Here, s(R;) is the conduction electron spin density at 
the impurity site R/ and Jk > is the antiferromagnetic 
Kondo exchange. The first two terms in Eq. (|130p are suf- 
ficient to generate an indirect RKKY interaction Irkky 
between the impurity spins. In some contexts a direct ex- 
change intera ction among the impurit y spins of strength 
Id can arise (|Jones and Varmal Il987f )1. so the last term 
has been added. The net effective exchange interaction 
between the spins is given by I c g = Id + /rkky and 
can be either ferromagnetic J e ff < or antiferromagnetic 
/off > 0. The properties of the model then depend solely 
on the ratio / e ff/7k, where Tk is the single-ion Kondo 
scale, and the details of the dispersion relation The 
model in E q. (11301) als o aris es in the Schrieffer- Wolff limit 
( Schrieffer and Wolfil . [1966) of the two- impurity Ander- 
son model, which in the notation introduced in Eq. @ 
reads 



H = 



1 = 1 a 



X £ kcL C ka + X X £ l°flho 
ka 1 = 1 

2 



i=i 



I — 1 ka 

+/ D S(Ri) ■ S(R 2 ) 



(131) 



The motivation to study such two-impurity models origi- 
nally arose in the context of heavy fermions. In these sys- 
tems, the competition between the local Kondo exchange 
and the intersite RKKY interaction is expected to lead 
to a phase transition between non-magnetic and mag- 
netically ordered ground states as a function of I e g/T K 
( DoniachT , ll977l h The nature of this quantum phase tran- 
sition remains an open que stion in heavy fermion physics 
v. Ldhnevsen et~al\ , 120061 ). It is therefore of some in- 



i=i 



-7 D S(R 1 )-S(R 2 ). 



terest to investigate the possibility of a transition in the 
two-impurity problem as this might shed light on the 
p hysics of heavy fer mions . 

IJones et ali ( 19881 ) have established that such a phase 
transition can occur under certain conditions, see be- 
low, in the particle-hole symmetric two-impurity Kondo 
model. This can be seen b y considering the s trong- 
coupling limits J e ff - * ±oo ( Affleck et all [l995). For 
I c g — > — oo the two spins combine to form a spin S = 1 
interacting antiferromagnetically with two conduction 
channels (characterized by even/odd parity) with, in gen- 
eral, energy dependent couplin g strengths J e (k),J (k) 
replaced in (jjones et all Il988l ) by constants (see be- 
low). The resulting two-stage Kondo effect progressively 
screens the S = 1 spin down to a singlet and leads to a 
Fermi-liquid ground state characterized by phase shifts 
<5 e , for electrons in the even/odd parity channels. The 
assumed particle-hole symmetry and the nature of the 
strong-coupling ground sta te ensures t h at th ese phase 
shifts will be exactly tt/2 (jMillis et all Il990f ). In the 
other limit, J c ff — > oo, the spins form an intersite sin- 
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glet S = and the Kondo effect is absent so that the 
phase shifts are exactly zero. Since the fixed points at 
leg = ±00 are both stable and characterized by different 
(constant) phase-shifts, it follows that there can be an 
unstable fixed point at some critical intermediate value 
of the intersite exchange, I c , at which the phase shifts 
change discontinuously. This phase transition has also 
been found in th e particle-hole symmetric two-impurity 
Ander son model (iPaula et all [T999: Sak ai and Shimizul . 



1992af ). IJones et al\ (|l988h estimates I C /T K « 2.2. 



The associated critical point has been characteri z ed us- 
ing conformal field theory ( Affleck and Ludwisi 1 1 99*21 ; 
lAffleck et all Il995h and bosonization ( Ganl . 19951 ) . and 
the physics is found to be similar to that of the two- 
channel Kondo model. In particular, the staggered sus- 
ceptibility, Xs(T), diverges logarithmically at low tem- 
perature and the residual entropy has the same va lue as 
in the two-channel Kondo model S(T = 0) = In \[2 (|Ganl . 
Il995l) . In contrast C{T)/T — 7 is predicted to remain fi- 
nite exactly at the critical point, in contrast to the behav- 
ior in the two-channel Kondo model. Close to the critical 
point, conformal field theory predict s 7 ~ (I P f? — I r )~ 2 



in agreement with numerical results ( Joned . Il990) . The 
level structure of the fixed point at quantum criticality 
agrees well with NRG calculations and i s rather com- 
plex, exhibiting a hidden SO (7) symmetry ( Affleck et all 
Il995l) . 

For generic situations, the natural energy dependence 
of J e j obtained from transforming the Kondo model 
(| 1 30[) or Anderson model (|131[) as described below in 
Eq. (|132[) . as well as a charge transfer term of the form 
t^2a(fiaf2<j~\~h j .c.) in the two-impurity Anderson model, 
breaks par ticle-hole symmetry and destroys the criti- 
cal point (|Affleck et all Il995t ISakai et all Il990fh A 
similar charge transfer term involving conduction elec- 
trons has the same effect in the two-impurity K ondo 
model (|Zarand et all l2006t Iziiu and Varmal . [2006). Po- 
tential scattering, if it does not induce charge transfer, 
breaks particle-hole symmetry bu t may not af f ect th e 
critical point , for a discussion see Afflec k^ al\ (|l995h : 
IZarand et all (j2006l ); IZhu and Varmal (|2006f) . Thus, in 
general the quantum phase transition discussed above 
will be absent in the two-impurity models (| 1 30|) and 
(|131[) . although signatures of it might still be ob- 
servable as crossover behavior in various properties. 
We note briefly here the case of Ising coupled spins 
IS' z (Ri)S' z (R2). In this case, the ground state for large / 
will be doubly degenerate as compared to that for small 
/ where the Kondo effect screens the individual spins to 
singlets, so a quantum phase transition separating these 
two different ground states arises and is found to be of 
the Kosteritz-Thouless type ( Garst et aZ.1 . 120041 ). 

In order to formulate Eq. (|130p or (|131[) as a linear 
chain problem for treatment with the NRG, one needs 
an orthonormal basis set. The local conduction electron 
states on the i mpurity sites in Eq . (1130 1 are not orthog- 
onal. Following I Jones and Varmal ( 1987 ). the Kondo ex- 
change part of Eq. (|130|) is rewritten in terms of orthonor- 



mal even (e) and odd (o) parity states for the conduc- 
tion electrons. This results in more complicated interac- 
tion terms; in particular, one will obtain two exchange 
couplings J e / {k), with J e {k) ^ J (k) in general, which 
will d epend on momentum or energy (jjones and Varmal . 
I1987D . The precise form of J e / (fc) will depend on the 
details of the band-structure of the conduction electrons. 
For free electrons i n D = 3 it can be ap proximated as 
(|Jones and Varmal . Il987t ISakai et all Il990f) 



Je/o{k) ~ Jk 



1± 



ikR 



kR 



where R — |Ri — R 2 |. IJones and Varmal ( 19871 ) used 
constant couplings J e / (k) Je/o{kp) to obtain for the 
interaction part of the Hamiltonian (|130|) 

-Hint = S (e) • [ject a a a pc ef3 + J cl a a af 3C of . 



a/3 

S (o) • Y [W J eJoct a a a pc o/3 + h.c 

a/3 

J D S(Ri) ■ S(R) 2 , 



(132) 



where S( e / ) := S(R X ) ± S(R 2 ). The conduction elec- 
tron Hamiltonian now consists of two decoupled linear 
chains with even and odd parity symmetry. By neglect- 
ing the energy dependence of the couplings an particle- 
hole symmetric m odel results. This is the form used by 
I Jones et~ai\ ( 19881 ) to investigate the phase transition dis- 
cussed above. The results of retaining the full energy 
dependence of the couplings, using for example the for- 
mulation of Sec. [Til will be described below. We note 
here that from the NRG point of view the two-impurity 
models (|131l) and (|130[) present a challenging task be- 
cause, as in the case of the two-channel Kondo model 
(|128p . the "impurity" now couples to two semi-infinite 
chains. Consequently, the Hilbert space grows by a fac- 
tor 16 in each NRG step. While this is still manage- 
able with modern computer resources, it is apparent that 
larger clusters or more complex situations quickly become 
too expensive to be treated with NRG with sufficient ac- 
curacy, although the flow of the many-body eigenstates 
can still be used to identify fixed points and thus qual- 
itatively describe the physics of more complicated sys- 
tems, like the two-cha nnel two-impurity Kondo model 
( Ingersent et al. , 19921) and the three- impurity Kondo 
model (jlngersent et all 120051 ). In high-symmetry situ- 
ations even a calculation of thermodynamical quantities 
has been performed rece ntly for three-impurity models 
(jZitko and Boncal l2006al) . However, for a reliable calcu- 
lation of dynamics or in situations with less symmetries 
in the system - e.g. in an external magnetic field - ad- 
ditional tools like the ones described in section IIII.A.4I 
allowing one to work with large A> 1 and so maintain 
low truncation errors should be useful. 

The generic two- impurity Anderson model (|131[) . in- 
cluding a charg e transfer term, has been st u died by Saka i 
and coworkers ( Sakai and Shimizul . IT992a.b; Sakai et all . 
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1990) using the NRG. Single particle and magnetic exci- 
tation spectra were calculated and in the case of particle- 
hole symmetry. ISakai et all ( 1993ah showed that on pass- 
ing through the transition, a peak in the impurity single- 
particle spectra sharpened at I c g = I c into a cusp 
and turned into a dip for I c s > I c . In the generic 
case, the regime with Kondo screening, |J e g| <C Tk, 
and the non-local singlet regi me, J e g 



necte d via a smooth crossover dCampo. Jr. and Oliveira , 



2004: 


Saka 


i and Shimizul 


Silva et al.. 


1996h. 



Results from thermodynamic calculations are shown 
in Fig. 1181 for the squared effective magnetic moment 
taken from lSilva et all ( 1996) for the two- impurity Kondo 
model with Id = 0. In these calculations I c s = /rkky 
and the energy dependence of J e / Q (k) is crucial to gen- 
erate the intrinsic RKKY exchange interaction Irkky- 
Using the result f or free electrons in three dimensions 
( Sakai et all , fl9 90) . an approximate formul a for the en- 
ergy dependence of the coupling constants is ( Silva et~all . 
1996) 



J e /o(e) = J K 1 ± 



sin [k F R(l + e) 
k F R(l +e) 



(133) 



with e G [—1, 1], kp the Fermi momentum of the conduc- 
tion states. For the derivation of Eq. (|133[) a linearized 
dispersion relation et ~ -j^(k — k F ) was assumed and 
D = 1 used as energy scale. 

Dependent on the value of k F R, different regimes can 
then be identified (see, for example, Fig. [T5|) : For k F R — > 




FIG. 18 Effective local moment fi 2 (T) := T ■ Xim P (T) for 
the two-impurity Kondo model (|131fl for the three r egime s 
described in the text. Figure taken from ISilva et all (1 19961 ). 
The arrows indicate the Kondo scale Tk = 1.4 x 10~ 4 , and 
the RKKY interactions for ferromagnetic (Jfm = —8 x 10 -3 )) 
and atiferromagnetic cases (Iafm = 3 x 10 -3 ). 



oo we have Irkky = 0, single-impurity physics domi- 
nates and no non-local magnet i c exch ange is generated, 
as expected ( Jones and Varmal Il987l) . For k F R = tt/2, 



the RKKY exchange Irkky = ^fm is ferromagnetic, 
with | Irkky | /Tk ^ 1 and a two-stage screening scenario 
arises: First, the system locks into an S — 1 state at high 
temperatures due to the RKKY interaction. In the in- 
termediate temperature regime this triplet is screened to 
a doublet via the even channel, which then is further 
screened to a singlet by the odd channel. Finally, for 
k F R = n, the RKKY exchange Irkky = Iafm is anti- 
ferromagnetic with Irkky /Ik 3> 1 and a non-local sin- 
glet is formed eventually. Similar results for entropy and 
specific heat of the two-impurity Kondo model, exhibit- 
ing a smooth change of physical properties with changin g 
Irkky, can be found in lCampo. Jr. and Oliveiral (2004). 

We note, here, that while two-impurity models with 
energy independent coupling constants are crude ap- 
proximations in the context of bulk Kondo impuri- 
ties and heavy fermions, these can, however, be re- 
alized in quantum dots. Correspondingly, they have 
been proposed to describe various extensions of sin- 
gle quantum dots and studied in this context with 
NRG by several groups jBoese et a/.l . l2002l Borda et al. 



2003; Hofstcttcr and Schoeller, 2002; Voita et al, 2002a: 



Zarand et all . 120061 : IZitko and Boncal l2006bt) over the 
past years. Since modern nanostructure technol- 
ogy permits a rather broad tailoring of such meso- 
scopic objects, the models discussed typically in- 
troduce additional interactions as compared to the 
conventional two-imp urity Anders o n mo del (II 3 1[) like 



capacitive couplings dBoese et al. 



2002; IBorda et all . 



12003 ; Hofstetter and Schoellerl. 2002) or direct hoppin g 
(|Dias da Silva et all . 120061 : IZitko and Boncal l2006bth 
Consequently, these extended models show a much 
larger variety in intermediate- and low-temperature 
fixed points than the bare model (|13ip . ranging 
from conventio nal Kondo effect over a two-stage 
Kondo effect ( Javaprakash et all . Il98lt IVoita etall . 
l2002al ), two-channel physics a s intermediate fixed- 



point ( Zitko and Bonca . |2006a 



transitions (IVoita et al. 



B) to quantum - phase 



2002a; Zarand et 



IZhu and Varmal . l2006t IZitko and Boncal l2006bfK 



m-pnase 
, 2006; 



2. Local criticality 

The term 'local criticality' has been first used in the 
context of phase transitions in certain h e avy f e rmion 
systems, such as CeCufi -gAua, ( Si et all . l200ll Il999t 
Iv. Lohnevsen et all . |2006[ ) . It has been argued that the 
quantum critical point separating the magnetically or- 
dered and the paramagnetic phases at T = is charac- 
terized by critical excitations which are local. This ob- 
servation raised a considerable interest in models which 
show such locally critical behavior: these are either lat- 
tice models studied within certain extensions of DMFT 
(see also Sec. IV.Bf or impurity models as discussed in this 
section. Such impurity models might not be directly con- 
nected to the locally critical behavior in heavy fermion 
systems, nevertheless, the insights gained in studying im- 
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purity models might be helpful in constructing theories 
for lattice systems (for a general discussion of the rela- 
tion between quantu m impurity ph ysics and the physics 
of lattice systems, see IBullal ( 20061 )). 

Let us focus here on the soft-gap Anderson m odel, 
originally proposed bv lWithoff and Fradkinl (1990). The 
Hamiltonian is the same as the one for the standard 
single-impurity Anderson model Eq. (JJ), but the hy- 
bridization function is assumed to have a power-law form 



fixed points is particularly complex in the asymmetric 
case, see, for example, the schematic flow diagrams of 



A(w) = A M r with r > -1 



(134) 



either valid over the whole frequency range or restricted 
to some low-frequency region. The competing mech- 
anisms leading to a quantum phase transition in this 
model are local moment formation (favored by increas- 
ing U) and screening of the local moments. For values of 
the exponent r > 0, corresponding to a soft-gap in A(w), 
there are less degrees of freedom available to screen the 
moment and a quantum phase transition occurs at some 
finite value of Ao. 

This quantum phase transition and the physical prop- 
erties in the whole parameter regime have been studied 
in detail with a var i ety of techniques (f or an overview, 



see iBulla and Voital (120031) ; IVoital 1)20061 ) and the intro- 
ductory parts in ILee et all (|2005h ). The NRG method 
has been particularly helpful to clarify the physics of the 
soft-gap Anderson model (and the related Kondo version 
of the model) as we shall briefly discuss in the following. 
The technical details necessary to apply the NRG to the 
soft-gap Anderson model have alread y been introduced 
in SecH see also lBulla eTaA (|l997bl ). 

Thermodynamic and static properties of the var- 
ious phases of the soft-gap Anderson and Kondo 
models h ave been p resent ed in Chen and Javaprakashl 
(ll995al). Ilngersentl (fl996l). IBulla et all (Il997bl). and 
iGonzalez-Buxton and Ingersent ( 19981 ). The most 
comprehensive review of these res ults is given in 
IGonzalez-Buxton and IngersentI ( 19981 ). This paper cov- 
ers the discussion of thermodynamic properties and the 
analysis of the various fixed points also for the under- 
screened spin-1 Kondo model and the (overscreened) two- 
channel Kondo model (both with a soft-gap in the con- 
duction electron density of states) . 

The key role of partic le-hole symmetry has been iden- 



ti fied in IngersentI (Il996f) and investigate d in more detail 



in lGonzalez-Buxton and Ingersent] (1998). As shown, for 
example, in Fig. 5 in this work, the line of quantum 
critical points separating the local moment (LM) and 
strong-coupling (SC) phases is restricted to < r < 1/2 
in the particle-hole symmetric case (for r > 1/2, only 
the LM phase exists). This is different in the asym- 
metric case where transition line extends up to r — ► oo. 
Particlc-hole symmetry also influences the physical prop- 
erties of the various fixed points. The symmetric SC 
fixed point, for example, shows a residual magnetic mo- 
ment of Ximp = r/(8k-BT) and a residual entropy of 
S'imp = 2rfcBln2, whereas both values are zero in the 
asymmetric SC fixed point. The appearance of unstable 



Fig. 16 in lGonzalez-Buxton and IngersentI (1998 

The impurity spectral function of the sym metric soft- 
gap An derson model was first investigated in lBulla et all 
( 1997bl ): the spectral function shows a divergence A(tu) oc 
\uj\~ t for both the SC and quantum critical phases 
whereas it goes as A(li) oc \uj\ r in the LM phase (for 
the behavior in the asymmetric case, see the discussion 
in Sec. HV.C.3[) . 

In the symmetric SC phase, the product F(uj) — 
c\oj\ t A(uj) (where the prefactor cancels the divergence in 
the spectral function) contains a generalized Kondo res- 
onance at the Fermi level with a pinning of F(lu = 0) 
(for a properly chosen constant c) and a width that goes 
to zero upon approaching the quantum critical point. 
This feature, together with the scaling properties and the 
low-energy asymp totics has been discussed in detail in 
IBulla et all (|2000h . based on both the results from NRG 
and fr om the local momen t appro ach (also described ear- 
lier in lLogan and Glossod ( 20001 )'). 

Dynamical properties at the quantum critical point 
are particularly interesting: llngersent and Si (2002) have 
shown that the dynamical susceptibility at the criti- 
cal point exhibits w/T-scaling with a fractional expo- 
nent, similar to the locally critical behavior in the heavy 
fermion systems mentioned above. This result also im- 
plies that the critical fixed point is interacting, in contrast 
to the stable fixed points (SC and LM) which both can be 
composed of non-interacting single-particle excitations. 

The interacting fixed point of the s ymmetric s oft-ga p 
model has been further analyzed in ILee et all (2005). 



The general idea of this work can be best explained with 
Fig. HU which shows the dependence of the many-particle 
spectra for the various fixed points on the exponent r. 




FIG. 19 Dependence of the many-particle spectra for the 
three fixed points of the particle-hole symmetric soft-gap An- 
derson model on the exponent r: SC (dot-dashed lines), LM 
(dashed lines), and the (symmetri c) quantum crit ical point 
(solid lines). (Figure adapted from lLee et all (|2005l )). 

For the limits r — > and r — > 1/2, the many-particle 
spectra of the quantum critical point approach those of 
the LM and SC fixed points, respectively. The devia- 
tions and splittings of the spectra at the quantum critical 
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point close to these limits can then be understood from 
a proper perturbational analysis using suitable marginal 
operators. Information on these operators can be ex- 
tracted from ep silon-expansion techniques, as shown in 
iLee et all (|2005l) . 

The case of negative exponents in the hybridization 
function, A(w) oc \ui\ r with — 1 < r < 0, where the soft- 
gap turns i nto a divergence at the F ermi level, has been 
analyzed in lVoita and B ulla (2002a|) in the context of the 
Kondo model with both ferromagnetic and antiferromag- 
netic values of J. The behavior of this class of models 
turns out to be rather complex, see the schematic flow 
diagrams of Fig. 1 in this work. A remarkable feature 
here is the appearance of a stable intermediate coupling 
fixed point with universal properties corresponding to a 
fractional ground-state spin. 

The case of a hard gap in the hybridization function, 
that is A(u>) — within a certain gap region around the 
Fermi level, can be viewed as the r — > oo-limit of the 
soft-gap case, provided the powerlaw is restricted to the 
gap region < E g /2, with E g the width of the gap. 

From a technical point of view, two different strate- 
gies have been developed to apply the NRG to the hard- 
;ap ca se. iTakegahara et al\ ( 19921 ) and lTakegahara et al\ 
1993h considered the case of a small but finite value of 
A(ui) in the gap region, A(w) = A for \uj\ < E g /2, and 
based their conclusions on the extrapolation A — > 0. In 
this approach, the standard NRG for non-constant hy- 
bridization functions as described in Sec. [TT] can be ap- 
plied. 

If, on the other hand, the value of A is set to zero 
from the outset , the NRG approach ha s to b e modified. 
As discussed in lChen and Javaprakash] (|1998), the loga- 
rithmic discretization of a A(uj) with a hard gap results 
in a discretized model which maps onto a chain with a 
finite number of sites, M, with E g of the order of A~ M . 
The iterative diagonalization then has to be terminated 
at site M. Thermodynamic properties at temperatures 
T < E g can nevertheless be com puted using the Hamil- 
tonian of the final iteration (see IChen and Javaprakash! 
(119981) where also a variety of correlations functions have 
been calculated for both the Kondo and the Anderson 
model with a hard gap). 

Certain features of the soft-gap case with finite r 
are also visible in the fully gapped case. As expected 
from the discussion above, there is no transition in the 
particle-hole symmetric case, but a transition exists as 
soon as one is moving away f rom particle-hol e sym- 
metry (IChen and Javaprakashl Il998t llngersenti Il996t 



ITakegahara et a/F 1992 ). This transition turns out to 
be of first order. 



3. Kondo effect in superconductors 

Let us now consider magnetic impurities in supercon- 
ducting hosts. In this case, the screening of the mag- 
netic moments competes with Cooper pair formation of 



the conduction electrons. We therefore expect a quantum 
phase transition from a screened phase to a local moment 
phase upon increasing the value of the superconducting 
gap, A, similar to the phase transitions in the soft-gap 
(and hard-gap) impurity models discussed in Sec. IIV.CT2l 
In fact, a relation between impurity models in supercon- 
ductors and those in metallic hosts with a soft or hard 
gap can be established as discussed below. 

The first applications of the NRG to magnetic im- 
purities in superconductors focused on the s-wave case 
(jSakai et all Il993bt ISatori et all Il992h . with the stan- 
dard Kondo Hamiltonian Eq. (|115p supplemented by the 
BCS pairing interaction 
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Several strategies have been developed to transform the 
Hamiltonian including the BCS-term onto a semi-infinite 
chain which can then be diagonalized iteratively in the 
usual way. ISatori et al\ |l992) performed a sequence 
of transformations, including a Bogoliubov and a par- 
ticular particle-hole transformation, to map the original 
model onto a Hamiltonian which conserves particle num- 
ber (this is somewhat easier for the numerical implemen- 
ta tion though not ab solutely necessary). The approach 
in lSakai et al. (1993b) leads to the same Hamiltonian, the 
difference here is that the Bogoliubov transformation is 
performed before the logarithmic discretization. In both 
cases, the semi-infinite chain contains a staggered poten- 
tial of the form 
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This term does not fall off exponentially as the other 
terms in the chain-Hamiltonian so that the NRG- 
iterations should be terminated a few steps after the 
characteristic scale ujn of the chain Hamiltonian iTv has 
reached the superconducting gap A. This procedure still 
allows to access the properties of the localized excited 
state within the energy gap whose position and weight 
can now be determined in the full parameter space (in 
co ntrast to pr e vious investigations, see the references 
m ISatori et al\ (|1992D 1. Figure [201 shows position and 
weight of the localized excited state as a function of Tk/A 
(Tk is determined from the corresponding Kondo model 
with A = 0). The position changes its sign when Tk 
is of the order of A (the precise value depends on the 
model), corresponding to a change of the ground state 
from doublet for small Tk/A to singlet for large Tk/A. 
This quantum phase transition can be characteri z ed as a 
level crossing transition (see Fig. 5 in lSatori et al\ (|l992h ) 
and is not connected to quantum critical behavior. 

These studies of impurities in s-wave superconductors 
hav e been later extended to mor e complex impurity mod- 
els. lYoshioka and Ohash i ( 1998) investigated the case of 
an anisotropic interaction between impurity and conduc- 
tion electron spin, with basically the same NRG approach 
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FIG. 20 Position E\> and weight (intensity) of the localized 
excited state as a function of 7k/ A for the Kondo model in an 
s-wave superconductor. At 7k /A ~ 0.2, the position changes 
its sign and th e weight jumps by a factor of 2, see also Fig. 2A 
in ISakai et all (|l993bl ). 



as in ISakai et all ( 1993bJ h The phase diagram of this 
model turns out to be much more complex than the one 
for the isotropic case. For example, two localized excited 
states with different energies appear in certain regions of 
the parameter space. 

lYoshioka and Ohashil ( 2000f) considered the Anderson 
version of the impurity model with coupling to an s- 
wave superconductor. From a technical point of view, 
this case is different to the corresponding Kondo model 
sinc e the sequence of tr ansformations used in, for exam- 
ple, ISakai et all (|l993bl ) now produce an extra impurity 
term of the form <5(/j/[ + h.c), so that the whole Hamil- 
tonian does no longer conserve charge (note that the pa- 
rameter 6 is zero for the particle- hole symmetric case). 
The results for the Kondo regime of this model are, as 
expected, th e same as those obtained pre viously, but the 
approach of lYoshioka and Ohashil (2000) also allows to 
study other parameter ranges of the model, such as the 
mixed valent regime. 

It is i mportant to n ote h ere th at the final Hamilto nian 
used in ISatori et all (|l992j ) and ISakai et all (|l993bl ) for 
the NRG iteration is the same as the one for an impurity 
in a non-superconducting host with a gapped density of 
states (which corresponds to the quasiparticle density of 
states of the superconductor) . In addition, the sequence 
of transformations also generates a potential scattering 
term. In the light of the results for the hard-gap impu- 
rity models (Sec. IIV.C.2[) this potential scattering term 
is essential to observe the quantum phase transition from 
a screened to an unscreened phase. 

The question now arises, whether the quasiparticle 
density of states can be used as the sole bath charac- 
teristic (possibly supplemented by a potential scattering 
term) in more general situations, such as impurities in 
unconventional superconductors. Before we address this 
issue, let us have a look at what happens when a simi- 



lar sequence of transformations as in the s-wave case is 
applied to impurity models in p- or d-wave superconduc- 
tors. 

iMatsumoto and Kogal (|200lL [2002) considered the 
Kondo model with a coupling of the impurity spin to 
superconductors with p x + ip y and d x 2_ y 2 + ial xy symme- 
try (with extensions to spin-polarized supercondu cting 
states investigated in lKoga and Matsumotol (l2002al) an d 
to S — 1 impurities in Koga and Matsumotol (2002b)). 
The quasiparticle density of states in these cases also 
shows a full gap, as for s-wave superconductors, but the 
sequence of transformations now results in a model with 
a coupling of the impurity to two angular momenta of 
the conduction electrons. NRG calculations for this two- 
channel model give a ground state which is always a spin 
doublet for arbitrary values of Xk/A, in contrast to the 
s-wave case, and no level crossing is observed. This is 
supported by calculations of the impurity susceptibility 
which show that the effective magnetic moment is always 
finite, although strongly redu ced with increasing Tk/A 
(|Matsumoto and KogaL bo02h . The authors of this work 
argued that the orbital dynamics of the Cooper pairs is 
responsible for the ground state spin. 

This interp re tation has been questioned in 
iFritz and Voital (|2005l) . where is was shown that, 
indeed, the local quasiparticle density of states of the 
superconductor is the only necessary ingredient in 
a number of cases, in particular for unconventional 
superconductors. A p plied to the model studied in 
IMatsumoto and Kogal ( 200lh . this means that the 
results of the NRG calculations for the effective two- 
channel model can also be understood from a single-band 
calculation where screening is absent for a hard-gap 
density of states and particle-hole symmetry. 

The results o f IFritz and Voital <|200.^ also have impor- 
tant consequences for the study of impurities in uncon- 
ventional superconductors with d x i_ y i symmetry. In this 
case, the mappings which have been used for the models 
discussed above result in an impurity model with cou- 
pling to infinitely many bands to which the NRG clearly 
cannot be applied. For certain geometries, however, it 
is sufficient to consider only the quasiparticle density of 
states which, for a point-like impurity, shows a soft-gap 
with exponent r = 1. 

This simplific a tion ha s already been used earlier in 
IVoita and Bullal ( 2002bl ) (at that time it has been ar- 
gued to be a reasonable approximation). The results 
of this work are therefore both valid for the soft-gap 
Kond o model and for impuriti es in d-wave superconduc- 
tors. IVoita and Bulk] (|2002bl) have motivated their in- 
vestigations with experimental results for non-magnetic 
impurities in cuprate superconductors which have been 
seen to generate magnetic moments. As discussed in 
this work, an effective model for this problem then takes 
the form of a Kondo model in a d-wave superconductor. 
Connections to experimental results can indeed be made 
within this framework. For example, the T-matrix T(w) 
displays a very narrow peak at finite frequencies with the 
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peak energy corresponding to the energy scale which van- 
ishes at the quantum phase transition from a screened to 
an unscreened moment. A very similar peak has been 
observed in STM-experiments. 

This work has been later extended in IVoita et all 
(l2002bl) . where the effects of local and global magnetic 
fields have been investigated. For the case of a local field 
h\ oc , the quantum phase transition for zero field persists 
for hioc 7^ 0, but for a global field, the quantum phase 
transition turns into a sharp crossover since the global 
field induces a finite spectral weight at the Fermi level. 

The investigations described so far are mainly applica- 
ble to impurities in the bulk or on the surface of a super- 
conducting host. A different geometry is realized in quan- 
tum dot systems (see Sec. IIV.A~3)l . For superconducting 
leads, such a setup introduces a new control parameter to 
the problem, that is the phase difference, $ = $l — $r, 
between the phases of the two superconducting leads. 
The resulting Josephson current, in particular the tran- 
sition from 0- to 7r-iunction b ehav ior, has been studi ed 
in detail in[Choi et all (|2004a| ) and lOguri eTall (|2004h . 

IChoi et al. (2004a) investigated various static and dy- 
namic properties for this geometry with identical s-wave 
superconductors as the two leads. For zero phase differ- 
ence, $ = 0, the local pairing correlation shows a sign 
change at Xk/A « 0.42. Physically, this is connected to 
the same quantum phase transition as described above 
since for $ = and identical leads the model can be 
mapp ed onto the same model as discussed in lSatori et all 
(1992). For finite phase difference (or for non-identical 
leads) the system remains a two-channel problem and 
the NRG analysis is more complicated. Nevertheless, de- 
tailed information on ground state properties such as the 
single-part i cle exci tation spectrum have been obtained in 
IChoi et all ( 2004al ) and interpreted as a phase-dependent 
form ation of And re ev bo und states. 

In lOguri et all ((20041) ■ the Hamiltonian of an An- 
derson impurity coupling to two superconducting leads 
has been considerably simplified by studying the limit 
|Al| 3> |ArJ in which the model can be mapped exactly 
onto a single-channel one with an extra superconducting 
gap on the impurity. Results for this limit show that the 
phase difference changes both the energy and the wave 
function of the bound state. In particular, the phase dif- 
ference appears to work against the screening of the local 
moment. 



D. Orbital effects 

1. Multi-orbital Anderson model 

The physics of the Kondo effect requires the existence 
of local magnetic moments, as realized, for example, in 
systems with open d or / shells, such as transition metal 
or rare-earth impurities in non-magnetic host metals. For 
such systems, the local Coulomb correlations and Hund's 
exchange determine the electronic structure of the impu- 



rity and they usually give rise to finite spin and orbital 
magnetic moments. Thus, a realistic description of such 
impurities in solids, requires taking both spin and or- 
bital magnetic moments into account. The same is true 
for the compounds of transition metal, rare-earth and ac- 
tinide elements, where the interplay of orbital and spin 
degrees of fre e dom g ives rise to very rich phase diagrams 
(jlmada et all . II998I ). Among the methods to theoreti- 
cally study the properties of these materials, the dynam- 
ical mean-field theory (see Sec. [V]) has become a stan- 
dard approach. Since in this approach one ends up with 
an effective quantum impurity problem which retains the 
full local orbital and spin structure of the original lattice 
system, the development of a reliable method to solve 
quantum impurity models with orbital and spin degrees 
of freedom is of crucial importance. 

In this section we, therefore, discuss the application of 
the NRG to situations where orbital and spin degrees of 
freedom are both present. Some orbital effects in quan- 
tum dots have been discussed in Sec. IIV.A.3I and will 
be discussed further below. The appropriate model is 
again a suitable extension of the single-impurity Ander- 
son model @ and is given by 
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where m labels the orbital degrees of freedom, and n 
fLafma, wi * n = \ E /L^/mfl- In addition to the 

intra-orbital Coulomb term U also occurring in @, the 
following interaction terms are present now: An inter- 
orbital Coulomb repulsion U' and an exchange term J. 
The exchange t erm we have s plit in accordance with stan- 
dard notation ( Imada et aZ.I . fl998h into a Heisenberg-like 



spin-exchange term (Hund's coupling) and an orbital ex- 
change term. To account for the proper combination of 
operators in the general exchange contribution, an addi- 
tional part proportional to J appears in the inter-orbital 
Coulomb term. For free atoms, rotational invariance usu- 
ally imposes U' = U — 2 J as constraint for the different 
Coulomb parameters ( Imada et a/.l .[l998). Further mod- 
ifications to the model (|f 35|) can be made, to take into 
account, for example, spin-orbit and crystal-field effects. 
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The mapping of (|135|) onto a linear-chain model, see 
Eq. of Sec. QTJ clearly leads to m semi-infinite con- 
duction chains coupled to the local Hamiltonian, which 
in turn means that at each step of the iterative diagonal- 
ization, the Hilbert space increases by a factor 4 m . For 
large m, this exponential increase in the number of states 
makes the NRG truncation scheme useless, because the 
number of states one can keep is way too small to al- 
low for a reasonable accuracy. Thus, calculations for the 
model (|135l) involving a full d shell (m = 5), or even a full 
/ shell (m = 7), and taking into account all interactions 
seem to be impossible. 

In practice, however, one is typically not interested in 
rotationally invariant situations, as described by (|135|) . 
but in situations where the impurity is embedded in the 
crystalline environment of a solid. The reduced point- 
group symmetry due to the crystalline electric field, then 
leads to a splitting of the orbital degeneracy. The en- 
ergy associated with this crystal-field splitting can be 
much larger than the temperatures one is interested in 
experiments, for example in 3d transition metals. Fur- 
thermore, the local Coulomb interaction can lead to a 
localization of electrons in the lower crystal-field multi- 
plets, as happens, for e xample, in the case of manganese 
in a cubic environment (I mada et aE[l998). In this case, 
these states form a localized spin according to Hund's 
rules. For manganese, for example, this results in a high- 
spin state (S = 3/2) of the threefold degenerate t 2g or- 
bitals, which couples ferromagnetically to the twofold de- 
generate e g electrons. Thus, the actual number of rele- 
vant orbitals, and thus the number of semi-infinite chains 
coupling to the local Hamiltonian, may be considerably 
reduced. Similar effects can be observed in the higher 
rare-earth elements, for example in gadolinium. 

In case the local point-group symmetry is reduced suf- 
ficiently, one may, in fact, be left with a localized spin S 
coupled to a single spin-degenerate, correlated orbital hy- 
bridizing with conduction states. This Kondo- Anderson 
model is given by © supplemented with the ferromag- 
netic exchange term —JhS ■ Sd, with Sd the spin-density 
of the correlated level in the Anderson m odel. Such 
a Kon do- Anderson model was studie d by ISakai et~al\ 
(|1996() and lPeters and Pruschkel (|2006f) and shown to ex- 
hibit different types of screening, ranging from conven- 
tional Kondo screening to two-stage screening and local 
singlet formation or two-channel Kondo effect. 



2. NRG calculations - an overview 

A first serious attempt to study ef fects of true orbita l 
degeneracy with NRG can be found in ISakai et all (|l989h . 
However, the authors did not study the full Hamiltonian 
(|I35[) . but a SU(N) version of it, using values of N rang- 
ing from two (i.e. the standard single-impurity Anderson 
model @) to five, representative for rare-earth ions like 
samarium or th ulium in solids under the influence of a 
crystalline field ( Allub and Aligial [l995; Shi mizu et all . 



I1990D . Since the SU(N) model has a large degeneracy 
of the individual levels, it allows for a considerable re- 
duction of the sizes of the individual Hilbert spaces in 
the diagonalization. This enabled the authors to use 
the NRG to calculate physical properties including dy- 
namical quantities and to study, for example, the de- 
velopment of the Kondo temperature or the behavior of 
the Abrikosov-Suh l resonance as f unctio n of degeneracy 
(Saka i et al .], 119891 ; IShimizu et al\ , [l990). Similar inves- 
tigations for the model (|135[) with fixed orbital degen- 
eracy to = 2 in the pre sence of a magnetic fi eld were 
presented in the work bv lZhuravlev et all (|2004 ). a com- 
pari son with STM experiments f or Cr(OOl) surface states 
by (|Kolesnvchenko et all , l2005h . and a detailed study 
of the dependence of the low-energy properties of the 
mu lti-orbital Anderson mode l for m = 2 with J > 
bv iPruschke and Bullal (|2005f) . Furthermore, from the 
fixed-point level structure, interesting information about 
quantities like "residual interac tions" in the heavy F ermi 
liquid state can be extracted (jHattori et al. I. l2005h . In 
particular for an f 2 ground state - as possibly realized 
in uranium compounds - a subtle enhancement of in- 



teror bital interactions can be observed (|Hattori et ali , 
I2005D . which can lead to superconducting correlations in 
a triplet channel when used as effective interaction in a 
model for heavy-fermion superconductivity. 

While these studies deal with the conventional 
Kondo effect in mul t i-orbi tal models, it was noted 
( Cox and Zawadowskil Il998l ) that for higher rare-earth 
and actinide elements the orbital structure in connec- 
tion with spin-orbit coupling and crystal-field effects can 
result in an orbital multiplet structure that leads to 
the two-channel Kondo effect (see Sec. IIV.BI). Multi- 
orbital models of that type w e re studied bv ISakai et all 
(119961 ) [Shimizu et all (Il998t). IShimizu et all (ll999aUbl). 
iKogal (|2000f) and lHattori and Mivakd (|2005f ) , covering a 
wide range of aspects possibly realized in actinide heavy- 
fermion systems. The authors could identify parame- 
ter regimes where non-Fermi liquid properties related to 
the two-channel Kondo effect can be observed and could, 
in addition, identify relevant symmetry breakings like 
crystal-field splittings or external fields that eventually 
lead to conventional Kondo physics below a tempera- 
ture scale connected to the energy scale of the symmetry 
breaking. 

As mentioned in the introductory remarks of this sec- 
tion, the crystal-field splitting is usually much larger than 
the relevant low-energy scales. However, this does not 
need to be true in genera l. Besides uranium-based com- 
pounds ( Kusunosa |2005| ) . a possible example where the 
crystal-field splitting can actually be of the orde r of th e 
Kondo scale is Cei_ 2; La a; Ni9Ge4 ( Scheidt et all , 120051) . 
For higher temperatures, the ground state seems to be 
a quadruplet, i.e. it can be described by a multi-orbital 
Anderson model with to = 2. The states building this 
quadruplet are obtained from spin-orbit coupled / states, 
which results in different ^-factors for its members. Inter- 
estingly, these different ^-factors in connection with the 
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small crystal-field splitting can lead to a behavior, where 
specific heat a nd susceptibility seem to have different low- 
energy scales (jScheidt et However, this dis- 
crepancy can be resolved by observing that the difference 
in g-factors leads to a "protracted" screening behavior for 
the specific heat in NRG calculations, while the system as 
a whole has only one, but strongly reduced Kondo tem- 
perature. In addition, NRG results for the low-energy 
spin dynamics show an anomalous energy dependence, 
which again is related to the diff erence in the q-factors 
enter ing the local susceptibility ( Anders and PruschkeL 
120061) . 

Conventionally, the Hund's exchange J appearing in 
(1 1 35[) is positive, i.e. mediating a ferromagnetic inter- 
action leading to Hund's first rule. However, there 
may be circumstances, for example a coupling to vibra- 
tional modes, Jahn- Teller distortions or crystal- field in- 
duced anisotropics ( De Leo and Fabriziol . 12004 ) . which 
can lead to an effective J < 0, i.e. antiferromagnetic 
exchange. In this case we encounter a situation simi- 
lar to the multi-impurity problem (see section HV.C.ip , 
where the exchange was generated by the RKKY effect. 
Under special conditions, an antiferromagnetic exchange 
then could lead to a quantum phase transition between 
Kondo screening and non-local singlet formation. Con- 
sequently, one may expect a similar transition for the 
multi-orbital model, too, when one varies J from the 
ferromagnetic to the antiferromagnetic re gime. Such a 
model was studied bv lFabrizio et al\ (120031) and, with an 
additi onal single-ion anisotropy. bv lDe Leo and Fabriziol 
(2004) for two orbitals, i.e. m = 2. The model indeed 
shows the anticipated quantum phase transition, which 
in this case is driven by the competition between the 
local antiferromagnetic exchange coupling and the hy- 
bridization to the band states. Furthermore, one can 
study the d evelopment of the spe c tral fu nction across this 
transition ( De Leo and Fabriziol . 12004 ) . One finds that 
the impurity spectral function on the Kondo screened 
side of this transition shows a narrow Kondo peak on 
to p of a broader reson a nce. As ha s also been observed 
by Pruschke and Bullal ( 20051 ) and iPeters and Pruschkd 
( 20061 ) this broad resonance is related to the exchange 
splitting J. The narrow peak transforms into a pseudo- 
gap on the unscreened side of the transition. 

iDe Leo and Fabrizid (|2005h have demonstrated that 
NRG calculations are possible even for m = 3, using 
the symmetries of the model to reduce the size of the 
Hilbcrt space blocks. They studied a realistic model for 
a doped Ceo molecule, taking into account the orbitally 
threefold degenerate t2 U lowest unoccupied molecular or- 
bitals. Again, coupling to vibrational modes can lead 
to Hund's coupling with negative sign. In this regime, 
one observes non-Fermi liquid behavior for half-filling 
n = 3, associated with a three-channel, 5 = 1 over- 
screened Kondo model. Interestingly, the critical sus- 
ceptibilities associated with this non-Fermi liquid appear 
to be a pairing in the spin and orbital singlet channel 
(|De Leo and Fabriziol [2005) . Using conformal field the- 



ory, the authors could deduce the residual entropy as 



S(T = 0) = -ln 



VE + i 

s/l- 1 



and also corresponding fractional values for the local 
spectral function p/(0) at the Fermi energy, which leads 
to non-unitary values in the conductance, and a non- 
integer power law for p/(uj) — p/(0). Away from half 
filling a quantum phase transition between Kondo screen- 
ing a nd local singlet formation a s function of filling oc- 
curs (|De Leo and Fabriziol [20051) . Multi-orbital models 
also arise in the context of quantum dots as described in 
Sec. IIV.A.3I and Sec. IIV.C.1I 



3. Selected results on low-energy properties 

In the following we shall present some selected results 
for the properties of the multi-orbital Anderson model 
(|135p with orbital degeneracy m — 2. 



a. Effect of Hund's coupling As a first example we want 
to discuss the influence of Hund's exchange J on the low- 
energy properties of the model Eq. (|135[) . Wc consider 
only J > 0, i.e. the usual atomic ferromagnetic exchange, 
and in addition use the constraint U' = U — 2,1 . For the 
conduction electrons we assume a band with a flat density 
of states pi 0) (e) = N F Q(D- |e|) and use A = t:V 2 N f as 
energy scale. Results for thermodynamic quantities are 
shown in Fig. [2U The calculations are for U / An = 167r 
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FIG. 21 Effective local moment T • x(T) (upper panel) and 
entropy (lower panel) for a two-orbital impurity Anderson 
model. Model parameters are U/Ao = 16tv at particle- hole 
symmetry. For comparison the results for one orbital are in- 
cluded (triangles). 

and particle-hole symmetry. NRG parameters are A = 5 
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and 2500 retained states per iteration. The triangles are 
results from a calculation with m = 1 for the same pa- 
rameters. For J — (circles) we find the behavior ex- 
pected for a SU{4) symmetry, i.e. an energy scale, T^ 1 , 



that is related to the Kondo scale at 



= 1, TV 



by 



= (T^) 1 /" 1 . However, even a small Hund coupling 
,/ = £7/100 leads to a dramatic reduction of Tk, which 
increases with increasing J. Let us note here, that for 
J > the effect of the orbital exchange term in (|135p is 
negligible, i.e. the results are indistinguishable if one does 
the calculation with and without this contribution. Very 
often, the orbital exchange term is neglecte d in theoreti- 
cal st udies of transition metal compounds ( Imada et all 
1998), an approximation which is supported by the above 
result. 



b. Crystal-field effects Recently, experiments showing 
unusual specific heat, magnetic susceptibility and re- 
sistivity data for Cei_ a; La 3 ;Ni9Ge4 have drawn a lot of 
attention, because this material has the "largest ever 
recorded value of the ele c tronic specific heat at low 
temperature" (|Killer et all l2004h of j(T) = AC/T w 



5JK _2 mol _1 . While the 7 coefficient continues to rise 
at the lowest experimentally accessible temperature, the 
magnetic susceptibility tends to saturate at low temper- 
atures. 

One possib l e sce nario ( Anders and Pruschkei l2006t 
IScheidt et all 120051) to account for the behavior of 
Cei_ a; La : ];NigGe4 is a competition of Kondo and crystal- 
field effects which leads to a crossover regime connecting 
incoherent spin scattering at high temperatures and a 
conventional strong-coupling Fermi liquid regime at tem- 
peratures much lower than the experimentally accessi- 
ble 30mK. The Hund's rule ground state of Ce 3+ with 
7 = 5 /2 is split in a tetragonal symmetry (jKiller et ali 
I2004D in three Kramers doublets. If the crystalline elec- 
tric field (CEF) parameters are close to those of cubic 
symmetry, the two low- lying doublets and , orig- 
inating from the splitting of the lowest r§ quartet, are 
well separated from the higher lying Tg doublet. Ignor- 
ing this Tq doublet, we can discuss two extreme limits. 
In a cubic environment, the CEF splitting vanishes and 
the low-temperature physics is determined by an SU (4) 
Anderson model. In a strongly tetragonally distorted 
crystal, on the other hand, the crystal-field splitting of 
the quartets is expected to be large. In this case, the 
low-temperature properties are determined by an SU (2) 
Anderson model. If, however, the material parame- 
ters lie in the crossover regime where the effective low- 
temperature scale T* is of the order of the crystal-field 

splitting <5cef = -ELm — E r (i), then the excited doublet 
7 7 . 

will have significant weight in the ground state so that 
the total magnetic response differs from a simple SU (4) 
Anderson model. 

Such a situation can be captured by a SU (4) Ander- 
son model with infinite U whose Hamiltonian is given by 
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FIG. 22 Comparison between y(T) = C(T)/T versus T and 

the susceptibility contributions of the two doublets versus T 

for £L ( i)/Ao = -8-5, <W/A = 0.015 and gl/gl = 2. The 
1 7 

contribution of the lower doublet, xi is much larger than the 
one of the upper doublet, \2- NRG parameters are A = 4 and 
1500 states kept in each step. 



(IScheidt et ali . l2005h 
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where \aa) represents the state with spin a and en- 
ergy E aa on the Ce 4/ shell, and Ckaa annihilates a cor- 
responding conduction electron state with energy tkaa- 
Note that locally only fluctuations between an empty and 
a singly occupied Ce 4/ shell are allowed. 

While the entropy and specific heat for the model Q136P 
can be calculated in the usual way, the Ce contribution to 
the susceptibility requires some more thought, because 
the spin-orbit coupled states have different ^-factors, 
which we label by g a . Thus , the total susceptibility is 
given by (|Scheidt et all , 120051 ) 
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While the g-factors are, in principle, determined by the 
CEF states of the multiplets, we view them as adjustable 
parameters and fix them together with E ( a ) by compar- 



ing with experiment ( Scheidt et all [20051) . 

The comparison between the temperature dependence 
of j(T) and x(T) is shown in Fig. [2U assuming a ratio 
of go/qi = 2 for a good fit to the experimental data 
(IScheidt et all l2005h . The ground state doublet dom- 
inates the magnetic response at low temperature and 
tends to saturate at temperatures higher than the 7- 
coeffi cient, consistent with the experiments ( Killer et all 
2004). We find this behavior only for CEF-splittings 
Sqef ~ T*(<5cef) while for much larger or much smaller 
values x{T) and 7(T) saturate simultaneously. 
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E. Bosonic degrees of freedom and dissipation 

In the models discussed so far, the bath consists of non- 
interacting fermionic degrees of freedom while the impu- 
rity is either represented by a fermion or a spin. This 
section deals with quantum impurity systems involving 
bosonic degrees of freedom. We distinguish between mod- 
els in which only a small number of bosonic degrees of 
freedom couples to the impurity and models where the 
impurity couples to a bosonic bath (corresponding to an 
infinite number of bosonic degrees of freedom) . 

As we will see below, the first case can be dealt with 
in the usual scheme, provided the subsystem consisting 
of impurity and bosons can be treated as a large impu- 
rity which is then coupled to the fermions. The second 
class, however, requires a different set-up for the NRG- 
procedure. 

Let us first consider the so-ca lled Anderson-Holstein 
model (jHewson and Meverl . 120021 ) in which the impurity 
is linearly coupled to a single bosonic degree of freedom 
(typically a phonon mode): 

H = H S1AM + A(6 f + b) .fife + <#tfb , (138) 

with -ffsiAM the Hamiltonian of the single-impurity An- 
derson model as in Eq. ^j. The coupling to the bosonic 
operators (tf and b) does not influence the mapping of 
the conduction electron part of the Hamiltonian to a 
semi-infinite chain. This means that the bosons enter 
the iterative diagonalization only in the very first step 
in which the coupled impurity-boson subsystem has to 
be diagonalized. The fact that only a limited number 
of bosonic states rib can be taken into account in this 
diagonalization imposes some restrictions on the param- 
eters A (the electron-phonon coupling strength) and wq 
(the frequency of the ph onon mode). As discussed in 
iHewson and Meverl ( 2002t h it should be sufficient to in- 
clude a number of rib ~ 4A 2 /wg bosonic states for the 
initial diagonalization. With an upper limit of n\> ~ 1000 
this means that the limit ljq — > (with fixed A) cannot 
be treated within this setup. 

Apa rt from this minor restriction, IHewson and Meverl 
(|2002f ) showed that the NRG (which is non-perturbative 
in both A and U) works very well for such type of im- 
purity models. In particular, both electron and phonon 
spectral functions as well as dynamic charge and spin 
susceptibilities can be calculate d with a hig h accu racy. 
As discu s sed in detail in iJeon et Mi 1)20031) and 



IChoi et all ( 2003al) . the calculation of the phonon 
spectral function needs some extra care and the au- 
thors introduced an impr oved method (as compared to 
IHewson and Meverl (|2002l )). The proper calculation of 
the phonon spectral function is important to discuss the 
softening of the phonon mode, see also the discussion in 
Sec. IV. CI in the context of lattice models with coupling 
to phonons. 

The low-energy features of the model Eq. (|138|) can 



be partly explained by an effective single-impurity An- 
derson model in which the coupling to the phonons is 
included in an effective interaction U c s- An explicit form 
of this interaction can only be given in the limit luq — > oo : 
U cS = U - 2X 2 /uj . Interestingly, IHewson et all ([20041 ) 
have shown that an effective quasiparticle interaction can 
be defined for any value of too. This can be accomplished 
with the renormalized perturbation theory by fitting the 
lowest lying energy levels obtained in the NRG calcula- 
tions to those from a renormalized Anderson model. 

These investigations represent a starting point for var- 
ious applications of the NRG to coupled electron-phonon 
system. For the investigation of transport properties 
of single molecule devices, for which the coupling to 
local phonons is a natural ingredient, similar models 
as Eq. CE3E 



(2004) and 



, have been investigated in Cornaglia et all 
Cornaglia and Grempell ( 2005af h Not only 



the coupling to the electron density as in Eq. (|138p . but 
also the change of the hybridization between molecule 
and leads due to the phonons has been shown to be im- 
porta nt for the conductance properties ( Cornaglia et all . 
120051) . 

Different physical phenomena can be expected in 
multi-orbital systems when the impurity degrees of free- 
dom couple to Jahn- Teller phonons. Such a model has 
been investigated in lHottal (J2006) and it was argued that 
within this model a new mechanism of Kondo phenomena 
with non-magnetic origin can be established. 

Two different strategies have been developed to study 
impurity models with a coupling to a bosonic bath, i.e. a 
bosonic environment with a continuous spectral density 
J(ijj). Let us discuss these strategies in the context of 
the spin-boson model 



H = -—<t x + 2<rz + 2^ 



(139) 

This model naturall y arises in the descrip tion of quantum 
dissipative systems ( Leggett et all 119871 ). The dynamics 
of the two-state system, represented by the Pauli ma- 
trices o~ Xt z, is governed by the competition between the 
tunneling term A and the friction term A; (a.; + aj). The 
ctj constitute a bath of harmonic oscillators responsible 
for the damping, characterized by the bath spectral func- 
tion 



J (w) 



UJ, 



(140) 



A standard parametrization of this spectral density is 

J(u) = 2iraul~ s lo s , < oj < uj c , s > -1 . (141) 

The case s = 1, known as Ohmic dissipation, allows a 
mapping of the spin-boson model onto the anisotropic 
Kondo model (for the definition of the Hamiltonian and 
the relation of its parameters t o tho se of the spin-boson 
model, see ICosti an d Kicffcr (1996)). The first strat- 
egy is then to apply the NRG to the anisotropic Kondo 
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model and to treat the fermionic conduction band in the 
usual way. Though restricted to the Ohmic case, such 
calculations have been shown to give very accurate re- 
sults for dynamic and thermodynamic quantities of the 
corresponding spin-boson and related models, as briefly 
discussed in the following. 

The focus of ICosti and Kiefferi (jl996h has been the 
equilibrium dynamics of the spin-boson model, in partic- 
ular the spin-spin correlation function, = ((c 2 ; o~ z )), 
at temperature T = 0. An interesting finding of this 
study, is that the spin relaxation function, x"(k-0/^, ex- 
hibits a crossover from inelastic to quasielastic behavior 
as a exceeds the value 1/3 signalling the onset to inco- 
heren t dynamics which occurs at a > 1/2 (|Leggett et all , 
1987). The accuracy of this approach has been shown via 
the comparison to exactly solvable limiting cases, such 
as the Toulouse point a = 0.5, and via the generalized 
Shiba relation. The issue of scaling and universality, 
concepts which are quite naturally connected to renor- 
malization group treatments of the Kondo problem, have 
been d iscuss ed in the context of the spin-boson model in 
ICostil (|l998h . Universal scaling functions have been cal- 
culated for thermodynamic quantities (the specific heat), 
the static susceptibility and the spin-relaxation function. 
Scaling as a function of temperature or frequency has 
been observed in the limit A — -> for fixed coupling 
strength a, this means that the scaling functions turn 

out to depend on the v alue of a. This is illustrated, 

for example, in Fig. 2 in ICostil ( 1998h which shows the 
temperature dependence of the specific heat for differ- 
ent a. In particular, a signature distinguishing weakly 
dissipative from strongly dissipative systems is found in 
7(T) — C(T)/T, which is found to exhibit a peak for 
a < 1/3 but is monotonic in temperature for a > 1/3. 
This, together with the above mentioned behavior in the 
spin-relaxation fu nction, is reminiscent of m easurements 
on Cei-^La^Akj ( Goremvchkin et all , \200^ i, but the ap- 
plicabil ity of an anisotrop ic Kondo model here is contro- 
versial ( Pietri et aLl . l200il ). Specific heats have also been 
calculated for more complicated tunneling models, such 
as th e ionic tunneling model w ith a spinless fermionic 
bath ( Ferreira and Libe"rol . l2000l ). which shows similar be- 
havior to the Ohmic spin boson model for a < 1 /4, and to 
an extension of this includ ing an assisted tunneling term 
( Ramos and Liberol . l2006h . The sp ectral density of the 
former has also been investigated ( Libero and Oliveiral 
Il990bh . 

The mapping of a bosonic bath to a fermionic one 
has been also expl o ited for various other problems. 
ICosti and McKenziei (|2003l ) have used the NRG to cal- 
culate the entropy of entanglement for the spin-boson 
model, a quantity which measures the entanglement be- 
tween the spin and the environment. Interestingly, the 
entanglement appears to be highest for a — > l - , where 
the system undergoes a quantum phase transition from 
a delocalized to a localized phase. 

The case of two bosonic baths which couple to differ- 
ent components of the impurity spin operator, has been 



discus sed in ICastro Neto et all ( 2003f) and iNovais et all 
(2005). The bosonic baths in these models can be 
mapped onto two independent fermionic baths, and a 
generalization of the anisotropic Kondo model is ob- 
tained. These models are of interest to study the effect 
of quantum frustration of decoherence. 

The second strategy to investigate impurity models 
with a coupling to a bosonic bath does not rely on a 
mapping to a fermionic impurity model. This approach 
- which has been te r med " bosonic" NRG - was intro- 
duced in Bulla et al. (2003) in the context of the spin- 
boson model (for full details see (|Bulla et a/.l .[2005)). Let 
us just mention briefly the main differences to the stan- 
dard (fermionic) NRG: the logarithmic discretization is 
now directly performed for the bosonic bath (for the spin- 
boson model, the bath spectral function J(uj) Eq. (|140|) is 
discretized); the subsequent mapping onto a semi- infinite 
chain is technically very similar to the fermionic case but 
the resulting tight-binding chain is built up of bosonic 
sites. This gives rise to additional difficulties in setting up 
the iterative diagonalization scheme because only a finite 
number of bosonic states JVb can be taken into account 
when adding one site to the chain. Furthermore, the set 
of iVb states should in general be optimized to give the 
best description of t he lowest lying ma ny-particle states, 
see the discussion in iBulla et all (|2005h . 

The first applications of the bosonic NRG focused on 
the spin-boson model Eq. Q139p , in particular the sub- 
Ohmic case with exponents < s < 1 in the parametriza- 
tion of the bath spectral function J(u>) Eq. (|14ip . The 
sub-Ohmic case does not allow for the mapping to a 
fermionic impurity model, in contrast to the Ohmic case. 
Furthermore, the bosonic NRG turns out to have cer- 
tain advantages over o ther approaches usua lly applied to 
the spin-boson model ( Leggett et all . [l987h as it is non- 
perturbative in both a and A. As an example of the 
success of the bosonic NRG we show in Fig.[^2]thc T = 
phase diagram of the spin-boson model with bias e = 
in the a-s plane for different values of the tunneling am- 
plitude A. 

The most remarkable feature of the phase diagram is 
the line of quantum critical points for < s < 1 which 
terminates for s — > l - in the Kosterlitz-Thouless transi- 
tion of the Ohmic case. The critical e xponents along this 
line have been discussed in detail in IVoita et all (|2005h 
where it was shown that the exponents fulfill certain hy- 
perscaling relations. 

There are still quite a number of open issues and pos- 
sible applications in the context of the sub-Ohmic spin- 
boson model which ca n be investigated with t he bosonic 
NRG, see for example iTong and Vo jta| (12006ft . In addi- 
tion, the NRG can be generalized to more complex im- 
purity models with a coupling to a bosonic bath. Two 
recent exampl es are the investigations of th e Bose-Fermi 
Kondo model ( Glossop and Ingersentl 120051 ) and of mod- 
els which ar e connected to electro n and exciton transfer 
phenomena ( Tornow et all , l2006al ibh . 

Technically the most challenging of these extensions 
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FIG. 23 Phase diagram of the spin-boson model for T — 
calculated with the bosonic NRG for various values of A. 
A line of quantum critical points separates the delocalized 
from the localized phase. (Figure adapted from iBulla et all 
(120031 )'). 



is the study of impurity models which couple to both 
fermionic and bosonic baths. The first successful treat- 
ment of such a model, the Bose-Fermi Kondo model with 
Ising-type couplin g between spin and bo s onic bath, ha s 
been discussed in iGlossop and Ingersentl (|2005l l2006afh 
Here the two baths are mapped on two semi-infinite 
chains, one for the bosonic and one for the fermionic 
degrees of freedom. Due to the competition between dis- 
sipation and screening, quantum phase transitions occur 
which turn out to be in the same universality classes as 
the transitions in the spin-boson model. This can be un- 
derstood through a mapping between the (Ising-)Bose- 
Fermi Kondo model and the spin-boson model where the 
bath spectral function has both an Ohmic and a sub- 
Ohmic component (the Ohmic part represents the cou- 
pling to the fermionic bath, while the sub-Ohmic part 
is the same as the one in the original model, see also 
iLi et all (|2005l )). On the other hand, the (Ising-)Bose- 
Fermi Kondo model with an Ohmic bosonic bath can 
also be ma pped onto the aniso tropic Kondo model, see 
for example iBorda et all (|2005t ) . 



V. APPLICATION TO LATTICE MODELS WITHIN 
DMFT 

The application of the NRG is restricted to quan- 
tum impurity systems with the impurity degree of free- 
dom coupled to a non-interacting bath. Therefore, the 
NRG cannot be directly applied to lattice models of 
interacting particles, such as the Hubbard model (see 
Eq. (fT42| in Sec. IV. Ap . Early attempts to extend Wil- 
son's concepts to su c h mo d els f a iled ( Brav and Chuil . 
119791 : IChui and Bravl . 119781 : Heel . Il979h . The reason 
for this failure was later found to be connected with 
bounda ry conditions between "syst e m" and "environ- 
ment " ( Noack and Manmanal . 120051 : IWhite and Noackl 
I1992D . and led to a novel scheme now adays known as 
density-matrix renormalization group (|Hallberd . 120061: 



ISchollw5ckl I2005T ). which today is a standard technique 
to study one-dimensional interacting quantum models. 

There exists, however, an approximation for correlated 
lattice models, where the interacting lattice problem is 
mapped onto an effective quantum impurity model, for 
which the NRG can be applied. Underlying this ap- 
proach is the dynamical mean- field theory (DMFT). The 
DMFT has its origin in the investigation of correlated 
lattic e models in the limit of infinit e spatial dimension- 
ality (|Metzner and VollhardtJ . H989). A proper scaling 
of the hopping matrix element t in models such as the 
Hubbard model (] 142[) leads to a vanishing of all non- 
local self-energy sceleton diagrams. The resulting purely 
local self-energy S(z) can be identified with the self- 
energy of an effective single-impurity Anderson model. 
In this sense we speak of a mapping of a lattice model 
onto an effective quantum impurity model, typically the 
single-impurity Anderson model as introduced in Eq. ^ , 
supplemented by a self-consistency condition, which de- 
termines the bath degrees of freedom of the effective 
quantum-impurity. Since the technical details of the 
DMFT are not the sub ject of this re v iew, w e refer the 
reader to the review by I Georges et all ( 19961 ). 

To investigate lattice models in the DMFT we there- 
fore need a technique (analytical or numerical) to cal- 
culate the full frequency dependence of the self-energy 
for a single-impurity Anderson model defined by ar- 
bitrary input parameters (s{, U, T, and a manifestly 
energy-dependent hybridization function A(w)). There 
are many methods besides the NRG available to cal- 
culate dynamic quantities for quantum impurity models 
and we shall not give a n overview here (for reviews see, 
for example , Sec. VI in Georges et all (|l996l ). Sec. Ill in 
iMaier et all (|2005l ). and bullal (|2006l )l. but rather concen- 
trate on the application of the NRG method to the Hub- 
bard model (see Sec. IV. A[) . the periodic Anderson and 
Kondo lattice models (see Sec. IV. Bl) . and lattice mod- 
els with coupling to phonons (see Sec. IV. CP within the 
DMFT approach. 

Before we discuss the results obtained for those mod- 
els, let us comment on pecularities of the NRG when ap- 
plied to DMFT calculations. The DMFT self-consistency 
specifies at each iteration an input hybridization function 
A(tt>), the form of which depends on the model under in- 
vestigation, its parameters and also on the history of the 
previous DMFT-iterations. The frequency dependence of 
A(w) has to be taken into account within the logarith- 
mic discretization scheme, exactly as described in Sec. |TT] 
and already employed in the NRG investigations of the 
soft-gap Anderson model, see Sec. IIV.C.2I 

Concerning the output, the quantity of interest is usu- 
ally the self-energy Sam of the effective single-impurity 
Anderson model, although in some cases, as for the Bethe 
lattice, the knowledge of the single-parti cle Green func- 
tion is sufficient for the DMFT iteration ( Georges et all , 
119961) . It has proven advantageous to calculate, within 
DMFT, the self-energy Sam via the ratio of a two- 
particle and a one-particle Green function, see Eq. 
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As discussed in lBulla et ail (|l998l ) (see also Sec. IIII.B.2|) 
the calculation of the self-energy via Eq. (|55j) significantly 
improves the quality of the results. This approach has 
been used in most NRG calculations within DMFT. 



A. Hubbard model 

The simplest model for correlated fermions on a lattice 
is the single-band Hubbard model with the Hamiltonian 



transition, when the quasiparticle peak has disappeared, 
an insulating solution with a preformed gap is realized. 
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Consequently the first applications of the NRG within 
DMFT focused on this model; in particular on the Mott- 
transition, which the Hubbard model displays in the half- 
filled paramagnetic case. These investigations and fur- 
ther generalizations are described in the following sub- 
sections. 



1. Mott metal-insulator transition 

Although the qualitative features of the Mott transi- 
tion have been correctly described very ear ly in the devel- 
opmen t of the DMFT (see the review by iGeorges et all 
(1996)), the NRG helped to clarify a number of conflict- 
ing stateme nts i n the litera t ure (s ee the discussion in 
IBullal I)1999T ) and lBulla et~al\ (|200lh h The NRG method 
appears to be ideally suited to investigate the Mott tran- 
sition because (i) the transition occurs at interaction 
strengths of the order of the bandwidth which requires 
the use of a non-perturbative method, and (ii) at T = 
the Mott transition is characterized by a vanishing en- 
ergy scale, T* — > 0, when approached from the metallic 
side. Thus, a method is needed that is able to resolve 
arbitrarily small energies close to the Fermi level. 

The first investigation of the Mott transition w i th the 
NRG has been performed bylSakai and Kuramotol ( 1994( ) 
(see also IShimizu and Sakail (fl995h h These calculations 
did not use an expression of the self-energy as in Eq. ([55]). 
but nevertheless a Mott transition and a hysteresis region 
have been observed, with critical v alues very close to the 
ones reported later in lBullal ( 19991 ). 

A detailed discussion of t he NRG calc ulations for the 
Hubbard model is given in IBullal (|l999l) for T = and 
in lBulla et all ( 200lh for finite temperatures. The main 
results are summarized in Fig. 1241 Spectral functions cal- 
culated with NRG for the half-filled Hubbard model in 
the paramagnetic regime for different values of U and 
T = are shown in Fig. IMk . Upon increasing U from 
the metallic side, the typical three-peak structure forms, 
with upper and lower Hubbard peaks at w w ±{7 /2 and 
a central quasiparticle peak at lo = 0. The width of 
this quasiparticle peak goes to zero when approaching 
the transition from below, U /* XJ C ^ ~ 1A7W, with W 
the bandwidth of the noninteracting model. Right at the 
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FIG. 24 (a) Spectral functions for the half-filled Hubbard 
model at T = for various values of U (similar data as in 
Fig. 2 in IBullal (|l999l ): (b) phase diagram for the Mott tran- 
sition. For a com parison with other methods see the corre- 
sponding Fig. 9 in lBulla et all (|200lh . 

We should remind the reader at that point that the 
NRG results for dynamic quantities have a certain fixed 
resolution on a logarithmic scale (c.f. the discussion in 
Sec. IIII.BI) . This means that structures close to oj = 
are much better resolved than structures at, for exam- 
ple, the band edges of the Hubbard bands. In contrast, 
the dynamical density-matrix renormalization group re- 
cently applied to the DMFT for the Hubbard model 
works w ith a fixed resoluti on on a linear scale, see for 
example iKarski et al. ( 2005). The structures close to the 
inner band edges of the Hubbard bands seen in these cal- 
culations (see Fig. 2 in IKarski et al\ ( 2005D ) cannot be 
resolved in present implementations of the NRG method. 

Figure l2~4"b shows the T-U phase diagram for the Mott 
transition, again only consi dering paramagnetic phases. 
As already observed earlier ( Georges et all fl996h . there 
are two transition lines, because the insulator-to- metal 
transition occurs at a lower critical value (C/ Ci i(T)) as the 
metal-to-insulator transition (U C ,2(T)). Within this hys- 
teresis region, both metallic and insulating solutions can 
be stabilized within the DMFT self-consistency. With 
increasing temperature, the hysteresis region shrinks 
to zero at a critical T c , above which there is only a 
crossover from metallic-like to insulating-like solutions; 
this crossover region is indicated by the dashed lines in 
Fig. [Mb. The NRG values for U cA/2 (T) have later been 
verified by a number of other, non - pertu rbat ive tech- 
niques (see for example iTong et all (I2001D or Potthoff 
(l2003h h 

As discussed above, most of the NRG calculations 
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within DMFT have been performed using Eq. (|85|) for 
the calculation of the self-energy. This quantity itself 
shows intere sting propertie s (see, for example , Fig. 3 in 
iBullal (|l999l) and Fig. 5 in iBulla et all (|200ll) h and al- 
lows to calculate th e {/-depende nce of the quasiparticle 
weight, see Fig. 1 in|Bulla| (|l999h . 

The Mott transition can also be induced by moving 
away from half-filling, provided the value of U is larger 
than the U c for the half-filled case. Unfortunately, no 
systematic NRG calculations have been published for this 
filling-induced Mott transition, despite the fact that the 
NRG can be easily extended to the Hubbard model away 
from particle-h ole symmet r y. On ly a few results for the 
phase diagram ( Ono et al\, 200lh and spectra l funct ions 
(jFreericks et all l2003t iKrug von Nidda et all |2003| ) are 
available in the literature. 

A nice feature of the DMFT is that it also allows for the 
calculation of physical quantities other than the single- 
particle Green function, in particular susceptibilities and 
also transport properties, both static and dynamic. This 
aspect of the DMFT has been intensively used already 
in the early application s (se e for example the r eviews 
bv lPruschke et all ( 19951 ) and lGeorges et al\ (1996)), cm- 
ploying different methods to solve the effective quantum- 
impurity problem. Howev er, apart from discuss i ons o f 



large U, ferromagnetic solutions can be stabilized. For 
intermediate U and finite doping, magnetic ordering ap- 



i i i 1 i 

the Ai g Raman response (|Freericks et all 12 001. 2003 1 



and calculations of th e resistivity ((Georges et all 12004 ; 
iLimelette et al , 2003 ) and the lo cal dynamic susceptibil- 
ity ( Krug von Nidda et all [20031 ). no detailed studies of 
such quantities for the paramagnetic phase of the Hub- 
bard model have been performed yet with the NRG. 



2. Ordering phenomena 

The Mott transition from a paramagnetic metal to a 
paramagnetic insulator is merely one of the many fea- 
tures in the rich phase diagram of the Hubbard model 
and its generalizations. In addition, various types of or- 
dering phenomena occur, such as charge, orbital (in case 
of multi-orbital models), and magnetic ordering, and - 
possibly - superconductivity. The NRG has been used 
in particular to study magnetic ordering phenomena in a 
wide range of parameters. 

For the investigation of symmetry broken phases 
within DMFT, the self- consistency equations have to be 
adapted appropriately ( Georges et all [T996). The ef- 
fective impurity models still have the structure of the 
single-impurity Anderson model so that the application 
of the NRG is str aightforward, see the discussion in 
IZitzler et al\ (|2002f) . This work also contains a detailed 
study of the magnetic phases of the Hubbard model at 
T = both at and away from half-filling. Right at half- 
filling and for a particlc-holc symmetric band-structure, 
the groundstate is always antiferromagnetically ordered. 
Upon doping, the situation is more complicated as shown 
in Fig.(2U For small values of U, phase separation within 
the antiferromagnetic phase is observed, while for very 
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FIG. 25 Ground state magnetic phase diagram for the Hub- 
bard model on a hypercubic lattice. The phases show anti- 
ferromagnetic (AFM) order, which for smaller U and doping 
8 > also shows phase separation (PS), and ferromagnetic 
order at large U. To display the whole interval [0, oo), the 
vertical axis was r escaled as U/(l + U) (see also Fig. 10 in 
IZitzler et all l|2002l )). 

pears to exist but its type could not be determined yet. 

In contras t to t his work, the NRG calculations in 
IZitzler et all (2004) concentrated on the antiferromag- 
netic phase at half-filling in a Hubbard model with frus- 
tration. As expected, the antiferromagnetic region in the 
T-U phase diagram is suppressed upon increasing frus- 
tration. However, the resulting phase diagram turns out 
to be significantly different from the one proposed fo r 
the frustrated Hubbard model in Georg es et al. I (|1996T) . 
where it was claimed that the main effect of frustration is 
to suppress Tn such that the first-order Mott transition is 
visible above the antiferromagnetic region. This contro- 
versial issue certainly calls for more detailed calculations 
(with NRG and other methods); after all, the similarity 
of the phase diagram for the Hubbard model in DMFT 
and the experimental one for the transition metal oxide 
V2O3 h as been claimed to be one of the early success of 
DMFT (jGeorees et aZ.Ul996l ). 

The optical conductivity in the antiferromagnetic 
phase of the Hubbard model at half-filling and 

zero temperature h as been studied in detail in 

iPruschke and Zitzlerl (|2003h . For small values of U , the 
antiferromagnetic phase shows signatures of a Slater in- 
sulator, while for large U a Mott-Heisenberg picture ap- 
plies. There is a smooth crossover between these two 
limiting cases upon variation of U, in contrast to the 
Mott transition in the paramagnetic phase. The evi- 
dence from the optical data has been supported by a de- 
tailed discussion of the local dynamical magnetic suscep- 
tibility, giving additional insight into the subtle changes 
in the physics of charge- and spin-degrees of freedom 
across the Mott metal-insulator transition ( Pruschkel 
120051 : IPruschke and Zitz ieil l2003l ). 

When the Hubbard model Eq. (|142p is supplemented 
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by a nearest-neighbor Coulomb repulsion V, a transition 
to a charge ordered state is observed upon increa s ing V . 
This transition has been studied in IPiet ig et al 
for the quarter-filled case. NRG calculations, together 
with results from non-crossing approximation and exact 
diagonalization show a phase diagram with a reentrant 
charge ordering transition, a feature which has also been 
observed in a variety of transition metal oxides. The 
NRG results in this work are restricted to T — 0, where 
the transition is of first order. It would be very interest- 
ing to extend the NRG calculations to a wider range of 
parameters, in particular to finite temperatures to study 
the change of the character of the transition which is 
continuous at higher T. 



3. Multi-band Hubbard models 

The application of the NRG to the investigation of 
multi-band Hubbard models within DMFT is still in a 
very early stage. This is because (i) the computational 
effort grows considerably with the number of orbitals and 
(ii) the DMFT requires a very high accuracy for the cal- 
culated dynamic properties. Furthermore, self-consistent 
solutions of the DMFT equations have to be obtained. 

The first (and so far the only) DMFT results for a two- 
b and Hubbard model using th e NRG have been presented 
in lPruschke and Bulla! ([2005) . In this work, two different 
strategies have been used to handle the complexity of 
the problem. The first one is to explicitly include the 
orbital quantum numbers in the iterative construction of 
the basis states. As for the impurity models discussed in 
Sec. HV.Di this additional quantum number significantly 
reduces the typical matrix size. However, this approach 
fails as soon as the Hamiltonian contains terms which 
break the orbital symmetry. 

The second strategy is an asymmetric truncation 
scheme: Instead of adding both orbital degrees of free- 
dom simultaneously, the Hilbert space is truncated after 
adding each orbital individually, which also leads to a 
significant reduction of the typical matrix size. This ap- 
proach works quite well in a wide range of parameters, 
but it appears to violate the orbital symmetry, if present. 
However, in the presence of a crystal-field splitting of the 
orbitals, such a strategy might be usable . 

The focus in iPruschke and Bul la (2005) was on the role 
of the Hund exchange coupling J on the Mott transition 
in the two-band Hubbard model. It was found that both 
the position in parameter space and the nature of the 
Mott transition depend on the value of J and the precise 
form of the coupling. For example, the replacement of a 
rotationally invariant Hund exchange by an Ising-like ex- 
change leads to a significant change in the physics of the 
Mott transition. Note that such features can be partly 
understood already on the level of the corresponding ef- 
fective impurity models, which underlines the importance 
to thoroughly investigate the properties of the impurity 
models appearing in the DMFT. 



4. Other generalizations of the Hubbard model 

Let us conclude this section with a brief overview of 
applications of the NRG to various other generalizations 
of the Hubbard model. 

The influence of correlations in a conduction band 
(modeled by a Hubbard model within DMFT) on the 
physics of the single-impurity Anderso n model has been 
investigated in IHofstetter et all (|2000h . The DMFT ap- 
proach allows to map this model on an effective impu- 
rity model with two coupled correlated sites, the first 
one corresponding to the original impurity and the sec- 
ond one coming from the DMFT treatment of the Hub- 
bard model. This two-site cluster coupl es to a free effec- 
tive c onduction band. As discussed in IHofstetter et all 
(120001) . correlations in the conduction band have a sig- 
nificant influence on the low-energy scale and also lead 
to a suppression of the Kondo resonance. 

Within the so-called Anderson-Hubbard model, disor- 
der effects can be incorporated via a random distribution 
of the local energies £j. This model has been studied in 
iBvczuk et all (|2004h within DMFT for binary alloy dis- 
order. The application of the NRG here is standard - 
two independent single-impurity Anderson models have 
to be considered at each iteration step. Nevertheless, the 
physics of this model is already quite interesting, in par- 
ticular the occurrence of a Mott tra nsition at non -i nteger 
filling. The DMFT treatment in IBvczuk et all (|2004l ) 
does not allow for true Anderson localization (as far as 
disorder is concerned, the DMFT is equivalent to the co- 
herent potential approximation and the main effect of the 
binary disord er is to split the band s). This deficiency has 
been cured in lBvczuk et al. I (l2005h where a generalization 
of the DMFT approach has been used, based on the geo- 
metrically averaged (typical) local density of states. This 
allows to study both Mott insulating and Anderson in- 
sulating phases, see Fig. I in this paper. The calculation 
has been performed using a continuous probability distri- 
bution, approximated by up to 30 different values of £j, 
so that in each DMFT step a corresponding number of 
independent single-impurity Anderson models have to be 
considered. All the NRG calculations for the Anderson- 
Hubbard model have been so far restricted to T = and 
to phases without long-range order. 

Recently, the NRG has been used within an extension 
of the stand ard DMFT. The ' DMF T+1 V approach as in- 
trodu ced in lSadovskii et al. I (120051) and iKuchinskii et all 
(2005) adds to the local self-energy a fc-dependent part 
Efc. Applied to the one-band Hubbard model, the effec- 
tive single-impurity Anderson model is still of the same 
type as the one appearing in standard DMFT, the only 
difference is in the structure of the self-consistency equa- 
tions. 

For details of the NRG calculations and the dis- 
cussion of the physic s of t he Falikov-Kimball model 
(lAnders and CzvchollL 120051) and the ionic Hubbard 
model (jJabben et all . 120051) we refer the reader to the re- 
spective references. Both papers show the usefulness of 
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the NRG approach to a wide range of correlated lattice 
models within DMFT, in particular for the calculation of 
dynamic quantities at low temperatures. 



B. Periodic Anderson and Kondo lattice models 

A variety of Lanthanide- and Actinide-based com- 
pounds can be characterized as heavy fermion systems 
with a strongly enhanced effective mass of the quasiparti- 
cles. These compounds contain well localized 4/ or 5/ or- 
bitals coupling via a hybridization to a conduction band 
consisting of s, p or d orbitals. The appropriate micro- 
scopic model for these materials is the periodic Anderson 
model (PAM) 
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When charge fluctuations of the / orbitals are negligi- 
ble, the low energy physics of the PAM can equally be 
described by the Kondo lattice model 
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The large effective mass in these models arises from a 
strongly reduced lattice coherence scale To; this is one 
of the reasons why the NRG is very well suited for the 
investigation of heavy fermion behavior when the PAM 
or the Kondo lattice model are treated within DMFT. 
The main difference for the NRG treatment (as com- 
pared to the Hubbard model) lies in the DMFT self- 
consistency. This means that the structure of the ef- 
fective impurity model is changed only via the effective 
hybridization A(w) (which may lead to complications as 
discussed later). 

The PAM with on-site hybridization Vij — V5ij and 
particle-hole symmetry on a hyp e rcubi c lattice has been 
discussed in Shi mizu and Sakai (|!995t) . iPruschke et all 
(|2000h . and IShimizu et all (|2000t ). In this case, a hy- 
bridization gap at the Fermi level appears in the spectral 
functions for both conduction and / electrons. This effect 
for the / spectral function is shown in Fig. [2U by the full 
lines in the main panel and left inset. Apparently, the 
Kondo resonance of the corresponding single-impurity 
Anderson model (dashed curves in the main panel and 
left inset in Fig. I26p . for which the hybridization func- 
tion is given by the bare density of states of the lattice, 
is split in the periodic model. The energy scale of the 
gap in the PAM (proportional to the lattice coherence 
scale To) depends exponentially on the model parame- 
ters, similar as for the width of the Kondo resonance in 
the impurity model (proportional to the Kondo temper- 
ature Tk). Further analysis shows that the lattice coher- 
ence scale Tp is enhanced over th e impurity scale Tk (for 
details see IPruschke efal\ (|2000f) ). 




FIG. 26 Comparison of dynamic properties for the particle- 
hole symmetric periodic Anderson model (solid lines) and 
the corresponding single-impurity Anderson model (dashed 
lines). Main panel: / electron spectral function; left inset: 
enlarged view of the region around the Fermi energy; right in- 
set: (effective) hybrid ization function. (Figure adapted from 
IPruschke etail (|2000l )). 



The right inset to Fig.l2l)lshows the hybridization func- 
tion of the effective impurity model after self-consistency 
has been reached (full line) in comparison to the same 
quantity entering the isolated impurity (dashed line). At 
first sight, the only difference seems to be the gap at the 
Fermi level. However, for the particle-hole symmetric 
case, one can show that A(w) has a pole at the Fermi 
level. At first sight, this pole appears to be a problem 
for the NRG as the logarithmic discretization explicitly 
excludes the point u> = 0, i.e. such a pole cannot be incor- 
porated in the mapping to the semi-infinite chain. The 
way out is to take the pole into account via an extra site 
which couples directly to the impurity, thus removing the 
pole from the hybridization function. 

Due to the appearance of the hybridization gap, the 
particle-hole symmetric PAM seems rather suitable to 
describe so-called Kondo insulators, but not the metallic 
heavy fermion behavior. There are various ways to drive 
the PAM into the metallic regime, two of which we will 
discuss in the following. One possibility is to use asym- 
metric parameters for the / electrons (e/7^— C//2) and 
to keep the conduction band symmetric so that n c i=a 1. 
More interesting (and also physically more relevant) is 
the opposite case, namely keeping ej = —U/2 and shift- 
ing the conduction band center-of-mass away from the 
Fermi level, so that n c is reduced from 1. This situation 
has been disc ussed in the con text of Nozieres exhaus- 
tion principle ( Nozieresl Il998l ) which states that, upon 
decreasing n c , there will not be enough conduction elec- 
trons available to screen the impurity spins. Collective 
screening then becomes effective only at a strongly re- 
duced lattice coherence scale. We do not want to go into 
the details h ere and refer the reader to the discussion 
in Sec. 5.4 of [Vidhvad hiraia and Loganl (|2004h , which is 
based on results for the PAM obtained with the local- 
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moment approach. 

There is one particular feature found in the DMFT cal- 
culations which at first sight seems to s upport Nozieres 
i dea, n amely that, as shown in Fig. 6b in lPruschke et all 
(2000), the effective hybridization function is strongly re- 
duced in a region close to the Fermi level. Since this 
quantity can be interpreted as being proportional to the 
conduction band density of states effectively seen by the 
/-states, it seems that indeed there are less conduction 
electrons available to screen the mom e nts of the / elec- 
trons. Figure 6a in iPruschke et~all (|2000h shows the 
corresponding /-spectral function, which consequently 
now displays metallic (heavy fermion) behavior. The 
corresponding lattice coherence s cale is now red u ced as 
compared to Tk, see Fig. 8 in IPruschke eial\ ()2000h 
where the dependence of both Tq and Tk on n c is plot- 
ted. However, in contrast to Nozieres original claim, i.e. 
To oc (Tk) 2 , a behavior To cx Tk is found, with a prefactor 
decreasing exponentially with decreasing n c . Again, the 
ability of the NRG to accurately identify exponentially 
small energy scales proves to be of great value here. 

Another route to metallic behavior in the PAM is to 
incorporate a dispersion of the / electrons of the form 



f j(jfii 



(145) 



<ij>,a 



The effect of such a dispersion term - in particular the 
closing of t he gap upon increasin g tf - has been studied 
in detail in [Sh lmizu et all f|2000h for both the particle- 
hole symmetric and asymmetric cases. The authors of 
this work also study the relation between charge and spin 
gaps in the dynamical susceptibilities and the hybridiza- 
tion gap in the spectral function. 

A metallic ground state of the particle-hole symmetric 
PAM can also be realized when the hybridization be- 
tween / electrons and c electrons is only between nearest 
neighbors: 



V : i,j nearest neighbors 
: otherwise 



(146) 



For T = 0, the PAM with nearest-neighbor hybridiza- 
tion shows a notable difference to the models discussed 
above: the low-energy scale T does no longer depend 
exponentially on U bu t vanishes at a finite critical U c 
(jHeld and Bull a. 2000). This behavior is reminiscent of 
the physics of the Mott-transition in the Hubbard model. 
The difference, however, is that in the PAM defined by 
(I143[) and (| 146[) the Mott-transition occurs only in the 
subsystem of the / electrons - a gap opens in the / 
electron spectral function while the c electron part still 
h as finite spectral weigh t at the Fermi level (see Fig. 3 
in iHeld and Bu lla ( 2000h ) so that the overall system re- 
mains metallic. 

The calculations described so far have all been re- 
stricted to T = 0. Finite temperature calculations for 
single-partic le and magnetic ex c itatio n spectra have been 
presented in lCosti and Maninil 12 002) for the Kondo lat- 
tice model. As in IPruschke et all ( 2000l) . one focus 



of these studies has been the variation of the spectra 
with conduction band filling n c . For this case it was 
found that the spectra exhibit two energy scales, one 
being the Kondo temperature Tk of the corresponding 
single-impurity Kondo model, the other one being the 
Fermi liquid coherence scale To which, for low carrier 
densities, n c <§; 1, is strongly reduced as compared to 
Tk, si milar to the observations made in IPruschke et all 
<|2000h for the PAM. A comparison of the temperature 
dependence of the photoemission spectra with experi- 
mental data on YbInCu4 showed good agreement. A 
ferromagnetic version of the Kondo lattice model with 
Coulo mb interactions in t he co nduction band was stud- 
ied bv lLiebsch and Costil (|2006h in the context of the or- 
bitally selective Mott phase of the two-band Hubbard 
model with inequivalent bands. The physics of this 
model is quite different to the us ual Kondo lattice mode l. 
In particular, one finds, as in iBiermann et al. (2005). 
non-Fermi liquid or bad metallic behavior, depending 
on whether the ferromagnetic exchange is isotropic or 
anisotropic, respectively. 

All these results have been obtained for the paramag- 
netic phases of the PAM or Kondo lattice model. Of 
course, the presence of localized moments implies the 
possibility for magnetic ordering, which is frequently ob- 
served in heavy fermion compounds, partly in close vicin- 
ity to superconducting phases. These issues have not 
been addressed yet with the NRG, but, as has been 
demonstrated for the Hubbard model, are of course acces- 
sible with this method and are surely a promising project 
for future NRG calculations for the PAM or Kondo lat- 
tice model within the DMFT. 

Magnetic quantum phase transitions in the Kondo lat- 
tice model have b een the focus o f calculations within the 
extended DMFT (|Si et aZ.L l200lh . for which the effective 
impurity model includes a coupling to both fermionic and 
bosonic baths. The NRG has been generalized to such 
types of impurity models, see Sec. IIV.E1 and recent ap- 
plications of the NRG within the e x tended DMFT are 
discus sed in IGlossop and IngersentI ( 2006bl ): IZhu et all 
(|2006l) . 



C. Lattice models with phonons 



Let us consider the Hubbard model Eq. (| 142[) and sup- 
plement it by a local coupling of the electron density 
to the displacement of Einstein phonons with frequency 
luq. This results in the Hubbard-Holstein model with the 
Hamiltonian 



ia <r/>tr i 

'Y,{ b l + b ) c l c - + ^H b l h ^ ( 14? ) 



-fj 



The limit U — > of this Hamiltonian gives the Holstein 
model, still a highly non-trivial model as discussed in the 
following. 
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Within DMFT, the model Eq. (fT47)) maps onto the 
Anderson-Holstein (impurity) model , to which the NRG 
metho d has been first applied by iHewson and Meyerl 
(2002), see Sec. HV.El The self-consistency equations are 
the same for both Hubbard- Holstein and the pure Hub- 
bard model, the only difference lies in the calculation of 
the self-energy for the effective impurity model which now 
contains an additional contribution from the coupling to 
the phonons. This contribution can also be calculated as 
a ratio of two correlation functions. 

From a technical point of view, there is no difference 
in the NRG-treatment of the Hubbard-Holstein model 
with either finite or zero value of U . Historically, the 
first applications of the NRG have been for the U = 0- 
case (the Holstein model , see iMever et all (|2002h and 
iMever and Hewsonl ( 20031 )1 and have already revealed a 
number of interesting results. As compared to other 
methods applied to the Holstein-model, the advantage 
of the NRG combined with the DMFT is that it is non- 
perturbative in both g and U and that it allows to study 
the case of a macroscopic electron density (in contrast to 
the few electron case). 

For the half-fi lled case, the investigations of 
IMever et al. I (|2002l) showed some unexpected properties 
for the transition from a metal to a bipolaronic insula- 
tor at a critical coupling g c . In contrast to the Mott 
transition in the Hubbard model, no hysteresis and no 
preformed gap is observed here (at least for small values 
of ujq) which indicates that the physics of the transition 
to the bipolaronic insulator might be completely different 
(whether it is connected to locally critical behavior is an 
interesting subject for future research). 

For large values of ujq, the physics of the transi- 
tion is getting closer to that of the Hubbard model 
( Mever and Hewsonl . 120031 ). This is because in the ujo~^ 
oo-limit the Holstein model can be mapped onto the at- 
tractive Hubbard model which has the same behavior as 
the repulsive Hubbard model when charge and spin de- 
grees of freedom are interchanged. 

The phase diagram of the Hubbard-Holstein model at 
half-filling, T = 0, and neglecting any lo ng-range or- 
dered p hases has been discussed i n detail in IKoller et all 
(l2004bl ). iJeon et all (|2004) . and iKoller et all (|2004afh 



The main features are summarized in Fig. [57]which shows 
the position of the phase boundaries between metallic, 
Mott-insulating, and bipolaronic insulating phases. The 
nature of these various transitions, together with the be- 
h avior of dynamic qua ntities has been discussed in detail 
m IKoller et all ((2004a') . Let us point out here the behav- 
ior of the phonon spectral function which shows a consid- 
erable phonon softening upon approaching the transition 
to the bipolaronic insulator. Such a softening is absent in 
the approach to the Mott insulator, for the simple reason 
that close to the Mott transition, where charge fluctua- 
tions are strongly suppressed, the phonons are effectively 
decoupled from the electrons. One of the interesting top- 
ics for future research are models which do show such a 
phonon softening even close to the Mott transition; this 



can possibly be accomplished by considering additional 
orbital degrees of freedom. 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0. 



bipolaronic ^ 




' / 

/ / metallic 


J Mott 






J insulator 







4 
U 



6 7 



FIG. 27 Phase diagram of the half-fi lled Hubbard-Holste in 
model for T = (figure adapted from IKoller et all (|2004al 1V 
The thick solid lines denote the phase boundaries while the 
dashed line corresponds to U c a = U— 2g 2 /u>o = 0. The dashed 
and do t-dashed lines are polaronic lines, see Koll er et all 
(I2004al 1. 

Recently, the Hubbard-Holstein model has been stud- 
i ed for the quarter -filled case and large values of U 
( Koller et al. I. l2005alh In this situation, a strongly renor- 
malized band of polaronic quasiparticle excitations oc- 
curs within the lower Hubbard band of the electronic 
spectral function. 

All the investigations of lattice models with electron- 
phonon coupling described in this section have been re- 
stricted to T = and to phases without long-range or- 
der. Generalizations to finite temperatures and ordered 
phases (such as superconducting and charge ordered 
phases) appear to be possible within the DMFT/NRG 
approach and will certainly give interesting results and 
new insights. Other possible generalizations are different 
types of lattice models, such as the periodic Anderson 
model and multi-orbital models, and models with a dif- 
ferent type of coupling between electrons and phonons. 



VI. SUMMARY 

Let us first summarize here the main purpose of this 
review: 

(i) To give a general introduction to the basic concepts 
of the numerical renormalization group approach 
(Sec. |H|) and to the general strategy for the cal- 
culation of physical quantities within this method 
(Sec. ED . 

(ii) To cover the whole range of applications over the 
last 25 years, following the seminal work of lWilsonl 
(I97 5a|) on the Kondo probl em and the work of 
iKrishna-murthv et all ( I980aj ) on the Anderson im- 
purity model (Sees. [TVl and IVf . 
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Apparently, the range of applicability of the NRG 
widened considerably, in particular over the last ten 
years. This can be easily seen in the list of references 
in which more than 50% of the entries are from the years 
starting with 2000. In physical terms, the NRG is now 
being used to study very different phenomena of con- 
densed matter physics: Typical correlation phenomena 
such as the Mott transition and heavy-fermion behavior, 
the physics of a two-state system in a dissipative envi- 
ronment, Kondo correlations in artificial atoms such as 
quantum dots, to name but a few. Of course, we ex- 
pect that there are still very many problems to which 
the NRG will be applied in the future, and we hope that 
this review will be helpful as a starting point for such 
investigations. 

Some of the concepts discussed in Sees. [TTI and Mil are 
fairly recent developments: For example the generaliza- 
tion of the NRG to quantum impurities coupled to a 
bosonic environment (see also Sec. IIV.E[) and the calcu- 
lation of time-dependent quantities (transient dynamics, 
see Sec. IIII.B.3I) . As only a few applications of these new 
concepts have been considered so far, one line of future 
research of the NRG is their generalization to a broader 
class of impurity models. 

We already discussed some open issues and ideas 
for further investigations in the various subsections of 
Sees. IIVI and [V] Let us mention here a few suggestions 
for further generalizations and applications of the NRG: 

• Application of the bosonic NRG to generalizations 
of the spin-boson models such as coupled spins in 
a dissipative environment. 

• Magnetic, orbital and charge ordering in lattice 
models within DMFT. 



• Application of multiple-shell techniques (Sec. IIII.Bl) 

to further improvement of the dynamics, particu- 
larly at the lowest temperatures. 

What are the main open issues of the NRG approach? 
As discussed in Sec. IIV.DI multi-site and multi-orbital 
models pose severe technical problems for the NRG, be- 
cause the Hilbert space increases dramatically with the 
number of orbitals. This limits, in particular, the ac- 
curacy in the calculation of dynamical quantities which 
in turn restricts the applicability of the NRG to multi- 
band models within DMFT (see Sec. IV. A[) or its exten- 
sions. Concerning dynamical quantities, another short- 
coming of the present implementations of the NRG is the 
rather poor resolution at high frequencies, for example 
features such as the band edges of upper and lower Hub- 
bard bands in the Hubbard model or the sharply peaked 
and highly asymmetrical spin-resolved Kondo resonance 
at high magnetic fields, as shown in Fig. [71 

A gradual improvement of the accuracy can, of course, 
be achieved by simply increasing the computational ef- 
fort, but for a real breakthrough (concerning multi-band 
models and the high-energy resolution) we probably need 
completely new ideas and concepts. 



From a conceptual point of view it will be very inter- 
esting to view the NRG in a bro ader context. One step i n 
this direction has been made in IVerstraete et all ( 20051 ). 
The authors of this work interpret the NRG iteration in 
terms of matrix product states, and establish a connec- 
tion to the widely used density matrix renormalization 
group method. 

Concerning the future prospects of the numerical 
renormalization group, let u s conclude wit h a remark 
from Wilson's original paper (|Wilsonl (|l975al) . page 777), 
about the prospects of the renormalization group in gen- 
eral: 

... However, most of the unsolved prob- 
lems in physics and theoretical chemistry are 
of the kind the renormalization group is in- 
tended to solve (other kinds of problems do 
not remain unsolved for long). It is likely that 
there will be a vast extension of the renor- 
malization group over the next decade as the 
methods become more clever and powerful; 
there are very few areas in either elementary 
particle physics, solid state physics, or the- 
oretical chemistry that are permanently im- 
mune to this infection. 
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